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studies  carried  out  for  the  Air  Force  by  General  Americ...  Transportation 
Corporation  and  Newmark-Hansen  Associates.  The  five  volumes  comprise 
RTD  TDR-63-30%  and  are  organized  as  follows: 

Vol.  I  Structure  Interior  Motions  Du*'  to  Air  Blast  Induced 
Ground  Shock 

Vol.  II  Structure  Interior  Motions  Due  to  Directly  Transmitted 
Ground  Shock 

Vol.  Ill  Response  Spectra  of  Single-Degree-of-Freedom  Elastic 
and  Inelastic  Systems  _ _ _  ...  . - 

Vol.  IV  Response  Spectra  of  Two-Degree-of-Freedom  Ela_  ;ic 
and  Inelastic  Systems 

Vol.  V  Response  Spectra  of  Multi-Degree-of-Freedom  Elastic 
Systems 

Volumes  I  and  II  are  authored  by  General  American  Transportation 
Corporation.  Volumes  III,  IV,  and  V  are  authored  by  Newmark-Hansen  and 
Associates.  Volumes  II,  IV,  and  V  will  be  published  early  in  1965. 

Acknowledgment  is  made  to  Captain  H,  Auld,  Captain  D.  H.  Merkle, 
and  Lt  J.  F.  Flory  of  AFWL  for  their  continued  cooperation  during  the 
course  of  the  project  and  to  Dr.  S.  Raynor  of  Northwestern  University  who 
serves  as  consultant  to  MRD  of  General  American  Transportation  Corporation. 


! 


RTD  TDR-63-3096,  Vol  I 


i  : 


*»u  d 


'O. 


ABSTRACT 

The  primary  purpose  of  this  report  is  to  provide  guidance  for  designers 
of  shock  isolation  systems  during  the  initial  phases  of  design.  Volume  1 
presents  methods  for  estimating  appropriate  free  field  waveforms  and  the 
influence  of  soil- structure  interaction  upon  interior  structure  motion.  Volume 
III  of  RTD  TDR-63-3096,  prepared  by  Newmark-Hansen  Associates,  presents 
methods  for  synthesizing  peak  relative  response  spectra  from  the  spectra 
characteristic  of  pulses  of  simple  shape. 

The  inherent  error  in  shock  isolation  design  is  at  least  +20  percent;  to 
reduce  this,  much  more  soil  test  data  than  is  now  available  will  be  required. 
Further,  based  on  purely  theoretical  arguments,  shock  isolation  per  se  can 
be  eliminated  for  much  equipment  used  in  hard  installations.  If  isolators  are 
required  they  should  be  designed  as  low  frequency  systems  that  impose  one  g 
acceleration  on  the  isolated  equipment.  Increase  of  the  acceleration  to  be 
tolerated  by  the  equipment  will,  in  the  great  majority  of  instances,  reduce 
neither  the  rattle  space  required  nor  the  isolator  cost. 

Methods  for  making  the  necessary  computations  are  given. 
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missile  have  been  established.  The  inside  diameter  of  the  silo  within  which 
the  miss? le  will  be  house '  may  depen!  more  strongly  on  the  rattle  space  required 
by  the  missile  and  its  isolation  system  than  upon  any  other  single  factor  save 
missile  envelope  diameter.  The  cost  of  the  missile  silo  and  associated  door 
and  operating  equipment  varies  somewhere  between  the  square  and  the  cube  of 
the  silo  diameter,  other  factors  being  equal.  Therefore  the  shock  isolation 
designer  is  placed  in  the  difficult  position  of  having  to  establish  space  re¬ 
quirements  that  are  adequate  for  isolation  but  held  to  a  minimum  in  order  to 
keep  costs  down  at  a  time  when  he  has  little  detailed  information  upon  which 
to  base  his  estimates.  The  primary  purpose  of  this  report  is  to  provide  methods 
that  will  assist  the  shock  isolation  designer  in  arriving  at  conclusions  to  this 
and  similar  problems. 

Two  general  avenues  of  approach  were  open  at  the  outset  of  the  work  herein 
reported: 

1.  Formulas  and  the  necessary  associated  compendium  of  data  could  be 
prepared  so  that  any  problem  likely  to  face  a  designer  would  have 
been  considered  and  graphs  presented  covering  suitable  ranges  of 
the  various  parameters. 

d.  Methods  could  be  developed  the  use  of  which  would  enable  the 

designer  to  consider  the  particular  problem  he  faced,  from  basic 
principles  and  to  expeditiously  carry  out  the  design  without 
reliance  upon  a  compendium  cf  previously  computed  data. 


Both  approaches  were  considered  on  the  project  but  it  became  evident  that 
the  latter  approach  vas  by  far  the  more  desirable.  Shock  isolation  design  for 
the  suspension  system  of  a  large  liquid  fueled  missile  and  the  design  of  isolators 
for  a  rugged ized  item  of  electronic  equipment  counted  on  a  wall  bracket  have 
little  in  common  other  than  the  name  "shock  isolation"  and  the  fundamental  laws 
of  mechanics  governing  the  design.  The  compendium  required  for  the  first  approach 
would  have  been  truly  massive.  B’urther,  it  developed  that,  within  engineering 
accuracy,  the  second  method  became  quite  practical. 

The  vague  concept  "shock  sensitivity  of  equipment"  in  the  past  has  received 
considerable  attention.  Nevertheless  at  the  present  time  it  is  practically 
impossible  to  clearly  define  the  shock  sensitivity  of  an  individual  item  of 
equipment  let  alone  to  formulate  a  general  sensitivity  criterion  with  sufficient 
accuracy  and  mathematical  precision  to  enable  deduction  of  shock  isolation 
criteria.  The  complexities  of  equipment  used  are  so  great  and  so  varied  that 
efforts  in  this  direction  have  been  unavailing.  However,  at  least  one  precise 
statement  fundamental  to  shock  sensitivity  cun  be  made:  damage  to  an  equipment 
element  is  induced  when  either  the  total  energy  communicated  to  the  element  or 
the  rate  of  flow  of  energy  delivered  to  the  element  exceeds  certain  critical 
values.  Often,  of  course,  it  is  a  practical  impossibility  to  distinguish  the 
critical  energy  component  and/or  its  critical  level.  Knowing,  however,  that  the 
damuging  agent  ultimately  is  energy,  regardless  of  how  described,  we  can  gain 
considerable  insight  into  the  overall  shock  isolation  problem  by  considering  the 
path  followed  by  the  energy  leaving  an  exploding  weapon  and  ultimately  arriving 
at  a  sensitive  element  of  equipment. 


Figure  1,1  shows  two  block  diagrams  representing,  in  a  general  way,  the 
propagation  of  energy  originating  at  the  detonation  of  a  nuclear  weapon  through 
the  blast  environment  to  shock  sensitive  equipment.  One  block  diagram  repre¬ 
sents  the  conditions  for  whicxi  we  have  da^a  obt&inec?  at  the  nuclear  weapons 
test  sites;  the  second  represents  the  situation  at  a  proposed  installation  site. 
Consider,  first,  the  test  site.  Eight  blocks  are  shown.  Of  these  we  have 
factual  information  for  three  blocks;  yield,  the  detonation  parameters,  range 
and  height  of  burst;  the  free  field  gage  outputs  and,  in  a  few  instances, 
response  gage  outputs.  We  have  partial  infomation  for  a  fourth  block,  the 
soil  characteristics  and  depth  at  which  measurements  were  taken.  It  should  be 
observed  that  the  output  ee  the  free  field  gage  is  not  what  might  be  termed 
the  true  free  field  underground  blast  motion  condition.  It  is  the  output  of 
a  gage  which,  by  its  very  presence,  modifies  the  free  field  that  would  exist  in 
its  immediate  vicinity  if  the  gag*1  were  not  there.  Further,  it  is  the  output 
measured  over  a  relatively  small  area,  perhaps  that  of  a  cylinder  a  few  inches 
in  diameter  or  less.  Portions  of  the  gage  output  would  be  substantially  the 
same  if  the  gagf  were  moved  significant  distances,  with  range,  soil  parameters 
and  depth  held  constant.  Other  portions  of  its  output,  however,  would  vary 
markedly  if  the  gage  were  moved  only  a  few  feet  or  even,  in  some  cases,  a  few 
inches.  These  latter  portions  can  be  accounted  for  in  any  design  procedure 
only  on  a  statistical  basis,  if  they  are  indeed  significant. 

At  a  weapons  system  installation,  it  should  be  noted,  the  energy  released 
by  detonation  of  the  bomb  is  modified  at  least  four  times  before  reaching  the 
portion  of  the  equipment  subject  to  damage.  The  energy  is  first  converted  from 
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the  pressure -volume  and  kinetic  energy  of  the  bomb  fragments  and  air-soil 
environment  into  the  energy  that  we  term  the  free  field  pressure.  This,  in 
turn,  may  or  maj  not  be  markedly  modified  by  any  structure  buried  in  the  free 
field.  Then  the  energy  is  further  greatly  modified  by  any  isolator  provided 
within  the  structure  between  the  structure  itself  and  the  equipment  mounted 
therein.  It  is  to  gain  this  modification  that  the  isolator  is  provided. 

Finally,  in  coursing  through  the  equipment,  that  portion  of  the  energy  passing 
through  the  isolator  is  again  further  modified  significantly  before  reaching 
the  particular  element  of  the  equipment  at  which  damage  might  be  caused.  Our 
problem  is  to  deduce  from  the  limited  information  represented  by  the  shaded 
blocks  in  the  test  site  block  diagram  a  series  of  criteria,  theories,  and  design 
formulas  which  will  enable  us,  within  an  acceptable  degree  of  accuracy,  to  follow 
the  energy  through  the  block  diagram  representing  conditions  at  an  installation 
site  and  arrive  at  meaningful  criteria  for  the  design  of  isolators  of  minimum 
size  and  cost  but  adequate  mechanical  properties  to  insure  that  the  v-quipment 
mounted  therein  survives  the  blast . 

Since  the  energy  successively  passes  through  some  four  major  modifications 
before  arriving  at  the  point  at  which  it  can  cause  damage,  it  is  evident  that 
neglect  of  one  of  the  four  can  introduce  more  error  in  our  overall  design  procedure 
than  might  be  e1 imira^e  j  by  refined  design  of  any  of  the  other  energy  conversions. 
In  particular,  the  block  diagram  shown  in  Figure  1.2  which  has  taciuly  been 
assumed  for  some  hard  structure  shock  isolation  designs,  can  be  a  poor  repre¬ 
sentation  of  the  overall  problem.  The  methods  of  extrapolating  test  site  gage 
records  to  conditions  au  an  installation  site  are  questionable  for  large  structures. 


The  point  of  departure  for  the  project  reported  herein  was  Block  3,  the 
free  field  motion  history  of  Figure  1.1.  The  energy  chain  from  the  bomb  burst 
to  free  field  is  an  extremely  complex  area  that  is  receiving  a  great  deal  of 
effort  by  a  number  of  investigators  at  the  present  time.  For  this  project  we 
have  made  use  of  information  presently  available  to  describe  generally  the 
free  -pield  motion  history  that  would  exist  under  any  particular  set  of  circum¬ 
stances.  However,  we  have  not  imposed  the  requirement  that  the  free  field  data 
be  of  any  particular  form  or  characterizable  by  any  specific  mathematical 
expression.  Tims,  as  further  information  is  developed  on  the  nature  of  the 
free  field  this  can  be  fit  into  the  shock  isolation  design  scheme  presented. 

A  substantial  portion  of  tne  effort  carried  out  on  the  project  is  concerned 
with  the  second  energy  change  represented  by  block  4  of  Figure  1.1,  the  soil 
structure  interaction.  Even  though  a  structure  buried  in  soil  remains  substantially 
rigid  during  ground  motion  its  presence  has  a  significant  effect  on  the  total  energy 
available  to  the  isolator  and  isolated  equipment  and  possibly  even  more  influence 
on  the  rate  of  delivery  of  this  energy.  Ihe  natural  vehicle  for  jar  vying  out 
these  computations  is  the  normal  mode  method  which  is  discussed  in  considerable 
detail  in  Appendix  A  of  this  report. 

In  the  following  subsections  particular  aspects  of  the  isolation  problem 
are  discussed  in  nonmathematical  terms. 


1  1  Normal  Modes,  Damping  and  Coupling 

1.1  1  Undamped  Modes 

The  vibration  of  an  undamped  linear  system,  that  is  a  system  in  which  all 
displacements  are  proportional  to  the  loads,  free  of  damping,  always  can  be 
expressed  as  an  infinite  series  of  vibrations  of  certain  special  configurations 
termed  "normal  modes".  This  is  true  no  matter  how  complex  the  overall  system  of 
vibration  is  and  is  independent  of  whether  the  vibration  is  free  or  forced; 
indeed  for  free  vibration  each  modal  motion  consists  of  a  particular  geometric 
configuration,  each  point  of  which  oscillates  harmonically  in  phase  with  all 
other  points  at  a  specific  frequency  characteristic  of  the  mode  .  Such  a  vibra¬ 
tion  would  continue  indefinitely.  The  total  ene.gy  of  each  m^-de  would  be 
independent  of  time,  and  there  would  be  no  interchange  of  energy  between  the 
various  excited  modes.  It  is  this  energy  independence  quality  that  ultimately 
defines  the  various  modal  configurations.  At  the  maximum  displacement  of  a  mode 
all  of  its  energy  would  be  potential,  stored  in  the  system  elements.  A  quarter 
period  later  all  of  the  energy  would  be  kinetic..  For  this  reason  this  energy 
Is  termed  ttr  ’dating  energy  of  the  mode. 

When  tne  mathematics  are  carried  out  the  amplitude  of  the  free  vibration  of 
a  normal  mode  is  found  to  satisfy  the  same  equation  as  a  simple  mass-spring  system, 
x.e-, 

x  +ai2x  =0  (1.1) 

n  n  n 
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where 


is  the  amplitude  of  the  nth  mode,  a  function  of  time 
O)  is  the  natural  frequency  of  the  nth  mode. 

If  a  distributed  system  of  forces  is  applied  to  the  structure,  vibration 
again  will  consist  of  the  infinite  series  of  normal  modes.  However,  in  this 
instance  each  mode  will  not  vibrate  harmonically,  rather  the  medal  amplitude 
will  be  given  by  solution  of  Equation  (1.2). 


x'  +  co  x  ■  f 
n  n  n  n 


(1.2) 


This  equation  is  seen  to  be  analogous  to  the  equation  of  forced  motion  for  a 
mass -spring  system  where  the  function  f  fills  the  role  of  forcing  function. 
It  will  be  found  that  fn  is  equal  to  the  rate  of  energy  delivery  to  the  nth 
mode  by  the  distributed  forces  applied  to  the  structure  and  multiplied  by  a 
constant  characteristic  of  the  mode. 


Since  the  energy  delivered  to  a  particular  mode  by  a  forcing  function 
generally  is  inversely  proportional  to  a  power  of  the  modal  frequency  on  the 
order  of  2-to-3  the  significance  of  the  high  frequency  modes  decreases  rapidly. 
Therefore,  though  me  thematically  accurate  solution  of  the  problem  requires 
consideration  of  a  infinite  number  of  modes,  engineering  accuracy  generally  can 
be  obtained  by  considering  two  or  three,  and  often  times  only  the  fundamental 
mode  is  of  real  significance.  An  exception  to  this  situation  occurs  when  a 
resonance  effect  takes  place;  however,  for  the  transient  response  -hat  is  of 
principal  importance  for  nuclear  blast  induced  shock,  resonance  plays  a  secondary 
role  because  the  forcing  function  (ground  shock  wave)  does  not  last  a  sufficient 
number  of  periods  to  build  up  a  strong  resonance. 
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An  aspect  of  normal  mode  theory  that  is  unfortunately  frequently  glossed 


over  in  the  literature  deserves  specific  attention.  Most  of  the  elementary 
,’oblems  discussed  in  texts  on  vibration,  e.g.,  strings,  beams,  membranes, 
and  plates  possess  a  characteristic  that  simplifies  analysis  but  is  not  intrinsic 
to  the  general  problem:  The  displacements  of  all  points  of  the  structures 
mentioned  are  parallel  to  each  other,  i.e.,  all  displacements  of  the  string  and 
the  beam  are  perpendicular  to  the  axes  of  these  members;  the  displacement  of 
all  points  on  the  surface  of  the  plate  and  membrane  are  perpendicular  to  the 
rest  planes  of  these  elements.  In  general,  this  situation  is  not  true.  Thus 
the  displacements  of  the  circular  cylinder  vibrating  in  its  bending  modes  consists 
of  a  radial  component  and  a  simultaneous  lateral  component  tangent  to  the  surface 
of  the  cylinder.  Though  the  motion  of  any  particular  point  on  the  surface  of  the 
cylinder  is  along  a  straight  line,  the  lines  at  different  points  on  the  structure 
are  not  parallel.  Further,  for  the  cylinder  there  is  no  point  at  which  the  scalar 
magnitude  of  the  motion  is  zero  in  a  manner  analogous  to  the  nodal  points  of 
vibration  of  strings  and  beams  or  the  nodal  curves  of  membranes  and  plates. 

1.1.2  Damping  and  Coupling 

It  can  be  shown  mathematically  that  the  normal  mode  solution  exists  for  a 
damped  structure  if  the  equations  of  motion  can  be  put  in  the  form  of  Sturm- 
Liouville  problem.  Within  this  limitation  the  only  case  which  is  of  physical 
significance  and  which  may  be,  from  a  practical  standpoint,  solves,  mathematically 
is  the  case  where  the  damping  force  is  proportional  to  velocity  at  every  point  of 
the  system  and  has  the  same  direction  as  the  displacement  vector  at  every  point 
of  the  system.  Furthermore,  the  scalar  proportionality  constant,  i.e.,  the 
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damping  constant  relating  the  damping  force  vector  and  the  velocity  vector, 
must  be  constant  for  the  entire  system.  Most  real  systems  do  not  comply  with 
this  criterion.  However,  it  often  can  be  used  as  a  fair  approximation  to  the 
real  system  1 it  must  be  appreciated  that  the  major  reason  for  making  the 
approxirt  hion  is  not  physical  but  mathematical;  the  problem  is  much  simpler  to 
solve  if  velocity  proportionate  damping  is  assumed. 

If  the  damping  of  the  system  is  not  velocity  proportionate  the  normal  mode 
method  of  solution  can  be  employed  though  in  this  case  it  loses  some  of  its 
physical  reality  and  becomes  essentially  a  mathematical  tool  of  some  complexity 
that  will,  however,  give  results  to  any  desired  degree  of  accuracy  if  sufficient 
computations  are  carried  out. 

The  influence  of  no.. velocity  proportionate  damping  can  have  serious  con¬ 
sequences  to  design;  therefore,  its  effects  will  be  discussed  at  some  length 
from  the  physical  standpoint  so  that  the  designer  may  gain  some  feel  as  to  what 
might  actually  take  place  in  a  particular  isolation  system.  The  crux  of  the 
matter  is  that  substantial  motions,  other  than  those  anticipated  by  an  elementary 
application  of  the  normal  mode  theory,  may  develop  after  passage  of  the  ground 
shock  wave.  These  might  require  considerable  additional  rattle  space  and  sub¬ 
stantial  modification  of  the  suspension  mechanism  over  those  required  if  the 
system  were  truly  undamped  in  the  mathematical  sense. 

The  physical  reasons  for  this  situation  can  most  easily  be  observed  by 
considering  the  simplest  of  multi-degree  of  freedom  systems,  a  pair  of  masses 
and  springs  in  series.  In  Figure  1.3  we  have  shown  two  masses  in  series  with 
two  springs  but  only  one  dashpot-  connecting  the  two  masses.  Evident!  the 
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Figure  1.3  ELEMENTARY  COUPLED  LOT  \R  SYSTEM 
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requirement  that  the  proportionality  constant  for  damping  be  the  same  for  all 
elements  is  not  met  by  the  system  shown  in  Figure  1.3.  Now  if  the  dashpot  were 
eliminated  the  system  shown  would  have  two  normal  modes.  For  each  mode, the 
contribution  to  the  displacements  x^  and  x^  would  always  be  proportional  though 
the  proportionality  constants  for  the  two  modes  would  differ.  With  the  dashpot 
in  place,  however,  during  motion  of  the  masses  a  for  a  is  developed  by  the  dashpot 
which  is  proportional  to  the  difference  between  the  velocity  of  mass  2  and  maos  1. 
The  energy  dissipated  is  equal  to  the  integral  of  the  product  of  this  dashpot 
force  and  the  differential  relative  displacement  of  the  masses.  However,  the 
dashpot  force  is  applied  to  both  masses.  Therefore,  if  we  assume  that  the  gross 
mefcion  of  the  dashpot  is  to  the  right  the  dashpot  system  will  remove  energy  from 
the  Mg,  kg  system  and  add  energy  to  the  M^,  system.  In  other  words,  in  addition 
to  dissipating  energy  as  heat,  the  damping  element  transmits  energy  from  one  spring- 
mass  element  to  the  other,  depending  on  motions  of  the  elements  connected. 

In  Appendix  A  of  this  report  the  analysis  of  the  general  system  is  presented. 
There  it  is  found  that  for  nonvelocity  proportionate  damping  the  damping  factors 
for  all  modes  couple  energv  into  all  other  modes.  Thus,  to  return  to  Figure  1.3 
if  one  mode  only  were  to  be  excited  the  dashpot  Immediately  would  begin  to  remove 
energy  from  thi.  mode  and  transmit  it  to  the  mode  that  was  not  initially  excited. 

Of  course,  while  doing  this  it  also  would  dissipate  some  energy  as  heat.  If  the 
vibration  continued  for  a  long  enough  period  of  time,  intimately  all  of  the  energy 
initially  supplied  to  the  excited  mode  would  have  been  dissipated  as  heat  or 
coupled  to  the  mode  that  was  not  initially  excited.  In  general,  this  transfer  of 
energy  requires  a  great  many  periods  of  vibration.  Therefore,  the  effect  is  most 
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pronounced  when  the  amount  of  damping  is  small  so  that  only  a  small  portion  of 
the  energy  is  dissipated  as  heat  during  each  period  of  vibration  and  the  system 
continues  to  vibrate  long  enough  for  the  damping  elements  to  transfer  the  energy 
from  uhe  excited  to  the  unexcited  mode.  These  are  precisely  the  conditions  that 
exist  in  a  real  system  that  we  normally  consider  undamped. 

The  several  modes  for  real  systems  may  differ  ma~*v.edly  from  each  other  in 
their  geometry  and  directions  of  motions.  For  instance,,  a  pendulum  supported 
on  springs  has  not  only  the  obvious  pendulum  motion  in  two  directions  and  the 
gross  vertical  and  lateral  motions  of  the  support  but  also  a  rotation  about  the 
axis  of  the  pendulum.  If  a  suspension  system,  such  as  used  for  some  missiles, 
consisting  essentially  of  a  pendulous  cage  containing  the  missile  and  supported 
by  springs  at  the  top  of  the  silo  were  not  sufficiently  damped, the  energy 
communicated  to  the  pendulum  by  the  ground  shock  wave  might,  after  a  period  of 
time,  appear  an  a  rotary  motion  of  the  entire  suspension  system  about  a  vertical 
axis,  a  situation  that  could  induce  severe  damage  to  the  missile  and  its  suspension 
system  if  not  provided  for  in  the  design.  Further  even  if  the  torsional  mode  were 
not  significantly  excited,  coupling  of  even  a  small  amount  of  the  energy  of  verti¬ 
cal  oscillation  to  the  horizontal  pendulum  mode  can  be  serious  because  a  small 
amount  of  circulating  energy  can  represent  a  wide  swiig  for  a  long  pendulum. 

The  modes  of  a  system  having  nonvelocity  proportionate  damping,  as  discussed 
in  the  previous  paragraphs,  are  accurately  described  as  coupled  modes.  The  damping 
couples  energy  from  one  mode  to  the  other.  Unfortunately,  the  term  "coupled  modes" 
is  sometimes  applied  to  an  entirely  different  motion  configuration.  This  may  be 
explained  as  follows.  The  vibration  of  a  general  system  is  a  vector  quantity,  i.e., 
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for  a  given  mode,  motion  at  different  parts  of  the  s/stem  may  take  place  in 
different  directions  even  though  the  motions  are  at  all  times  in  phase  with  each 
other.  For  any  particular  vector  mode,  this  motion  can  be  resolved  into  scalar 
components  in  any  suitably  chosen  set  of  coordinate  directions.  The  scalar 
components  are  sometimes  referred  to  as  coupled  modes.  These  cannot  exist 
independently  of  each  other  and  individually  do  not  have  the  simple  energy 
relations  characterizing  true  modes.  In  this  report  when  modal  components  are 
considered  they  are  so  identified. 

If  any  damping  at  all  is  present,  during  each  cycle  of  vibration  a  certain 
fraction  of  the  circulating  energy  is  removed  from  each  mode.  Normally,  damping 
is  considered  to  be  a  conversion  of  the  mechanical  circulating  energy  into  heat 
by  one  mechanism  or  another,  dashpots,  internal  friction,  etc.  A  more  general 
definition  of  damping  is  any  mechanism  that  removes  some  of  the  circulating  energy 
from  the  system  during  the  course  of  a  period  of  vibration.  Under  this  definition 
the  divergent  wave  of  energy  emanating  from  a  structure  vibrating  in  a  material 
medium  is  a  form  of  damping.  The  energy  carried  by  the  wave  does  not  return  to 
the  structure.  For  instance,  a  missile  silo  vibrating  in  soil  will  radiate  a 
substantial  energy  wave.  The  energy  carried  by  this  wave  does  not  return  to  the 
silo.  Therefore,  this  aspect  of  soil  structure  interaction  introduces  a  signifi¬ 
cant  amount  of  damping.  In  Appendix  A  it  is  shown  that  for  the  lower  order  modes 
of  reasonably  proportioned  structures  the  damping  is  greater  than  critical. 
Therefore,  for  these  modes  there  is  negligible  coupling  of  energy  between  modes 
induced  by  departure  of  the  real  damping  from  velocity  proportionality  and  approxl 
mations  can  be  made  for  computational  purposes  without  thereby  neglecting 
significant  motion  components. 


1.2  Shock  Spectrum  . 


The  peak  relative  response  spectrum  provides  a  useful  tool  for  shock  isola¬ 
tion  design  under  certain  circumstances.  The  companion  volume  of  this  report, 
prepared  by  Newmark-Hansen  &  Associates  is  devoted  exclusively  to  the  several 
aspects  of  shock  spectra  as  they  are  generally  termed.  Here  we  point  out  simply 
the  general  ideas  behind  the  shock  spectrum  approach  and  call  to  the  attention 
of  the  reader  how  certain  information  developed  by  use  of  the  techniques  presented 
in  this  report  can  be  used  as  a  starting  point  for  application  of  the  Newmark- 
Hansen  report. 

A  shock  spectrum  is  simply  a  plot  of  peak  displacement,  peak  velocity,  or 
peak  acceleration  of  a  mass -spring  system  against  the  system's  natural  frequency 
for  a  specific  forcing  function.  If  one  has  available  a  shock  spectrum  for  a 
forcing  function  input  of  interest  then  the  peak  amplitude  of  each  of  the  modes 
of  any  linear  vibrating  system  characterizable  by  normal  modes  and  subjected  to 
the  forcing  function  can  be  taken  fro*,  the  spectrum  if  the  frequencies  of  vibration 
of  the  modes  can  be  estimated.  Further,  if  the  damping  of  the  system  of  interest 
is  sufficiently  small  so  that  each  cycle  of  vibration  is  approximately  sinusoidal, 
the  modal  motions,  displacement,  acceleration,  and  a  so-called  pseudo  velocity, 
that  very  nearly  matches  actual  peak  velocity  can  all  be  taken  from  the  same  graph, 
if  it  is  plotted  on  quadruple  coordinate  logarithmic  paper. 

£?ince  the  spectrum  is  merely  a  plot  of  peak  values  of  response  of  mass¬ 
spring  systems  versus  frequency,  it  is  applicable  strictly  only  to  linear, 
uncoupled  systems.  Further,  if  damping  is  to  be  considered,  the  percentage  of 
critical  damping  used  in  development  of  the  spectrum  must  be  the  same  as  the 
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percentage  of  damping  of  the  equipment  to  be  designed  with  the  aid  of  the  spectrum. 
Also,  in  the  mathematical  sense  the  peak  relative  response  spectrum  is  not 
applicable  to  couple  ..ystems  or  to  nonlinear  systems.  In  actual  practice  the 
errors  introduced  by  employing  spectra  to  design  systems  for  which  they  are  not 
strictly  applicable  often  are  small  enough  to  neglect . 

In  this  report  methods  are  developed  for  determining  the  motion  on  the 
interior  of  a  protective  structure  in  terms  of  the  free-field  motion.  For  large 
structures  there  can  be  considerable  divergence  between  these.  The  Newmark- 
Hansen  report  presents  approximate,  but  quite  accurate,  methods  for  constructing 
spectra  from  certain  cardinal  characteristics  of  the  input  motion.  Thus,  by  a 
two-step  process, developing  first  the  interior  structure  motion  and  then  the 
spec  urn  it  is  possible  to  construct  a  peak  relative  response  spectrum  for  the 
design  of  shock  isolator^  reflecting  in  the  spectrum  the  soil-structure  inter¬ 
action,  structure  orientation  and  free-field  ground  shock  wave  chara~  -  ristics. 

The  shock  spect'-um  also  provides  a  useful  tool  for  evaluation  of  the  effects 
introduced  by  simplifications  of  an  analysis.  Figu  e  1.4  shows  a  comparison  of 
the  spectra  resulting  from  a  step  functior  velocity  input,  and  a  decaying  expo¬ 
nential  velocity  input  of  the  same  peak  amplitude.  A  glance  at  the  spectra  is 
all  that  is  required  to  show  that  if  the  decay  constant,  ct,  has  a  numerical 
magnitude  of  1  sec  1  or  greater,  no  more  than  10  percent  error  is  introduced  by 
a  step  wave  approximation  at  isolator  frequencies  of  5  cps  or  greater. 
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1.3  Equivalent  Drop  Method 

The  equivalent  drop  method,  is  particularly  useful  during  the  early  stages 
of  design  for  making  upper  limit  estimates  of  rattle  space  requirements  and 
establishing  for  a  considerable  variety  of  equipment  whether  or  not  shock 
isolators  will  be  necessary. 

The  soil  particle  motion  at  depth  in  real  soils,  due  to  detonation  of  a 
large  nuclear  weapon  may  be  quite  complex;  it  will  have  at  least  two  mutually 
perpendicular  components  of  motion  and  possibly  three.  However,  r'n  particulate 
soils  subjected  to  air  blast  the  vertical  component  of  motion  generally  is 
considerably  larger  than  either  of  the  horizontal  components.  For  the  remainder 
of  this  discussion,  we  will  focus  our  attention  on  the  vertical  motion  component 
only,  realizing  that  the  conclusions  drawn  apply  to  the  smaller  horizontal 
components  also. 

In  general,  when  engulfed  by  the  ground  shock  wave  the  soil  particle  vertical 
velocity  component  very  rapidly  increases  from  zero  to  some  peak  value  and  then 
much  more  slowly  decreases  back  to  zero.  In  the  subseismic  region,  generally 
below  the  100  psi  overpressure  level  for  particulate  soils,  the  vertical  parti< 
velocity  msy  oscillate  to  some  degree.  Further,  even  in  the  superseismic  -sgion 
if  bedrock  is  near  the  surface  subjected  to  air  blast,  the  ground  shock  w&v  : 
reflected  upward  from  bedrock  may  impart  a  small  oscillatory  component  to  the 
vertical  particle  velocity  though  generally  this  is  of  not  more  than  one-half 
period  duration.  The  area  under  the  particle  velocity  versus  time  curve  ^ives 
the  net  vertical  displacement  of  the  soil  particle. 
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These  complex  motions  can  be  conservatively  idealized  by  considering  the 
vertical  particle  velocity  to  be  a  step-wave.  If  it  were  not  for  the  reflection 
from  bedrock  a  surface  particle,  under  this  assumption,  would  continue  to  move 
downward  indefinitely.  However,  if  the  reflected  wave  proceeding  upward  from 
bedrock  is  assumed  to  be  equal  in  amplitude  to  the  incident  wave,  the  downward 
motion  of  the  surface  particle  would  be  stopped  at  the  instant  of  arrival  of 
the  reflected  wave,  though  the  permanent  displacement  of  the  surface  thereby 
obtained  would  be  considerably  greater  than  the  permanent  displacement  actually 
resulting  from  detonation  of  a  real  bomb.  Evidently,  this  idealization  may 
considerably  over-estimate  soil  particle  motion.  However,  a  step  function  of 
particle  veloci+v,  or  velocity  jump  as  it  is  generally  termed,  enables  very  rapid 
estimates  of  rattle  space  and  shock  input  to  equipment. 

The  peak  relative  response  spectrum  of  a  velocity  jump  when  drawn  on  the 
conventional  four  coordinate  logarithmic  paper  is  just  a  horizontal  line  at 
pseudo-velocity  equal  to  the  velocity  jump  (see  Figure  lA).  Now,  if  an  object 
were  suspended  a  certain  distance  above  an  absolutely  unyielding  floor  and  the 
bottom  surface  of  the  object  parallel  to  the  floor  were  made  of  absolutely  un¬ 
yielding  material,  and  the  object  then  allowed  to  drop  and  impact  against  the 
floor  the  shock  delivered  to  the  object  would  be  equivalent  to  a  velocity  jump 
equal,  in  magnitude  to  the  velocity  obtained  by  tne  object  falling  under  gravity 
just  at  the  instant  of  impact.  The  shock  spectrum  of  this  impact  then  is  an 
upper  envelope  of  all  shock  spectra  having  peak  values  on  the  velocity  axis 
equal  to  the  velocity  jump.  Thus,  from  the  standpoint  of  shock  damage  to  equip¬ 
ment  it  is  possible  to  conservatively  equate  overpressure  levels,  bedrock  depth, 
and  soil  parameters  to  an  equivalent  drop  of  the  equipment  to  a  hard  surface. 
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(The  horizontal  component  of  shock  can  be  similarly  handled  except  that  the 
equipment  should  be  envisioned  as  turned  sideways  when  dropped.)  Though  the 
method  is  obviously  crude  it  has  the  important  advantage  that  the  somewhat 
intangible  quantities  significant  to  more  sophisticated  shock  isolation  analysis 
are  converted  into  a  parameter  for  which  engineers  have  considerable  intuitive 
feel.  Further,  the  numerical  values  of  the  equivalent  drop  and  velocity  jump 
are  surprisingly  small.  Thus,  at  the  100  psi  overpressure  level,  in  soft  soil 
having  a  modulus  of  elasticity  of  only  20,000  psi  at  the  100  psi  level,  the 
equivalent  drop  is  only  1  inches.  It  could  be  ''included  therefore,  that  most 
machinery  and  ruggedized  electronic  equipment  wound  require  no  shock  isolation 
for  these  conditions.  Only  very  sensitive  elements,  e.g.,  a  liquid  fueled  missile, 
would  be  unable  to  resist  a  4  inch  drop  under  gravity  to  a  hard  surface. 

Even  for  a  sensitive  missile  isolated  by  a  spring  system,  the  equivalent 
drop  provides  a  useful  index  for  approximate  determination  of  the  relations  between 
rattle  space,  suspension  frequency,  and  acceleration. 

The  development  of  the  formulas  for  the  equivalent  drop  approach,  including 
a  relaxation  of  the  step  function  form  of  the  parti c] e  velocity  to  forms  mor 
nearly  duplicating  actual  conditions  are  contained  in  Reference  1. 

Graphs  of  the  equivalent  drop  versus  input  parameters  are  included  in 


Section  7. 


SECTION  2 


CONCLUSIONS  AND  RECOMMENDATIONS 


2.1  Conclusions 

2.1.1  At  the  present  state  of  knowledge  highly  refined  shock  isolation 
design  is  not  justified. 

2.1.2  The  inherent  error,  due  to  soil  parameter  uncertainties,  in 
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determination  of  rattle  space  for  a  shock  isolation  system  is  at  least  -  20 

+ 

percent  and  in  some  instances  may  he  as  great  as  -  50  percent. 

2.1.3  The  largest  source  of  inherent  error  in  shock  isolation  calculations 
is  due  to  uncertainty  in  the  stress-strain  behavior  of  real  soils  at  the  over¬ 
pressure  levels  of  interest.  Uncertainties  in  the  values  of  the  propagation 
velocities  of  compression  waves,  which  ultimately  hinge  on  the  stress-strain 
behavior  of  soil,  contribute  tue  greatest  inherent  errors  to  computed  rattl 
space  requirements. 

2.1.4  Theoretical  deductions  indicate  that  at  "typical"  real  sites  the 
second  largest  source  of  inherent  error  results  from  uncertainty  (or  neglect) 
of  bedrock  elevation.  (For  some  conditions  this  is  the  largest  error.) 

2.1.5  Theoretical  deductions  indicate  that  field  tests  carried  out  at 
Frenchman  Flat  and  on  the  Pacific  i  oving  Ground  cannot  provide  significant 
information  or  the  influence  of  bedrock  at  more  typical  sites.  This  is  due 
to  the  peculiar  geology  at  these  test  sites. 


2.1.6  The  third  most  important  source  of  inherent  error  is  due  to  un¬ 
certainties  in  energy  dissipation  (wave  attenuation)  characteristics  of  soil. 

2.1.7  A  ham  structure  functions  as  an  effective  filter  of  high  fi equency 
ground  shock  components  regardless  of  assumed  soil  characteristics.  The  maximum 
frequency  that  can  he  transmitted  to  the  interior  is  the  larger  of: 

1.  The  rigid  body  frequency  in  soil. 

2.  The  first  deformations!  mode  frequency  in  soil. 

(Kote:  The  second  of  these  may  he  much  higher  than  the  first 
deformations!  frequency  in  vacuum. ) 

2.1.8  To  within  less  error  than  the  inherent  error  of  the  free  field 
motion,  the  interior  struct  tire  motion  can  be  considered  to  he  the  sum  of  two 
components : 

1.  A  rigid  body  motion  equal  to  the  free  field  parti cle  motion 
having  a  rise  modified  to  reflect  soil-structure  interaction; 

2.  A  deformational  motion  quasi -static ally  computed  but  having 
a  rise  modified  to  reflect  soil-structure  interaction;  the 
quasi-static,  deformation  is  equivalent  to  the  deformation  the 
struct' -re  would  undeigo  at  any  instant  if  the  ground  shock 
wave  were  statically  applied  to  the  soil- structure  complex. 

This  motion  is  a  function  of  time  hut,  since  it  is  computed 
quasi -static ally,  it  is  independent  of  soil  and  structure 


inertia. 


2.1.9  Theoretical  considerations,  not  known  to  have  been  verified  experi¬ 


mentally,  indicate  that  shock  isolators  can  be  omitted  for  much  equipment  in 
hard  shelters,  even  at  high  (1000  psi)  overpressure  levels. 

2.1.10  If  shock  isolation  is  required  the  minimum  energy  (and  generally 
minimum  cost)  vertical  system  will  result  if  the  frequency  is  chosen  so  that 
one  gravity  unit  acceleration  is  applied  by  the  system  to  the  isolated  equipment 
at  maximum  departure  of  the  equipment  from  equilibrium. 

2.1.11  Tnc  reduction  in  rattle  space  achieved  by  manipulation  of  suspension 
frequency  and  increase  of  tolerable  equipment  acceleration  is  only  a-  fraction  of 
the  error  inherent  in  determination  of  rattle  space. 

2.1.12  It  is  believed  that  negative  spring  systems  having  natural  frequencies 
two  or  three  orders  of  magnitude  lower  than  conventional  (~  1  cps)  systems  can  be 
developed.  A  negative  spring  is  defined  as  any  device  naving  a  negative  force 
versus  displacement  characteristic.  The  main  problem  to  be  overcome  in  systems 
using  negative  springs  is  adjustability  under  varying  magnitude  and  distribution 
of  static  load  to  maintain  constant  static  deflection.  Such  systems  would  be 
substantially  independent  of  cross  coupling  effects  between  modes,  would  not  ring 
significantly  after  excitation,  and  would  have  only  slightly  greater  energy 
capacity,  for  given  stroke,  than  linear  spring  systems. 

2.1.13  Air  springs  arid  similar  systems  having  concave  upward  force  versus 
displacement  curves  at  a  given  stroke  will  either 

Have  much  larger  energy  capacity  than  a  linear  spring  system  of 
the  same  stroke  and  peak  acceleration,  or 
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2.  Subject  the  equipment  to  a  considerably  greater  acceleration 
than  a  linear  spring  system  of  the  same  stroke  and  energy 
capacity. 

2.1.14  Theoretical  considerations  indicate  that  sharp,  short  pulses  of  stress 
having  rise  times  measured  in  microseconds  can  be  generated  within  structure  and 
equipment  by  tension  and  co’  ression  waves  propagating  through  materials  having 
concave  stress  strain  curves.  These  pulses  can  damage  small,  delicate  elements 
such  as  vacuum  tube  heaters  and  filaments.  They  can  be  prevented  from  reaching 
equipment  by  mounting  the  equipment  on  plates  or  brackets  in  such  a  way  that  all 
energy  propagated  through  the  bracket  must  appear  at  some  point  as  a  bending  wave 
(as  opposed  to  a  direct  tension  or  compression  wave). 

2.1.15  Peak  relative  response  spectra  for  interior  structure  motion  can  be 
synthesized  approximately  by  combination  of  spectra  due  to  simple  pulses.  The 
errors  introduced  by  She  approximations  are  considerably  smaller  than  the  errors 
inherent  in  the  input  (ground  motion)  data. 

2.1.16  Shear  waves  can  be  neglected  in  shock  isolation  design  without 
introduction  of  significant  error.  Neglect  of  the  shear  waves  is  slightly 
conservative . 

2.2  Recommendations  to  Shock  Isolation  Designers 

2.2.1  The  input  motion  to  the  isolator  will  be  the  least  accurately  known 
data  used  in  the  design.  Accordingly,  this  aspect  of  the  problem  deserves  as 
much  meaningful  refinement  as  can  be  carried  out.  (See  Section  4) 
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2.2.2  If  equipment  must  be  isolated  to  a  few  g  then  the  minimum  energy, 

and  probably  minimum  cost,  vertical  system  will  result  if  the  frequency  is  chosen 
so  that  -  one  g  acceleration  is  developed  at  maximum  stroke.  For  these  and  other 
low  frequency  systems  peak  displacement  only  need  be  computed.  (See  Section  4) 

Tiie  concept  of  the  "shock  response  spectrum"  is  a  useful  one  in  the  design 
of  shock  isolation  equipment,  cor  a  detailed  discussion  of  its  applicability  as 
an  index  of  equipment  shock  damage  the  reader  is  referred  to  Reference  7. 

2.2.3  A  small  component  of  circulating  energy  represents  a  wide  swing  for 
a  long  pendulum.  The  circulating  energy  of  a  vertical  system  will  be  equal  to 

the  product  of  maximum  stroke  and  equipment  weight.  Thus,  if  even  a  small  fraction 
of  this  is  coupled  to  the  pendulum  the  latter  can  be  excited  to  destructive  ampli¬ 
tudes.  Therefore,  both  systems  should  be  damped  using  Coulomb  (not  viscous) 
dampers  and  the  ratio  of  vertical  to  pendulum  frequency  should  be  as  large  as 
practical,  greater  than  three  and  an  irrational  number.  Even  frequency  multiples 
should  be  avoided,  particularly  2.  An  accurate  nonlinear  analysis  of  the  resulting 
system  should  be  carried  out. 

2.2.4  In  many  installations  isolators  can  be  omitted  from  much  equipment. 
Section  7  gives  a  simple  method  for  determination  of  a  convenient  conservative 
criterion  (equivalent  drop  distance)  for  determination  of  whether  isolators  are 
necessary. 

2.3  Research  Recommendations 

Though  the  complex  problems  of  wave  propagation  and  soil-structure  inter¬ 
action  are  as  yet  not  understood  completely,  it  is  evident  that  existing  approximate 
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techniques  introduce  much  less  error  nito  shock  isolation  calculations  than  does 
the  present  dirth  of  knowledge  of  the  stress-strain  beha.ior  of  real  soils  at  high 
pressure  levels.  Even  static  test  data  at  500  to  1000  psi  is  extremely  scarce.  I 

j 

Conceding  at  the  outset  that  dynamic  effects  may  modify  static  data  somewhat  • 
and  further  uhat  the  almost  infinite  variety  of  soils  will  preclude  ever  obtaining 

X 

a  compendium  of  precise  data  it  still  can  be  argued  that  practically  any  factual 
data  would  place  the  designer  in  a  better  position  than  he  now  occupies. 

If  a  designer  had  available  a  series  of  stress-strain  r^rves  with  numbers 
for  real  soils  he  c'-uld  at  least  use  some  judgement  in  estimating  a  curve  for 
other  soils  and  have  a  degree  of  confidence  in  the  values  of  the  slope  of  his 
approximation  at  pressures  varying  between  zero  and  peak  overpressure . 

Accordingly,  it  is  strongly  recommended  that  the  following  data  be  obtained 
from  static  and  laboratory  dynamic  tests  for  a  variety  of  real  soils: 

2.3.1  Stress  versus  strain  curves  from  confined  compression  tests  in 
the  range  between  zero  and  500  psi. 

2.3.2  Slopes  of  the  hysteresis  loops  across  the  entire  range  of  pressures. 

2.3"3  Percyclic  (or  other)  damping  factors  for  single  high  pressure  pulses. 


SECTION  3 


INHERENT  ACCURACY 
OPTIMUM  SYSTEM  CHARACTERISTICS 


If  all  of  the  input  parameters  of  a  shock  isolation  system  design  were 
known  with  three  figure  accuracy,  sophisticated  analysis  and  optimization  could 
he  carried  out  and  nicely  engineered  suspension  systems  d rve_  ped.  Regardless 
of  the  accuracy  that  might  he  desirable  however,  the  inherent  uncertainties  in 
input  parameter  values  limit  the  accuracy  that  can  he  obtained. 

In  this  section  the  problems  of  input  error  and  system  optimization  will 
he  examined.  These  will  he  used  to  establish  limiting  accur'  y  to  he  striven 
for  and  to  delineate  the  mathematical  complexity  necessary  to  achieve  this, 
accuracy. 

3.1  Wave  Reflection  in  Real  Soils 

At  most  locations  within  the  hounds  of  continental  United  States,  bedrock 
occurs  within  one  or  two  hundred  feet  of  the  surface,  though  actual  depths  at 
specific  locations  may  'retry  from  zero  at  Limestone,  Maine  to  several  hundred 
feet  in  parts  of  North  Dakota. 

If  bedrock  lies  below  the  surface  a  distance  on  the  order  of  one  hundred 
feet,  theory  predicts  that  a  strong  reflected  wave  will  he  developed  in  the 
particulate  overburden  and  that  neglect  or  consideration  of  the  reflected  wave 
introduces  an  uncertainty  factor  greater  at  the  lower  overpressures  than  that 
of  any  other  parameter  even  including  the  square  of  the  wave  propagation  velocity. 
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Considerable  evidence  has  been  presented  however  that  the  predictions  of 
elastic  theory  do  not  appear  to  occur  in  real  soils.  Sauer  (Reference  9)  points 
out  that  experimental  data  obtained  from  nuclear  field  tests  do  not  exhibit  1 

significant  reflections  of  stress  waves.  Newmark  and  Haltivanger  (Reference  6)  ( 
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state  further  that  laboratory  tests  on  soil  with  differences  in  seismic  velocity 
and  density  of  adjacent  layers  have  yielded  similar  results.  The  authors  of  ; 

Reference  7  have  concluded,,  based  on  the  data  presented  by  Musk  at  and  Meres  ' 

(Reference  1 ) ,  that  the  energy  carried  by  reflected  waves  is  relatively  unimportant, 
on  the  order  of  11  to  15  percent  of  the  incident  wave  energy.  Thus  in  spite  of 
the  predictions  of  elementary  elastic  theory  the  tendency  in  shock  isolation 
design  has  been  to  ignore  reflections  as  being  a  theoretical  complication  of 
considerrble  magnitude  yet  one  that  influences  results  well  below  the  inherent 
error  in  the  problem  and  one  that  has  not  been  evidenced  in  actual  test  data. 

In  spite  of  the  evidence  the  author  of  this  nort  does  not  agree  with  this 
philosophy  and  has  not  adopted  it  in  preparation  of  the  report. 

In  Nevada  ground  motion  data  have  been  obtained  for  only  KT  weapons.  At  the 
overpressure  levels  cl  interest  (>  100  psi)  the  positive  phase  duration  is  on  the 
order  of  l/lO  second  and  the  compression  wave  propagation  velocity  is  about  1500 

ft  per  second.  Thus,  neglecting  all  attenuation  factors  the  wave  length  of  the 
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entire  positive  phase  is  on  the  order  of  150  feet. 

At  the  Nevada  Test  Site  bedrock  occurs  at  a  depth  of  650  feet.  The  bedrock 
seismic  velocity  is  about  10,000  ft.  per  second.  Above  the  bedrock  is  an 
approximately  45 0  foot  layer  of  sand  and  gravel  having  an  average  seismic  velocity 
on  the  order  of  3000  to  3500  ft.  per  second]  above  this  is  a  200  foot  layer  of 


loess  having  an  average  seismic  velocity  of  about  2000  feet  per  second.  Since 
the  ratio  of  acoustic  impedances  of  the  regions  above  and  below  the  minus  200 
ft.  elevation  is  less  than  two,  a  sizeable  reflection  at  the  interface,  even  if 
well  defined,  is  not  expected.  Thus  the  only  sizeable  reflection  that  could  be 
expected  in  Nevada  is  that  from  bedrock.  But  the  bedrock  reflection  must  travel 
a  distance  of  about  9  or  10  initial  wavelengths.  Thus  even  a  conservative  esti¬ 
mate  of  attenuation  indicates  that  the  reflected  wave  would  have  an  amplitude 
less  than  2  percent  of  the  initial  amplitude. 

It  is  not  surprising  that  no  bedrock  reflected  war2  is  observed  in  Nevada; 
none  should  be  expected. 

In  the  Pacific  the  surface  layer,  not  more  than  about  20  feet  thick  has  a 
seismic  velocity  varying  from  800  to  4000  ft.  per  second  at  water  level.  Below 
the  water  level  a  layer  about  2500  feet  thick  has  a  average  seismic  velocity  of 
about  800  ft.  per  second.  Then  below  this  is  a  layer  5  to  10  thousand  feet  thick 
having  a  seismic  velocity  of  about  11,000  ft.  per  second. 

At  the  100  psi  level  the  positive  phase  duration  of  a  multi-megaton  bomb 
is  on  the  order  of  2  seconds.  At  higher  overpressures  the  positive  phase 
duration  increases  but  the  time  from  peak  overpiessure  to  10  percent  of  peak 
overpressure  decrease  .  For  these  large  durations  the  20  foot  surface  layer 
is  just  a  skin  separating  the  aero-shock  from  the  ground  shock  wave  in  the 
800  ft/second  layer.  Thus  the  wave  length  of  the  entire  positive  phas  ~  i  the 
2500  ft.  thick  "surface"  layer  is  about  1600  feet;  the  initial  wave  length  of 
the  portion  of  a  500  psi  wave  between  500  and  50  psi  would  be  only  about  500  ft. 
for  a  multi-megaton  bomb. 


Thus  the  reflected  waves  must  travel  a  distance  between  three  and  ten 
initial  wave  lengths  before  appearing  near  the  surface.  The  low  seismic  velocity 
exhibited  is  due  primarily  to  the  structured,  porous  nature  of  the  coral  island, 
and  therefore  the  inelastic  damping  of  the  material  is  very  high.  If  geometric 
attenuation  and  a  percyclic  damping  factor  of  0.5  are  assumed,  even  at  three 
initial  wave  lengths  propagation  the  reflected  wave  will  be  attenuated  to  about 
10  percent  of  its  initial  peak  value  and  wi1!  arrive  at  the  surface  at  least  six 
seconds  after  blast  wave  arrival. 

We  conclude  therefore  that  at  the  Pacific  Proving  Ground  evidence  of  "bedrock” 
reflecJ  'd  waves  should  not  be  expected,  particularly  at  the  higher  overpressures. 

The  Muskat,  Meres  paper  uannot  be  taken  as  evidence  that  significant  bedrock 
reflections  will  not  occur ;  these  authors  did  .xot  consider  that  problem.  Muskat 
and  Meres  have  limited  their  data  to  an  interface  acoustic  impedance  ratio  of 
2.6  ('y  a  in  their  paper).  The  elementary  one  dimensional  elastic  theory  indicates 
a  reflection  having  only  20  percent  of  the  incident  wave  energy  at  this  ratio. 

At  a  soil  rock  interface  occurring  within  one  to  two  hun.  '  "eet  of  the  surface 
the  ratio  would  ‘be  more  nearly  seven  to  ten.  Ac  an  acoustic  impedance  r;  ,  of 

7  the  reflected  energy  is  55  percent  of  the  inc.de  u:srgy.  Further,  the 
reflected  pressure  is  75  percent  of  the  incident  pressure.  Surface  displacement 
is  proportional  to  the  pressure. 

Another  argument  against  reflections  in  soil  that  is  sometimes  advanced  is 
that  the  interface  between  particulate  soil  and  rock  generally  is  not  sharp  but 
consists  of  a  five  or  ten  foot  layer  of  irregular  rock  surface,  boulders,  ce  ented 
soil,  etc.  The  transition  region  presents  formidable  difficulties  to  analysis 
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so  it  may  safely  have  imputed  to  it  any  set  of  mystical  properties  desired.  The 
reflected  wave  above  the  transition  region  consists  of  the  sum  of  an  infinite 
set  of  multiple  reflections  and  refractions  taking  place  within  the  transition. 

Here  we  will  make  no  attempt  to  consider  the  problem  mathematically  but 
will  present  a  simple  physical  discussion  that  clearly  illustrates  that  the 
physical  effect  we  generally  choose  to  characterize  as  a  mathematical  reflection 
must  occur  in  real  soils  and  will  have  significant  influence  on  surface  motions. 

Figure  3.1  shows  schematically  a  column  of  real  soil  resting  upon  a 
relatively  rigid  bedrock.  Between  the  two  is  a  relatively  thin  but  finite 
transition  region. 

If  a  pressure  is  applied  to  the  surface  of  the  sell  and  maintained,  a  wave 
is  propagated  into  the  soil  column.  Physically,  the  wave  represents  the  distance 
within  the  soil  column  that  is  under  compression,  though  the  latter  need  not  he 
constant  within  the  wave.  The  compression  causes  the  column  to  shorten  and  the 
surface  pressure  delivers  energy  to  the  soil  equivalent  to  the  product  of  surface 
pressure  and  surface  displacement.  The  compression  of  any  element  of  soil  causes 
all  elements  above  to  move  downward.  Since  the  process  is  continuous  all  elements 
of  soil  acquire  a  velocity  (particle  velocity)  which  need  not  be  equal  for  the 
various  elements.  The  particle  velocity  results  in  the  soil  column  acquiring 
kinetic  energy. 

In  most  real  soils  compression  of  the  soil  induces  some  heat.  The  sum  of 
heat  energy,  kinetic  energy  and  strain  energy  in  the  column  at  any  instant  is 
just  equal  to  the  energy  delivered  by  the  surface  pressure  acting  through  the 
distance  the  surface  has  displaced  at  that  instant. 
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When  the  wave  front  passes  through  the  transition  region  into  the  bedrock 
a  number  of  things  can  happen.  Here  we  are  concerned  only  with  the  fact  that 
some  wave  does  propagate  into  the  rock.  The  pressure  associated  with  this  wave 
can  be  either  higher  or  lower  than  the  applied  surface  pressure.  However,  we 
have  postulated  that  the  rock  be  stiffer  than  the  soil;  therefore  at  a  given 
unit  strain  a  unit  volume  of  rock  contains  more  strain  energy  than  an  equivalent 
volume  of  soil.  •Thus  if  a  unit  volume  of  rock  compressed  as  much  as  the  soil 
above  had  been  compressed,  energy  would  be  created,  an  impossibility.  We 
conclude  therefore  that  regardless  of  the  detailed  characteristics  of  soil  and 
rock  if  the  latter  is  stiffer  than  the  soil  the  unit  compression  will  be  less 
than  in  the  soil.  This  argument  is  independent  of  the  rock-soil  transition 
characteristics . 

Now  if  the  unit  compression  in  the  rock  is  less  than  that  in  the  soil  the 
veloc  by  gained  by  the  rock  is  less  than  that  already  established  in  the  soil. 

The  result  is  that  the  soil  piles  up  against  the  rock,  the  pileup  commencing  at 
the  interface  and  proceeding  to  the  surface.  (This  is  the  reflected  wave). 

During  the  pileup  the  soil  particle  velocity  is  reduced  to  the  rock  surface 
(particle)  velocity  and  the  difference  in  kinetic  energy  is  converted  into  soil 
strain  energy  and  heat.  The  additional  strain  in  the  soil  induced  by  the  pileup 
causes  the  surface  to  further  displace  and  allows  the  applied  pressure  to  deliver 
moi.3  energy  to  the  column.  However,  the  stress  induced  in  the  soil  is  greater 
than  the  applied  pressure  so  when  the  pileup  reaches  the  surface  it  immediately 
moves  upward  the  surface  on  which  the  pressure  is  applied,  causing  the  soil  column 
to  expand  and  the  surface  to  receive  an  upward  motion  component. 
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At  this  point  the  surface  pressure  is  insufficient  to  maintain  the  compres¬ 
sion  in  the  pile  so  the  latter  relaxes  starting  at  the  suiface  and  proceeding 
downward.  This  is  the  surface  reflected  tension  wave. 

Now,  it  should  he  evident  that  a  finite  transition  region  between  the  rock 
and  soil  can  affect  only  the  initial  details  of  the  process  but  not  the  grcjs 
results . 

If  the  rock  were  so  stiff  that  it  could  be  considered  with  tolerable  error 
to  be  infinitely  rigid  then  nearly  all  of  the  kinetic  energy  in  the  soil  would 
be  converted  into  strain  energy  and  heat.  The  "nearly"  qualification  is  used 
because  for  the  most  real  materials  all  of  the  particles  in  the  soil  would  not 
^simultaneously  be  reduced  to  zero  velocity.  Evidently  however,  from  physical 
reasoning,  the  residual  velocities  remaining  in  particles  after  engulfment  by 
the  reflected  wave  must  be  small.  Thus  for  a  first  estimate,  the  surface 
displacement  induced  by  the  reflected  wave  would  be  equal  to  the  displacement 
necessary  to  accommodate  as  strain  energy  the  algebraic  sum  of  three  energy 
components: 

1.  The  kinetic  energy  left  in  the  soil  by  passage  of  the 
initial  compression  wave  (positive), 

2.  The  energy  delivered  to  the  soil  by  the  surface  pressure 
acting  through  the  additional  surface  displacement  induced 
by  the  reflected  wave  (positive), 

3.  The  energy  lost  as  heat  (negative). 

The  real  blast  wave  decays  with  time  and  consequently  detailed  analysis 
becomes  quite  complex.  However,  we  crn  make  a  fair  estimate  of  the  additional 
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surface  displacement  by  considering  the  stress  strain  curve  of  the  soil.  In 
order  to  simplify  the  discussion  a  nonlinear  elastic  soil  will  be  considered, 
i.e.,  one  that  does  not  generate  heat  during  compression.  (Such  a  soil  probably 
does  not  exist.)  Since  the  heat  generating  process  (hysteresis)  produces 
permanent  displacement  this  restriction  will  not  impose  unrealistic  results. 

Figure  3.2  schematically  shows  the  stress  strain  curve  of  a  hypothetical 
nonlinear  elastic  soil.  If  pressure  and  velocity  variations  are  not  too  great 
the  same  curve  can  be  taken  as  a  reasonable  approximation  of  the  average  stress 
versus  average  strain  in  a  finite  column  of  soil. 

Now,  consider  a  column  of  soil  loaded  at  a  constant  pressure,  p.  The 

displacement  of  the  surface  is  proportional  bo  the  soil  strain.  Thus  in  - - 

Figure  3.3  the  abscissa  noted  is  proportional  to  the  initial  surface  displacement 
for  a  soil  column  of  fixed  length. 

The  energy  delivered  to  the  column  by  the  surface  pressure  is  equal  to  the 
product  of  pressure  and  displacement.  This  is  proportional  to  the  area  of 
rectangle  OADE  in  Figure  3.3.  Now  the  strain  energy  in  the  soil  is  proportional 
to  the  area  under  the  curve  ODA.  Therefore,  since  there  is  no  heat  loss  in  the 
nonlinear  elastic  soil  the  kinetic  energy  in  the  column  is  proportional  to  the 
shaded  area  ODE.  We  emphasize  again  that  these  relations  would  be  strictly  true 
only  if  the  pressure  pulse  traveling  in  the  soil  were  a  step  wave  having  a  clean 
shock  front.  Then  the  compression  would  be  constant  within  the  loaded  column  of 
soil  and  the  relations  described  here  would  be  valid.  In  fact,  however,  a  step 
pulse  is  not  conserved  in  a  nonlinear  elastic  soil.  Therefore  these  relations 
are  at  best  reasonable  approximations  to  the  true  conditions  within  the  soil. 
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Figure  3.2  STRESS-STRAIN  CURVE  OF 
HYPOTHETICAL  NONLINEAR  ELASTIC  SOIL 


Figure  3.3  SURFACE  DISPLACEMENT  INDUCED  BY 
CONVERSION  OF  SOIL  KINETIC  ENERGY  TO  POTENTIAL  ENERGY 
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Upon  reflection,  the  kinetic  energy  is  converted  into  strain  energy  causing 

the  surface  to  displace  an  additional  amount  A  B,  and  the  total  energy  delivered 

to  the  column  by  surface  pressure  is  then  proportional  to  the  area  of  rectangle 

OBCE.  If  the  velocities  of  all  particles  simultaneously  become  zero  then  this 

energy  would  have  to  be  equivalent  to  the  strain  energy  proportional  to  the  area 

ODFCBA.  Thus,  for  these  conditions  the  shaded  area  ODE  would  have  to  be  equal 

to  the  stuppled  area  DCF,  and  the  ratio  of  reflection  induced  surface  displacement 

AB 

to  the  primary  wave  induced  surface  displacement  would  be  equal  to  the  ratio 

For  a  linear  elastic  medium  this  ratio  is  just  unity.  Evidently  however, 
for  the  hypothetical  soil  illustrated  in  Figure  3.3  this  ratio  can  be  considerably 
greater  than  •unity.  For  other  stress  strain  curves  it  might  be  less.  Computations 
are  further  complicated  in  a  real  case  because  some  heat  may  be  developed,  a 
portion  of  this  heat  being  reversible  and  the  remainder  irreversible.  Neglect  of 
the  heat,  however,  is  conservative  for  design  because  its  neglect  implys  greater 
kinetic  and  strain  energy  than  actually  occur. 

There  is  no  doubt  that  when  a  decaying  surface  pressure  wave,  hysteretic 
medium  and  geometric  attenuation  are  considered,  the  ratio  of  reflected  wave 
induced  to  primary  wave  induced  surface  motion  can  vary  widely  from  the  over¬ 
simplified  conditions  of  Figure  3.3.  Our  point  is  that  whatever  the  "exact" 
value  of  the  ratio,  neglect  of  this  factor  in  many  real  instances  would  introduce 
an  error  or  uncertainty  greater  than  the  uncertainty  in  any  other  parameter 
influencing  the  ground  motion  problem. 

3.2  Suspension  Frequency  Range 

Consider  the  principal  characteristic  of  the  vertical  motion  of  a  particle 
in  the  free  field.  The  particle  is  motionless  until  the  wave  front  of  the  ground 
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shock  wave  arrives.  Then,  due  to  the  rise  time  of  the  wave  the  particle 
experiences  a  rapid  increase  in  velocity  to  a  peak  value.  The  time  required 
for  this  to  occur  is  on  the  order  of  a  few  milliseconds.  Subsequently  the  wave 
continues  to  bedrock  and  reflects,  the  reflected  wave  ultimately  engulfing  the 
particle.  If  the  wave  had  a  flat  top  and  reflection  at  bedrock  approached  that 
of  perfect  reflection  from  an  infinitely  rigid  surface,  the  velocity  of  the 
particle  would  remain  constant  during  this  interval.  1‘hen  at  engulfment  by  the 
reflected  wave  the  velocity  of  the  particle  would  decrease  to  zero. 

because  of  the  departure  of  real  soils  from  the  idealized  elastic  assumptions 
the  imperfect  reflection  at  bedrock,  the  decay  of  the  wave  subsequent  to  peak 
-  pressure,  dissipation  the  wave  -energy  as  it  travels,  and,  the,  fact,  that  not  one 

wave  but  a  series  of  waves  are  actually  involved,  the  general  character  of  the 
particle  velocity  history  is  somewhat  modified.  But  its.  overall  character  can 
be  represented  by 

1.  A  rapid  rise  of  a  few  milliseconds. 

2.  A  nearly  constant  plateau  followed  by  a  slow  decay,  the 
total  occupying  a  time  period  on  the  order  of  l/4  second 
for  large  bombs,  and  deep  (100-200  feet)  bedrock. 

3.  A  rapid  decay  through  zero  to  a  small  negative  value. 

4.  A  decaying  negative  velocity  of  long  duration. 

The  directly  transmitted  ground  shock  component  is  superimposed  on  this. 

It  should  be  appreciated  that  the  motion  of  the  selected  particle  is 
physically  due  to  the  compression  of  the  soil  between  the  particle  and  bedrock 
plus  any  bedrock  motion.  Thus  since  a  uniform  wave  advancing  at  nearly  constant 


velocity  compresses  the  soil  at  a  uniform  rate  the  velocity  of  any  selected 
particle  above  such  a  wave  would  be  nearly  constant. 

From  the  above  we  can  see  that  the  simplest  approximation  to  free  field 
particle  velocity  is  a  velocity  step  function.  This  can  be  made  more  realistic 
by  imposing  a  short  rise  time,  a  decay  rate  dependent  on  bomb  size  and  possibly 
soil  characteristics,  and  finally  terminati  the  velocity  pulBe  at  the  time 
required  for  the  ground  shock  wave  to  make  the  round  trip  from  particle  to  bedrock 
and  back  to  particle.  (These  assumptions  would  imply  permanent  displacement  of 
the  surface . ) 

If  a  structure  of  finite  size  were  immersed  in  the  soil  the  motions  observed 

at  a  point  on  the  structure  interior  would  differ  from  .  .ose  in  the  free  field. _ 

A  reasonably  proportioned  structure  (even  very  unreasonably  heavy  structures)  will 
be  lighter  than  the  soil  displaced}  therefore  it  may  be  regarded  as  a  captive 
bubble  in  the  soil.  When  engulfed  by  the  ground  shock  wave  the  bubble  deforms 
and  reflected  waves  are  radiated  from  its  surface.  Depending  on  structure  rigidity 
as  compared  with  soil  rigidity  the  reflected  waves  may  either  increase,  or  decrease 
the  interface  pressure. 

Now  it  should  be  evident  that  at  some  time  subsequent  to  engulfment  the 
captive  bubble  will  be"riding"  with  the  surrounding  free  field.  Thus  the  net 
overall  effect  of  soil-structure  interaction  is  to  increase  the  rise  time  of 
structure  interior  velocity.  Roughly,  this  rise  time  increase  is  equal  to  the 
time  required  for  the  wave  to  engulf  the  structure,  a  few  milliseconds  for  silos 
and  one  to  two  tenths  of  a  second  for  large  command  centers.  Evidently,  from  the 
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standpoint  of  shock  isolation  the  soil-structure  interaction  becomes  more 
important  with  increasing  structure  size. 

Consider  now  the  basic  problems  to  be  faced  in  designing  an  isolator.  We 
will  not  complicate  the  discussion  by  considering  damping,  inelastic  effects,  etc. 

The  first  and  possibly  most  important  questions  to  be  answered  are  "Is  any 
shock  isolation  required;  does  the  structure  provide  enough  for  this  equipment?" 
These  questions,  we  will  find,  can  be  the  most  difficult  to  answer;  here  we 
assume  that  isolation  is  required. 

The  ultimate  isolator  would  have  zero  natural  -frequency  thus  supporting 
the  suspended  equipment  in  an  absolutely  motionless  state  while  allowing  the 
supporting  structure  ter  mo  ve~  around  It.  At  the  present  ^statebl^techhoTogy  such 
isolators  are  very  nearly  achievable  at  not  unreasonable  cost.  However,  the 
rattle  space  required  by  such  an  isolator  is  the  maximum  occupied  by  a  reasonable 
design  (an  unreasonable  design  could  increase  it).  If  some  acceleration  of  the 
suspended  equipment  could  be  tolerated  the  rattle  space  could  be  decreased. 

Since  Lecrease  of  rattle  space  may  be  reflected  as  savings  in  cost  of  the  entire 
installation  a  knowledge  of  the  tradeoff  between  isolator  rattle  space  required 
and  acceleration  of  equipment  to  be  tolerated  can  be  very  usexML. 

The  tradeoff  criteria  need  not  be  highly  accurate.  Decrease  of  rattle  space 
from  12  inches  to  11  inches  would  not  ultimately  appear  as  a  dollar  savings; 
decrease  from  12  inches  to  6  inches  however  might  justify  a  reduction  in  structure 
size  and  be  reflected  as  a  real  dollar  value. 

The  peak  relative  response  spectrum  of  the  forcing  function  described  above 
fairly  closely  follows  che  conventional  trapezoidal  spectrum  used  in  missile 
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shock  isolation  design.  On  the  left  hand,  or  displacement  leg  of  the  spectrum, 
changes  in  isolation  frequency  will  vary  the  acceleration  suffered  by  the  suspended 
equipment  but  the  rattle  space  required  would  remain  invarient.  On  the  horizontal, 
or  velocity  branch  of  the  spectrum,  the  product  of  displacement  and  acceleration 
remains  invarient  with  changes  of  frequency. 

Reference  to  Figures  2.17  to  2.22  of  the  Newmark-Hansen  report  shows  that 
the  "velocity  plateau"  of  spectra  for  simple  pulses  of  duration  t^  and  having  the 
general  character  of  the  ground  shock  induced  velocity  extends  across  the  range 
of  t^f  between  about  O.U  to  1.0  and  has  a  maximum  pseudo  velocity  about  1-1/2 
times  the  maximum  velocity  of  the  input  pulse. 


Thus  if  rattle  space  is  to  be  decreased  by  employing  a  stiffer  isolation 


system  the  minimum  isolator  frequency  that  could  be  used  would  be  on  the  order  of 
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min  t . 


cycles  per  s^uond 


(3-D 


The  peax  absolute  acceleration,  a,  in  gravity  units  would  be  given  by 

a  „  ^  z  f  2 
g  max  min 


(3.2) 


where  zn|ay  is  peak  relative  displacement. 

Now  peak  soil  particle  velocity  is  given  by 
pa 

V  =  — 

Pc 

where 


(3-3) 


Pa  is  peak  overpressure  and 

Pc  is  the  average  acoustic  impedance  of  the  soil. 


The  maximum  relative  displacement,  z  ,  is  equal  to  the  maximum  soil  displace- 
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ment  for  the  conditions  being  considered.  Thus 


=  F*  t ,  - 


fip  ta 
™ar  d 


(3.^) 


max  '  d  pc 

where  /3  is  a  factor,  somewhat  less  than  one,  so  chosen  that  /SVt^  equals  the 
area  under  the  velocity  curve. 


Also,  the  duration,  t^,  is  given  approximately  by 
d  c 

where  D  is  the  depth  of  soil  between  the  isolator  and  bedrock. 
Combining  these  approximations  yields 


(3.5) 


^50/3  [Pa(psi)J 

“(k)  =  [-y(pcf)]  [D(ft7! 


(3-6) 


where  7,  the  unit  weight  of  the  soil  has  been  substituted  for  pg. 


Now,  it  should  be  noted  that  the  three  parameters  remaining  in  Equation 
(3.6)  are  those  most  accurately  known,  the  wave  velocity,  c,  having  cancelled 
out.  Therefore  Equation  (3.6)  gives  an  excellent  estimate  of  the  minimum 
acceleration  that  must  be  withstood  if  designs  are  to  be  made  on  the  horizontal 


branch  of  the 

shock  spectrum. 

If  we  take 

7 

=  110  lbs/ft3 

D 

=  100  ft 

P 

=  0.85 

we  obtain 

Pa  100  psi  300  psi 

500  psi 

1000  psi 

a  3.5 

g  10.5  g 

17.4  g 

3^.8  g 
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We  conclude  immediately  that  the  experience  already  obtained  on  hard 
installations,  which  indicates  that  shock  isolation  systems  are  almost  univer¬ 
sally  designed  for  frequencies  on  the  displacement  leg  of  the  spectrum,  is 
indicative  of  design  conditions  for  fut  ue  higher  pressure  installations.  Only 
small  ruggedized  equipment  is  likely  to  be  installed  at  the  higher  g  levels. 

Thus  for  most  work  we  will  have  to  use  low  frequency  systems  designed  for  rattle 
space  equal  to  the  entire  soil  motion. 

Accepting  that  the  systems  used  will  he.ve  to  be  low  frequency  maximum  input 
displacement  systems,  we  ask  ourselves  if  there  is  an  optimum  frequency  on  the 
displacement  leg  of  the  spectrum.  It  is  fairly  simple  to  show  that  there  is. 

For  a  first  order  estimate  the  cost  of  a  spring  suspension  system  is  propor¬ 
tional  to  the  weight  of  the  springs  provided.  The  weight  of  the  springs  in  turn 
is  proportional  to  their  energy  capacity. 

Wow  the  energy  capacity  of  linear  springs  is  equal  to  one  half  the  product  of 
the  total  stroke  of  the  springs  and  the  force  developed  at  total  stroke.  The 
total  stroke  is  equal  to  the  equilibrium  stroke  ("static  deflection",  or  distance 
the  spring  compresses  when  the  equipment  to  be  isolated  is  placed  upon  it)  plus 
the  working  stroke. 

At  the  equilibrium  position  the  force  exerted  by  the  springs  is  equal  to  the 
suspended  weight,  i.e.,  in  proportion  to  one  g  acceleration;  at  total  stroke  the 
force  is  proportional  to  ng  acceleration,  n  being  the  maximum  design  acceleration 
of  the  system  including  gravity. 
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Denote 


u 

=  total  spring  energy 

X 

=  equilibrium  stroke 

z 

--  working  stroke 

w 

=  weight  of  isolated  equipment 

Ft 

=  spring  force  due  to  total  stroke 

F 

X 

=  spring  force  due  to  equilibrium  stroke 

k 

=  spring  constant 

Then 

W  = 

1! 

if 

(3.7) 

Ft  =  k(x  +  z) 

(3.8) 

Ft  ...  x  +  z 

F  x 

X 

(3.9) 

u  =  |  Ft  (x  +  z)  =  |  k  (x  +  z)2 

(3-10) 

TT  1  2 

U  wxn 

(3.H) 

Now. 

,  for  a  given  system  w  and  z  are  invarient  but  x  depends 

on  the  value 

n  chosen 

.  Then  for  minimum  energy 

dU  ~  1  r  ^  ^  ~\ 

—  =  0  =  ^  w  [n  4  2nx  -r-  J 
dx  2  dx 

(3.12) 

A  „  dn  A 

0  =  n  +  2  —  x 
dx 

(3.13) 

where 

A 

n 

is  the  value  of  n 

A 

and  x 

is  the  value  of  x  satisfying  Equation  (3.13). 

^7 


From  Equation  (3.9) 


dn  _  1  (x  +  z) 

dx  ~  x  2 

x 

Substituting  Equations  (3.9)  and  (3-1*0  into  (3.I3) 

0  „  *  +  -z  +  2  .  2  &  +  z) 
x  x 

Solving  yields 

/\ 

X  =  z 

A 

and  n  =  2 


(3-1*0 


(3-15) 


(3-16) 

(3-17) 


Thus  the  minimum  energy  vertical  linear  spring  system  will  be  obtained  if 
the  spring  cor '  bant  is  chosen  so  that  the  equilibrium  stroke  is  equal  to  the 
working  stroke.  The  total  acceleration  including  gravity  on  the  isolated  equip¬ 
ment  during  passage  of  the  blast  wave  then  will  vary  between  0  and  2  g. 

Since  the  equilibrium  energy  of  a  lateral  isolation  system  can  be  made 
zero  the  energy  capacity  of  the  lateral  system  can  be  made  indefinitely  low  by 
decreasing  the  system  frequency,  and  cost  will  be  purely  a  function  of  the 
mechanical  devices  elected  by  the  designer. 

3.3  Available  Accuracy 

The  uncertainty  in  input  parameters  limits  the  accuracy  with  which  shock 
isolation  computations  can  be  carried  out. 

If  criteria  are  available  for  establishing  the  limits  of  accuracy  available 
they  can  be  used  to  establish  the  degree  of  complexity  that  will  be  meaningful 
in  carrying  out  the  computations. 


It  is  shown  in  Section  3.2  that,  if  isolators  are  required,  in  the  majority 


of  instances  low  frequency  systems  are  indicated.  In  fact,  the  minimum  energy 
(and  prod  ally  minimum  cost)  systems  have  the  natural  frequency  of  oscillation 
so  chosen  that  the  active  stroke  of  the  isolator  is  equal  to  the  equilibrium 
stroke . 


The  period  of  vibration  of  a  linear  mass  spring  system,  T,  is  given  by 

(3.18) 

where  x  is  the  equilibrium  stroke,  i.e.,  the  static  displacement. 


T  =  2^/f 


For  the  one  g  system  naving  equilibrium  stroke  equal  to  working  stroke 

X  -  _  (3.19) 


max 


Substituting  Equation  (3.5)  into  (3.4)  yields 


x  =  z 


max 


3V 

p=2 


Tnerefore 


2jSpffD 

T“ 


(3-20) 


(3.21) 


Dividing  this  by  the  approximate  duration  of  the  velocity  pulse,  Equation  (3-5) 


T_ 
t . 


(3-22) 


We  note  that  again  the  more  inaccurate  parameters  have  cancelled  out  and  that 
Equation  (3.22)  therefore  gives  an  excellent  estimate  of  the  ratio  T/t^. 

Choosing 

j3  =  O.85 

7  =  110  lbs/ft3 

D  *  100  ft,  200  ft 
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for 


Pa  100  psi  300  psi  500  psi  1000  psi 

T  D  =  100  ft  k.9  8.5  11  15. k 

*4  D  =  200  ft  3.5  6  7. 7  12 

We  conclude  that  most  systems  will  be  designed  by  impulse  methods  requiring 
for  their  input  only  the  peak  displacement  of  the  soil.  The  exceptions  would 
occur  at  the  lower  overpressure  levels,  systems  in  structures  large  enough  to 
considerably  modify  the  input  motion  by  soil-structure  interaction  and  possibly 

1  systems  in  which  isolation  devices  were  omitted. 

* 

t 

,  Therefore  the  limiting  accuracy  available  for  design  of  most  isolation 

A  systems  can  be  determined  by  considering  the  accuracy  with  which  the  peak  input 
'  displacement  can  be  determined. 


First  consider  the  approximate  Equation  3.20. 

z  =  - — 

max  2 

pc 

The  relative  error  Az/ z  is  given  by 

A2=Aj6+AD_4_p_  2Ac 
z  fi  D  p  c 


(3.20) 


(3.23) 


Now  the  average  value  of  p  will  be  known  to  within  two  or  three  percent 
accuracy  if  soil  samples  are  taken.  The  depth  D  can  be  measured  to  within  a 
couple  of  feet  by  seismic  means  but  most  likely  will  vary  from  point  to  point 
at  a  particular  site  resulting  in  a  larger  average  error.  Probably,  if  the 
average  value  could  be  determined  to  -  5  percent  (100  ft  -  5  ft)  we  would  have 
to  be  satisfied. 
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•  ^  f 


The  average  wave  velocity,  c,  reflects  the  nonlinearities  of  the  stress- 
strain  curve  of  the  soil,  the  soil  hysteresis,  and  in  Equation  3.20  all  of  the 
other  approximations  to  the  "elastic  constants"  of  soil.  (We  note  that  the 
denominator  of  Equation  3*20  is  equivalent  to  the  modulus  of  elasticity  of  the 
soil,  possibly  the  least  accurately  determinable  soil  characteristic.)  One  is 
hard  put  to  make  even  a  good  estimate  of  the  accuracy  with  which  c  can  be 
determined.  However,  the  seismic  velocity,  c^,  can  be  measured  to  two  significant 
figures.  Wilson  (Reference  12)  suggests  that  the  velocity  of  the  peak  pressure 
be  taken  as  l/2  to  3/4  of  the  vertical  seismic  velocity.  Accepting  Wilson’s 
estimate  we  could  take  for  the  compression  wave  velocity  5/8  -  l/8  of  the  seismic 


velocity  of  the  soil,  an  inherent  error  of  -  20  percent. 

If  the  overpressure  level  is  in  the  vicinity  of  the  plateau  of  the  soil 
stress  strain  curve  the  reflected  wave  will  travel  faster  than  the  incident  wave 
due  to  the  precompression  Induced  by  the  incident  wave.  Because  of  the  large 
inherent  error  it  does  not  seem  reasonable  to  attempt  to  distinguish  between 
these  velocities  In  computation.  Therefore  :,.t  is  suggested  that,  lacking  actual 
test  data  the  average  compression  wave  velocity  be  taken  as  equal  to  66  -  15 
percent  of  the  uphole  seismic  velocity  recognizing  that  the  error  given  is 
approximate . 

In  circumstances  where  the  reflected  wave  does  not  add  to  the  incident  wave 
it  is  suggested  that,  lacking  test  data,  the  compression  wave  velocity  be  taken 
as  60  -  20  percent  of  the  uphole  seismic  velocity. 
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At  the  high  overpressure  levels  (>  1000  psi)  some  accuracy  in  the  calcula¬ 
tions  can  be  gained  by  making  use  of  the  fact  that  by  the  time  the  reflected 
wave  arrives  the  overpressure  will  have  decayed  to  a  small  fraction  of  its  peak 
value. 


Maximum  surface  displacement  will  occur  at  the  instant  of  peak  reflected 
pressure  arrival  except  in  the  rare  instances  where  soil  attenuation  is  so  great 
that  a  point  is  reached  where  rate  of  attenuation  of  pressure  just  balances  the 
rate  of  increase  of  impulse  delivered  to  the  soil  by  the  aeroshock  wave. 


Excepting  the  latter  condition  the  maximum  displacement  is  given  by 


x 


(3.24) 


where 


x  is  the  peak  displacement 

p  is  the  average  soil  density 

c  is  the  average  compression  wave  velocity 
I.  is  the  impulse  given  by 

IE 

Ia  =  /C  ?'(t)  dt  (3.25) 

o 

D  is  the  distance  from  structure  to  bedrock 
and  p(t)  is  the  pressure  in 'the  soil  including  the  effects 
of  attenuation  with  distance  traveled. 


Now  if  2D/c  is  a  time  sufficiently  long  that  it  represents  a  major  fraction 
of  the  positive  phase  duration  the  impulse  I  will  be  very  insensitive  to  D  and  c 
From  the  elementary  differentiation  formulas  the  relative  error  in  x  is  given  by 


(3.26) 


^2 


* 


I^  being  presumed  accurately  known  when  formula  3.24  is  used. 

Now  the  density  of  a  soil  is  probably  its  most  accurately  known  quantity; 
the  ratio  A p/  p usually  can  be  taken  as  no  more  than  5  percent  if  the  soil  at  the 
site  is  known.  The  wave  propagation  velocity  however  is  known  with  considerably 
less  accuracy,  particularly  at  higher  overpressures  and  for  the  more  inelastic 
soils.  Indeed,  for  real  soils  it  can  only  be  an  approximate,  average  estimate  of 
a  variable  quantity  of  considerable  complexity.  An  estimate  of  -  20  percent  for 
Ac/ c  generally  would  be  quite  good. 

Thus  the  percent  error  in  displacement  resulting  from  use  of  the  simple 

+ 

formula  and  inherent  input  data  error  is  on  the  order  of  -  25  percent. 


Some  en,  ineers  might  question  the  accuracy  of  the  formula  because  it  is 
based  on  elastic  theory  and  is  being  used,  with  appropriate  parameter  approxi¬ 
mations,  for  an  inelastic  medium.  They  might  argue,  and  quite  reasonably,  that 
for  the  inelastic  medium  there  is  not  justification  for  taking 

E  =  pc2  (3.27) 


which  is  done  in  the  development  of  Equation  3.24. 
to  use 


x 


E 


Instead,  they  might  prefer 


(3-28) 


which  is  more  primative  than  Equation  3.24,  and  estimate  an  average  effective 
value  of  E  independent  of  propagation  velocity. 


The  relative  error  u  x  resulting  from  Equation  3.28  is 

Ax  _  Ac  A_E 

x  c  E 


(3.29) 
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Now  there  is  probably  no  more  difficult  quantity  to  measure  and  estimate 
for  an  inelastic  soil  than  the  effective  modulus  of  elasticity.  Even  with 
laboratory  tests  accuracy  better  than  about  -  20  percent  is  -unlikely.  Thus  by 
the  use  of  Equation  3.28  the  error  in  x  would  be  -  4o  percent,  50  percent  higher 
than  would  be  obtained  by  use  of  Equation  3.24. 


SECTION  4 


DETERMINATION  OF  FREE  FIELD  INPUT 


In  Section  3  It  is  shown  that  an  error  on  the  order  of  -  20  percent  must 
be  tolerated  in  the  determination  of  input  to  shock  isolation  systems.  Indeed, 
in  many  real  situations  even  this  accuracy  could  bo  obtained  only  with  the  aid 
of  an  extensive  field  and  laboratory  test  program. 

It  also  is  pointed  out  in  Section  3  that  if  shock  isolation  systems  are 
v  ,ed  they  will  be  low  frequency  systems.  In  relatively  small  structures  the 
influence  of  soil-structure  interaction  on  system  design  will  be  nominal  with 
a  low  frequency  system.  In  subsequent  sections  of  this  report  it  is  shown  that 
the  accuracy  obtainable  with  existing  simplified  soil-structure  interaction 
theories  is  considerably  better  than  the  input  error  that  must  be  tolerated. 

Now,  since  the  actual  design  of  an  isolator  to  resist  an  established  input 
can  be  carried  out  with  scientific  precision,  if  desired,  the  principal  source 
of  error  in  the  entire  design  problem  is  the  error  in  the  input. 

Inasmuch  as  errors  on  the  order  of  -  20  to  30  percent  will  be  involved  n'n  a 
shock  isolator  design,  a  strong  case  could  be  made  for  adopting  some  simple 
convention  and  be  done  with  it.  It  can  be  argued,  quite  reasonably,  that  if  we 
must  guess  the  input  we  may  as  well  guess  the  result. 

Two  counter  arguments  can  be  advanced: 

1.  As  a  consequence  the  possibility  that  shock  isolation  can  be 
eliminated  from  a  particular  system  could  not  be  considered. 


2.  Such  an  approach  would  e  actively  stop  further  progress  in 
the  field  of  nuclear  shock  isolation. 

Therefore  it  is  believed  to  be  necessary  for  this  report  to  present  methods  that 
are  rational  and  as  accurate  as  practical  within  the  present  state  of  the  art 
for  determination  of  the  free  field  ground  shock  input  characteristics.  Since 
the  accuracy  level  of  a  given  design  will  reflect  the  accuracy  of  the  input,  this 
phase  of  the  shock  isolation  design  problem  deserves  a  greater  portion  of  the 
designer's  time  than  any  other  aspect  of  the  overall  problem. 

A  considerable  program  in  the  general  area  of  nuclear  blast  induced  ground 
motions  is  currently  being  carried  out  by  a  number  of  government  agencies  and 

contractors.  Already  a  substantial  literature  on  the  subject  has  developed. _ 

However,  to  date  a  single  consistent  theory  covering  the  several  elements  of 
ground  motion  ary*  real  soils  has  not  been  forthcoming.  Even  a  digest  of  the 
theoretical  and  experimental  data  already  developed  would  be  of  substantial  volume 
and,  because  of  gaps  in  the  data  and  internal  contradictions  would  be  of  limited 
value  to  designers. 

In  view  of  this  situation  the  estimation  procedures  contained  in  this 
section  must  be  regarded  as  tentative,  to  be  replaced  with  more  accurate  proce¬ 
dures  as  the  results  of  research  effort  become  available.  However  it  is  believed 
that  since  a  substantial  error  must  be  tolerated  in  any  event,  a  single  internally 
consistent  approach  that  covers  the  major  influences  significant  to  the  shock 
isolation  problem  is  desirable.  In  particular,  if  the  effects  of  soil-structure 
interaction  and  the  possibility  that  isolators  might  be  omitted  are  to  be 
considered  in  any  detail,  ground  shock  input  wave  forms  are  required.  Since 


practically  all  of  the  published  methods  for  making  ground  shock  estimates  are 
directed  solely  toward  estimation  of  the  three  numerical  input  parameters  required 
for  construction  of  a  trapezoidal  shock  spectrum,  considerable  departure  from 
these  methods  has  been  necessary. 

The  difficulties  inherent  in  the  theoretical  prediction  of  ground  waves  in 
real  soils  ultimately  hinge  on  the  nonlinear,  inelastic  behavior  of  the  soils. 
Propagation  velocity,  energy  attenuation,  wave  diffusion,  and  permanent  displace¬ 
ment  all  are  primarily  influenced  by  these  characteristics. 

A  procedure  for  determining  the  significant  characteristics  of  the  input 
wave  that  requires  a  minimum  of  soil  data  is  presented  subsequently.  T  '-a 
required  are 

1.  the  stress-strain  curve  of  the  soil  (confined  compression  test) 

2.  the  seismic  velocities  of  the  soil 

3.  the  soil  density 

k.  the  percyclic  damping  factor 

5 •  bedrock  depth . 

It  is  assumed  that  the  dissipation  process  conforms  with  the  percyclic 
damping  theory.  This  theory  postulates  that  a  sinusoidal  wave  decays  exponentially 
at  a  certain  rate  per  cycle.  Thus  the  fraction  of  the  amplitude  that  remains 

after  propagating  through  a  distance  of  one  wave  length,  as  compared  with  the 

-  5 

initial  amplitude  of  the  wave,  is  a  constant  and  is  denoted  by  e  .  The  theory 
assumes  that  the  speed  of  propagation  and  the  percyclic  decay  rate  for  sinusoida] 
waves  are  independent  of  wave  frequency  and  stress  level  and  dependent  only  on 
material  properties.  For  a  waveform  that  is  not  sinusoidal,  the  waveform  is 


represented  by  a  Fourier  integral  of  sinusoidal  components,  each  of  which  is 
assumed  to  obey  the  pe "cyclic  damping  law  independently.  Each  component,  having 
a  different  wave  length,  decays  a  different  amount  in  traveling  a  given  distance 
because  it  travels  a  different  number  of  wavelengths.  The  waveform  at  any  later 
time  is  represented  by  a  Fourier  integral  of  the  decaying  sinusoidal  components. 

The  percyclic  damping  theory,  though  based  on  artificial  assumptions,  is 
found  to  agree  well  with  experimental  data  and  has  found  acceptance  by  leading 
investigators  in  the  field  of  stress  wave  propagation.  Values  cf  the  logarithmic 
detriment^ ^ranging  approximately  from  0.2  to  0.5  have  been  reported  for  granular 
media.  Some  typical  values  are  given  in  Section  8. 

The  design  procedure  has  been  simplified  somewhat  by  making  the  assumption 
that  bedrock  is  infinitely  rigid  and  perfectly  reflecting.  The  combined  effects 
of  these  assumptions  is  slightly  conservative .  That  is,  the  motions  predicted 
are  somewhat  greater  than  would  be  predicted  if  the  exact  physical  characteristics 
of  the  bedrock  were  taken  into  account . 

Regardless  of  the  energy  dissipation  characteristics  of  the  soi],  the 
momentum  of  the  ground  shock  wave  is  conserved  as  it  proceeds  through  the  soil. 
Ultimately  all  of  the  momentum  enters  the  underlying  bedrock.  If  the  effects 
of  geometric  dispersion  of  the  wave  are  neglected,  the  displacement  of  the  bedrock 
due  to  the  momentum  can  be  computed  from  the  formula 


where 

y  is  the  bedrock  displacement 
I  is  the  impulse  of  the  ground  wave. 


Though  the  waveform  may  he  considerably  changed  by  passage  through  the  soil 
and  rock,  since  its  momentum  is  conserved  its  impulse  is  conserved  and  is  equal 


to  the  impulse  of  the  overpressure  wave. 

At  the  1000  psi  level,  for  a  bedrock  weighing  l60  pounds  per  cubic  ft.  and 
having  a  seismic  velocity  of  10,000  ft  per  second,  a  one  megaton  bomb  would 
produce  a  rock  surface  displacement  of  about  1.9  inches,  geometric  dispersion, 
being  neglected;  a  20  megaton  bomb  would  produce  a  rock  surface  displacement  of 
about  5*2  inches. 

However,  in  the  bedrock  the  wavelength  of  the  positive  phase  from  the  one 
megaton  bomb  would  be  greater  than  12,000  feet.  The  initial  root  of  the  product 
of  the  two  principal  radii  of  curvature  of  the  shock  front  at  the  1000  psi  region  ' 
would  be  about  7,000  feet,  dependent  on  the  soil  wave  propagation  velocity. 
Therefore,  geometric  attenuation  would  reduce  the  surface  displacement  to  about 
l/*4.  of  the  value  computed  above,  or  about  l/2  inch  for  the  one  megaton  bomb  and 
1-1/2  inches  for  the  20  megaton  bomb.  These  figures  are  very  small  fractions  of 
the  total  displacement  of  soil  surface  with  respect  to  bedrock.  Further,  the 
soil  displacement  computed  as  a  consequence  of  assumed  100  percent  reflection 
-from  bedrock  more  than  compensates  for  these  small  errors. 

4.1  The  Input  Wave  Estimate 

Determination  of  the  input  waves  is  carried  out  in  5  steps: 

1.  The  overpressure  wave  is  estimated. 

2.  The  compression  wave  propagation  velocity  is  estimated. 

1.  The  attenuation  of  the  compression  wave  is  estimated. 


4.  The  compound  wave  resulting  from  addition  of  primary  wave  and 
bedrock  reflection  is  estimated. 

5.  The  vertical  and  horizontal  free  field  input  waves  are  estimated. 

The  above  steps  are  considered  in  order  in  Sections  4.1.1  through  4.1.5. 

Where  desirable  a  brief  discussion  of  the  theory  is  included  though  in  the  main 
the  theory  has  been  described  in  the  references  and  is  not  repeated. 

The  reader  is  cautioned  that  the:  methods  presented  here  are  approximate 
Studies  now  under  way  will  provide  more  accurate  data  on  subsurface  waveforms . 
Because  of  the  importance  of  input  data  accuracy  these  data  should  he  used  for 
applicable  conditions  when  they  become  available. 

4.1.1  The  Overpressure  Wave 

In  a  series  of  classified  reports  Brode  has  presented  theoretically  determined 
overpressure  waveforms  resulting  from  a  surface  burst.  In  Reference  2  Brode 
presents  overpressure  data  in  unclassified  form.  Figure  4.1  was  constructed  from 
the  unclassified  data,  some  extrapolation  being  necessary  at  the  very  high  over¬ 
pressure  levels.  The  curves  on  Figure  4.1  enable  a  5  point  fit  of  the  overpressure 
wave  between  the  limits  of  peak  overpressure  and  10  percent  of  peak  overpressure. 

The  dashed  curve  gives  the  subt  argent  diration,  the  time  at  which  the 
overpressure  has  decayed  to  l/e  times  peak  overpressure.  If  the  initial  portion 
of  the  blast  wave  were  a  true  exponential,  the  subt argent  duration  would  be  equal 
to  the  time  interval  subtended  by  the  tangent  to  the  overpressure  curve  at  peak 
overpressure . 
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Within  the  intervals  between  the  points  given  in  Figure  4.1  the  waveform 
can  be  taken  as  an  exponential,  the  decay  constants  differing  from  interval  to 
interval.  The  notation  and  formulas  are  given  in  Figure  4.2. 

It  should  be  noted  that  the  first  derivative  of  the  approximate  curve 
jumps  at  the  transition  points.  Therefore  accelerations  should  not  be  computed 
from  the  approximation. 

4.1.2  Wave  Velocities 

r. 

The  average  modulus  of  elasticity  and,  equivalently,  average  wave  velocities 

f 

A 

|  of  a  real  soil  contribute  the  greatest  factor  of  uncertainty  to  shock  isolation 
design.  The  material  in  this  section  is  presented  in  an  effort  to  enable  the 

» 

*  designer  to  improve  his  estimates.  It  is  emphasized  however  that  no  tests  have 

:  been  conducted  to  prove  or  disprove  the  method  though  data  reported  for  other 

| 

purposes  do  not  contradict  the  approach. 

| 

|  In  order  to  carry  out  the  estimation  procedure  the  data  previously  listed 

i 

are  required. 

Even  if  these  data  are  only  estimates  themselves  some  improvement  in  the 
|  accuracy  of  the  input  can  be  achieved  over  that  obtainable  by  simply  taking  the 

j  compression  wave  velocity  as  l/2  to  3/^  of  an  estimated  seismic  velocity. 

|  Figure  4.3a  schematically  shows  a  stress-strain  curve  for  a  sandy  soil. 

» 

!  Now,  a  pressure  wave  initially  of  the  same  waveform  as  the  overpressure  wave 
|  would  be  dispersed  in  propagating  through  the  medium,  that  is  the  peak  would  be 

|  attenuated  and  would  lag  behind  the  "toe"  of  the  wave  which  would  be  of  zero 

amplitude  and  propagate  at  seismic  velocity  c^.  Between  the  toe  of  the  wave  and 
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Figure  4.2  APPROXIMATION  OF  BLAST  WAVE  3Y  EXPONENTIAL  INCREMENTS 


Figure  4.3  CONFINED  COMPRESSION  STRESS  VERSUS  STRAIN 
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the  maximum  the  pressure  would  he  continually  increasing.  Therefore,  each  element 
of  the  rise  would  propagate  at  the  instantaneous  velocity  c*  defined  as  shown  in 
Figure  4.3a. 


The  relation  between  instantaneous  wave  propagation  velocity,  stress,  strain, 
and  density  is  given  as  follows: 

(4.2) 


c*2  =  I_  da 

p  de 
r  o. 


where 


c*  =  infinitesimal  disturbance  propagation  velocity  with 

respect  to  particle  motion  at  a  pressure 
Pq  *  density  at  zero  stress 

a  =  stress 

e  ss  engineering  strain 

This  is  derived  from  continuity  and  momentum  conservation  only  and  is  independent  of 
stress  strain  laws.  It  is  valid  everywhere  except  at  a  finite  discontinuity  of 
stress,  i.e.,  a  shock  front.  It  should  be  noted  that  c*  Is  the  propagation 
velocity  measured  in  a  coordinate  system  moving  at  the  local  particle  velocity. 
Therefore,  to  obtain  the  wavefront  velocity  in  a  fixed  coordinate  system  the 
particle  velocity,  v,  should  be  added.  This  quantity  is  given  by 


v  =  ; 


da 

pc* 


(4.3) 


where 

P 


c* 


density  at  stress  a 

wave  velocity  with  respect  to  particle  velocity  at  stress  a 


Evidently,  the  particle  velocity  computation  is  lengthy  for  a  nonlinear  soil. 


Since  it  is  small  in  comparison  with  c*,less  than  5  percent  generally,  the 
computation  of  c  can  he  greatly  simplified  with  tolerable  error  by  approximating 
c  =  c*. 

Evidently,  since  the  pressure  peak  is  continuously  attenuated  it  will 
propagate  at  varying  velocity.  For  the  conditions  shown  in  Figure  4.3a  the 
peak  would  propagate  with  continuously  increasing  velocity.  It  should  be  noted 
that  although  the  peak  of  the  wave  is  continuously  attenuated  as  it  propagates 
through  the  ground,  at  any  given  point  in  the  soil  the  pressure  continuously 
rises  to  the  peak  without  going  through  a  hysteresis  loop.  Therefore,  the  peak 
of  pressure  always  propagates  at  the  velocity  defined  by  the  stress-strain  curve 
envelope  ABC.  -  - 

A  plot  of  the  ratio  c/c^  can  be  constructed  by  measuring  the  slopes  of  the 
stress-strain  curve  and  taking  the  square  roots  of  the  ratios  of  these  to  the 
slope  at  the  origin.  Since  the  slopes  orily  are  measured  and  the  plotted  results 
represent  the  roots  of  these,  the  error  in  the  ratio  c/c^  thus  obtained  is  less 
than  half  of  the  error  in  stress  at  a  given  strain  as  given  by  the  stress-strain 
curve.  The  dashed  curve  ADF  of  Figure  4.4  would  yield  a  velocity  ratio  curve 
practically  identical  with  that  resulting  from  curve  ABC. 

Since  the  seismic  wave  velocity  can  be  measured  to  within  -  5  percent  or 
better  accuracy  the  compression  wave  velocity  can  be  estimated  by  this  technique 
to  within  considerably  better  limits  than  the  -  20  percent  resulting  by  taking 
an  average  ratio. 


Three  pressures  are  identified  in  Figure  4.3a: 

1.  the  peak  overpressure, 

2.  p  peak  attenuated  pressure  at  a  buried  structure, 

8. 

3.  p^  peak  attenuated  reflected  pressure  arriving  at  the  structure. 

When  these  three  pressure  values  have  been  estimated  the  corresponding  instantaneous 
compression  wave  velocities  can  be  taken  from  a  graph  similar  to  Figure  4.3b.  An 
average  value  then  can  be  chosen  for  computation  purposes. 

In  order  to  estimate  wave  attenuation  by  the  method  of  Section  4.1.3  the 
compression  wave  propagation  velocity  is  required.  Thus  the  method  becomes  one 
of  guess  and  check.  An  original  guess  corrected  once  should  be  sufficient  to 
arrive  at  consistent  propagation  velocity  mid  wave  attenuation  faetors.  - 

The  arguments  leading  to  the  estimation  procedure  given  above  are  valid  as 
long  as  the  entire  ground  shock  wave  lies  on  a  portion  of  the  soil  stress-strain 
curve  that  is  everywhere  convex  upward. 

If  the  stress-strain  curve  of  the  30il  is  concave  upward  the  propagation 
phenomena  are  qualitatively  changed.  Inasmuch  as  all  soils  exhibit  a  concave 
stress  strain  curve  at  high  enough  pressure  this  problem  merits  some  consideration. 

Figure  4.5  illustrates  schematically  a  generalized  stress-strain  curve 
ABCDEF.  The  curve  is  convex  from  A  to  the  point  of  inflection,  D,  and  concave 
beyond  D. 

In  general  we  should  expect  that  if  a  wave  of  amplitude  p^  were  propagating 
into  such  a  material  a  shock  front  would  quickly  form.  If  the  pressure  maximum 
were  preceded  by  a  rising  pressure  the  stresses  within  the  interval  of  the  rise 


Figure  4.5  CONCAVE  STRESS -STRAIN  CURVE 
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time  would  travel  at  lower  velocities  than  that  at  the  peak  and  quickly  be 
overtaken.  Since  the  overpressure  wave  has  an  abrupt  rise  at  the  higher  over¬ 
pressures  we  conclude  that  the  sharp  rise  is  conserved  in  soil  if  the  air 
induced  ground  pressure  is  above  the  point  of  inflection  D.  We  repeat  that 
all  real  soils  exhibit  a  concave  region  similar  to  DEF  at  pressures  above  some 
critical  level. 

Looking  at  curve  ABCDEF  one  is  tempted  to  conclude  that  a  decaying  wave 
having  a  sharp  rise  would  be  stretched  out  in  propagating  through  such  a  soil 
e.g.  pressure  p,-  in  Figure  4.5  closely  following  peak  pressure  p^  would  move 
at  a  slower  speed  and  therefore  the  interval  p^-p^  v*<  uld  grow  with  time.  However, 
the  decay  of  the  wave  also  propagates  according  to  the  peculiar  laws  of  such 
stress-strain  curves  and  actual  physical  results  differ  greatly  from  this 
intuitive  conclusion. 

Kolsky  (Reference  3)  and  others  have  investigated  the  phenomena.  The 
results  of  these  investigations  predict  wave  behavior  that  differ  markedly 
from  the  intuitively  anticipated  behavior  prese  ted  in  the  previous  paragraph. 

The  principal  physical  effects  that  result  in  the  qualitative  peculiarities 
of  the  propagated  wave  are  two: 

1.  The  velocity  of  propagation  of  a  shoes'  front  is  not  given  by  Equation 
(4.2);  the  derivative  is  meaningless  at  the  shock  front  discontinuity. 

2.  The  velocity  of  propagation  of  the  decaying  portion  of  the  wave  follow¬ 
ing  the  shock  front  is  given  by  Equation  (4.2)  but  the  uerivative  should  be 

de 

taken  on  a  hysteresis  loop  (FG  of  Figure  4.5)  rather  than  on  the  envelope. 


TO 


The  accurate  analysis  of  these  phenomena  is  exceedingly  complex  and  is  still 
under  study  by  a  number  of  investigators.  However,  if  two  simplifying  assumptions 
are  ma  s  a  method  for  closely  estimating  resulting  waveforms  is  forthcoming.  One 
simplifying  assumption  reasonably  matches  reality.  It  is 

1.  The  hysteresis  loops  are  shrunk  to  their  centerlines  and  their 
centerlines  are  all  parallel  to  the  tangent  at  zero  stress  and 
strain.  Therefore  the  propagation  velocity  on  a  hysteresis 
loop  is  equal  to  the  seismic  compression  velocity. 

The  second  assumption  may  depart  considerably  from  reality.  It  is 

2.  The  stress-strain  curve  (ABCDEF  of  Figure  4.5)  is  the  Hugoniot 
of  the  soil.  At  very  high  pressures  this  assumption  has  been 
found  to  be  grossly  in  error  for  materials  having  definite, 
definable  properties.  The  error  decreases  as  the  pressure  is 
lowered. 

The  influence  of  the  second  assumption  will  be  considered  first. 

An  elementary  manipulation  of  the  continuity  and  momentum  conservation 
equations  results  in  Equation  4.4  for  piopagation  of  a  shock  front,  if  a  shock 
front  exists. 

~  (l-ea)2  (ae-ai) 

u*  -  ~1T~  ^ 

where 

Os.  and  Ci  are  the  stresses  on  the  high  and  low  pressure  sides  of 
the  shock  front  respectively. 


e2  and  e1  are  the  engineering  strains  on  the  high  and  low  pressure 


sides  of  the  shock  front,  respectively, 

Pq  is  the  density  at  zero  stress, 

u*  is  the  shock  velocity  measured  in  a  coordinate  system 

moving  with  the  particle  velocity  ahead  of  the  shock. 


Again,  only  a  small  error  is  introduced  if  the  shock  velocity  measured  in  a 

2 

fixed  system  is  approximated  by  U*.  Further,  the  term  (1-62)  generally  can  be 
considered  unity  without  serious  error.  Actually,  this  '  is  equivalent  to 

po 

— — .  an<^-  as  equal  to  1  if  changes  in  soil  density  due  to  pressure  are  small 
*2 

enough  to  be  neglected.  Thus  a  slightly  inaccurate  but  reasonable  app-rmri nw  +4  nn 
for  the  shock  velocity  is  given  by 


reasonable  approximation  for  the  shock  velocity  is  given  by 

P  ,  (02 -ai) 

U  ”  PQ  (e2-£iJ 


(*•5) 


Now,  an  arbitrary  curve  can  not  be  taken  as  the  locus  of  all  points 
( 0Q ,  €q),  (oi,  £i),  (b2,  £2)  for  which  Equation  4.4  is  satisfied.  The  correct 
curve,  termed  the  Hugoniot,  along  with  Equation  4.4  is  satisfied  also  satisfies 
the  equation  of  conservation  of  energy.  If  no  heat  were  generated  by  the 
passage  of  the  shock  front  the  static  stress  strain  curve  and  the  Hugoniot  would 
he  identical.  However,  in  any  real  material  heat  is  generated  by  passage  of  a 
shock  front.  Therefore,  for  real  materials,  the  curves  differ. 


However,  beyond  the  near  crater  region  the  heat  generated  by  passage  of  the 
ground  shock  wave  is  small.  Therefore,  the  assumption  that  the  stress  strain 
curve  is  the  Hugoniot  is  not  unreasonable.  The  reader  should  not  be  deluded 
into  believing  that  this  convenient  argument  indicates  that  differences  between 
stress  strain  curves  and  the  Hugoniot  actually  are  negligible.  Possibly  the 
best  counter  argument  is  that  dynamic  stress  strain  curves  lie  between  the  static 
curve  and  the  Hugoniot.  Nevertheless,  it  is  a  far  better  approximation  of  a 
real  soil  than  the  homogenious,  isotropic,  elastic  medium  of  the  theory  of 
elasticity. 

Now  consider  a  wave  having  an  abrupt  rise  in  pressure  from  p^  to  p^  (see 
Figure  4.5).  In  view  of  Equation  4.5  the  propagation  velocity  of  the  dis¬ 
continuity  (shock  front)  would  be  c&  not  cf.  Also  we  note  that  portions  of  the 
wave  at  pressures  below  p^  would  be  propagating  faster  than  the  shock  front. 

This  brings  out  the  first  qualitative  peculiarity  of  wave  propagation  through  a 
medium  having  a  concave  stress  strain  curve.  A  shock  front  is  conserv  .  but  is 
preceded  by  an  ever  increasing  precursor.  Also  we  note  that  if  the  shock  front 
decayed  (to,  say,  p^  of  Figure  4.5)  the  maximum  amplitude  of  the  precursor  would 
increase.  If  the  pressure  decreased  to  pcr  precursor  and  post  shock  wave  would 
merge  and  the  shock  would  vanish. 

Now  consider  the  decaying  portion  of  the  wave  behind  the  shock  front.  It 
follows  a  hysteresis  loop  and  therefore  is  moving  at  higher  velocity.  Thus  the 
decaying  of  the  wave  overtakes  the  shock  front  and  the  latter  is  continuously 
attenuated.  A  consequence  of  the  assumption  that  the  wave  velocity  on  any 
hysteresis  loop  is  equal  to  the  seismic  velocity  is  that  the  rarefaction  wave 


behind  the  shock  front  moves  at  constant,  seismic  velocity  regardless  of  the 
position  of  the  shock  front  on  the  stress  strain  (Kugoniot)  curve.  Actually, 
laboratory  tests  show  that  the  seismic  velocity  does  not  remain  constant  as 
confining  pressure  is  increased.  Generally,  the  seismic  velocity  increases  (e.g. 
Reference  11,  page  57). 

In  Section  4.1.3  it  is  shown  that  the  higher  velocity  rarefaction  wave  rapidly 
attenuates  peak  pressure.  Therefore  the  error  introduced  by  setting  the  "hysteresis 
velocity"  equal  to  the  seismic  velocity  has  the  net  result  of  decreasing  predicted 
attenuation,  an.  error  on  the  side  of  safety. 

In  summary,  the  shock  propagation  velocity  in  a  soil  having  a  concave  upward 

stress  strain  curve  should  be  as  determined  from  the  slope  of  the  secant  drawn  on 

the  stress  strain  curve  between  the  pre-shock  and  post  shock  stresses,  analogous 

to  line  BF  and  velocity  c  of  Figure  4.5.  Since  the  peak  stress  is  attenuated 

8. 

by  the  rarefaction  wave  the  value  of  propagation  velocity  will  vary  with  time. 

4.1.3  Compression  Wave  Attenuation 

The  percyclic  damping  theory  is  based  on  artificial  theoretical  assumptions. 

Its  three  salient  virtues  are: 

1.  It  checks  experimentally  obtained  data  with  excellent  accuracy. 

2.  It  is  simple  to  apply,  to  both  theoretical  and  practical  problems. 

3.  The  damping  constants  for  at  least  a  few  real  soils  have  been 
measured  and  reported  in  the  geophysical  literature.  (These  are 
given  in  Section  8.) 


(4 


The  writer  of  this  report  believes  that  designers  will  agree  with  him  that 
number  3  above  is  the  theory's  most  sterling  virtue. 

The  elements  of  the  theory  are  given  in  Section  4.0  of  this  leport. 

Weiner  (Reference  10 )  has  investigated  the  influence  of  the  theory  on  the 
one  dimensional  propagation  of  simple  pulses.  The  data  he  reports  are  directly 
applicable  to  the  ground  shock  problem.  Figures  4.6  through  4.9  are  taken  from 
Reference  10.  In  these  figures  "n"  is  the  number  of  initial  pulse  widths 
propagated. 

The  dimensionless  equation  for  pulse  form  after  propagation  of  n  initial 
pulse  widths  given  by  Weiner  for  the  initially  triangular  pulse  is 


2  ^2. 


r(s,B)  .  i  (1*23)  tan'1,  "  *  ^  loge 

tr)  +  (4s  “1)  /$n\  + 


+  (2s-l)2 


(4.6) 


(^)  +  (2s+l)‘ 


In  the  above  formula 


r(s,n) 


distance  from  center  of  pulse  measured  in  initial 
pulse  widths 

pulse  propagation  distance  in  initial  pulse  widths 
percyclic  damping  constant  (per  cycle) 

dimensionless  ratio  of  wave  amplitude  to  initial  amplitude 
at  centerline 


A  plot  of  the  ratio  of  peak  amplitude  to  initial  peak  amplitude  as  a  function 
of  pulse  propagation  distance  is  given  in  Figure  4.10. 


?S 


Number  of  Pulse  Widths  Traveled 


Figure  h.Q  DAMPING  OF  TRIANGULAR  HJISE 


Now,  we  note  that  the  percyclic  damping  theory  ie  developed  for  materials 
having  a  constant  wave  propagation  velocity  and  we  intend  to  apply  it  to  nonlinear 
media  having  variable  propagation  velocities.  Objections  certainly  can  be  raised 
but  we  point  out: 

1.  A  more  applicable  theory  of  comparable  simplicity  does  not 
appear  to  be  available. 

2.  The  percyclic  damping  theory  has  given  good  results  in  geo¬ 
physical  work  where  the  media  also  ere  nonlinear. 

In  the  following  paragraphs  a  method  for  estimating  compression  wave 
attenuation  for  soils  having  convex  upward  stress  strain  curves  is  developed, 
based  on  the  percyclic  damping  thc^-^y.  The  attenuation  of  peak  pressure  in  soils 
having  concave  upward  stress  strain  curves  is  so  rapid  that  a  separate  damping 
computation  would  add  no  accur  acy  to  the  waveform  determination. 

Figure  4.11  shows  a  typical  overpressure  wave  constructed  from  Figure  4.1. 

(A  1  ME  bomb  and  500  psi  have  been  used  for  the  construction).  The  wave  is 
equivalent  to  the  sum  of  the  triangvlar  wave  and  the  "difference"  wave  showr. 

Now  the  difference  wave  is  of  long  wave  length  and  any  high  frequency  Fourier 
components  which  it  contains  are  of  low  amplitude.  Therefore  it  will  be  attenuated 
only  slightly  in  passing  through  the  soil  and  the  attenuation  that  does  occur  will 
affect  only  the  point  where  the  trailing  edge  of  the  triangular  wave  joins  it. 

The  triangular  wave,  on  the  other  hand,  has  a  short  wave  length,  and  contains  a 
large  portion  of  its  energy  in  high  frequency  components.  Therefore  it  will  be 
attenuated  rapidly.  We  make  the  approximation  that  all  of  the  attenuation  occurs 
in  the  triangular  component  of  the  wave. 


81 


The  duration  of  the  triangular  wave  has  been  set  arbitrarily.  However,  the 
value  chosen  makes  the  central  angle  of  the  peak  of  the  triangle  equal  to  the 
central  angle  of  the  peak  of  the  blast  wave.  Since  the  high  frequency  components 
of  the  peak  are  most  strongly  attenuated  the  subtangent  duration  appears  to  be 
the  best  choice  for  the  triangular  base. 

To  apply  the  theory: 

1.  Determine  wave  path  length  from  surface  to  structure  and  surface 
to  bedrock  to  structure.  The  angle  between  the  paths  and  vertical 
is  equal  to  the  Mach  angle  of  the  compression  wave. 

2.  Determine  the  percyclic  damping  constant.  Test  data  should  be 
used  if  available;  otherwise  estimate  with  the  aid  of  data  in 
Section  8. 

3.  Determine  the  sub tangent  duration  of  the  overpressure  wave 
from  Figure  4.1, 

4.  Determine  the  subtangent  length  by  multiplying  the  subtangent 
duration  by  the  compression  wave  velocity.  (This  will  be  a 
guess  and  check  procedure  in  connection  with  determination  of 
the  compression  wave  velocity  per  Section  4.1.2). 

5.  Determine  n  for  the  direct  path  and  the  reflected  path. 

6  D. 

6.  Compute  ■—  for  direct  path  and  reflected  path  and  obtain 
attenuation  factors  from  Figure  4.10. 

7.  Determine  the  rise  time  by  computing  delay  from  time  of 
arrival  of  toe  of  wave  (moving  at  velocity  c^)  to  time  of 
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arrival  of  peak  (moving  at  velocity  c).  If  L  is  the  path 


length  and  t  the  rise  time 

*,  - 

The  curve  of  the  rise  can  be  interpolated  between  the  curves 
of  Figure  4.12  (Figure  4.12  has  been  constructed  from  Figure 
4.9  hy  continuing  the  maximum  slope  of  Figure  4.9  curves  to 
the  "s"  axis). 

According  to  the  elastic  theory  lateral  motion  of  a  soil  particle  should 
stop  shortly  after  engulfment  by  the  shear  wave.  Field  test  records  obtained  in 
Nevada  however  do  not  exhibit  this  abrupt  termination  of  the  horizontal  velocity. . 

This  behavior  can  be  rationalized  if  it  is  assumed  that  the  rate  of  energy 
dissipation  of  shear  waves  is  considerably  higher  than  that  of  compression  waves. 
According  to  the  elastic  theory  the  shear  wave  stress  is  considerably  lower  than 
the  compression  wave  stress  in  the  superseismic  region.  Thus,  if  the  shear  wave 
were  attenuated  at  a  greater  rate  than  the  compression  wave  its  intensity  rapidly 
would  be  reduced  below  the  level  of  the  wave  anomalies  introduced  by  the  nonlinear 
characteristics  of  the  soil  and  the  argument  based  on  elastic  theory  loses  its 
force . 


Now,  precise  numbers  are  not  available,  but  it  is  known  that  shearing 
distortion  of  particulate  soil  dissipates  energy  much  more  rapidly  than  pure 
compression.  Thus,  though  the  Nevada  test  site  can  hardly  be  considered  typical 
the  we'-ght  of  evidence  indicates  that,  lacking  test  data  to  the  contrary,  the 
shear  waves  can  he  neglected  in  shock,  isolation  design. 


4.1.4  The  Input  Wave 

The  input  wave  determination  for  convex  and  concave  stress  strain  diagrams 
are  considered  separately.  Specific  examples  are  used  to  illustrate  the  methods. 

Case  i  Stress  Strain  Curve  Convex  Problem 

A  structure  is  to  he  placed  35  feet  below  surface;  bedrock  is  120  feet 
below  surface.  The  design  overpressure  is  300  psi  and  the  design  bomb  size  is 
10  MT.  The  seismic  velocity  is 

C;L  =  2700  ft/sec. 

The  average  soil  density  is  120  lbs/ft^.  The  soil  is  a  lightly  cemented  sand. 

The  percyclic  damping  factor  is  estimated  (Section  8)  to  be  0.20. 

The  average  stress-strain  curve  of  the  soil  is  estimated  to  be  as  shown 
in  Figure  4.13.  [The  soil  chosen  for  this  example  might  almost  be  classified 
as  a  soft  rock.  It  has  been  chosen  to  illustrate  how  a  slightly  curved  stress- 
strain  curve  can  modify  ideal  elastic  results.] 

Develop  the  estimated  pressure  time  curves  at  the  structure  for  both  incident 
and  reflected  waves. 

Step  1  -  Con:  cruet  the  overpressure-time  curve  (Figure  4.l4) 

From  Figure  4.1 

At  300  psi  and  1  MT  At  300  psi  and  10  MT 

P  =  150  psi  tj-g  =  31  ms  t^Q  =  67  ms 

P  =  110  psi  t  ,  ---  55  t  =  118 
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At  300  psi  and  1  MT 


At  300  psi  and  10  MT 


P  =  90  psi 

*30  *  T5 

*30  •  161 

P  =  60  psi 

*20  - 

*20  ' 

P  =  30  psi 

*10  *  360 

t10  =  776 

Step  2 

Construct  the  c/c^  curve  corresponding  to  the  stress-strain  curve. 
(Figure  4. 13b) 

Step  3 

Estimate  the  propagation  velocity.  From  Figure  4.13b  the  range  of  c / c ^ 
between  100  and  300  psi  is  about  O.85  to  O.65.  Accordingly  estimate  c  to  be 
c  =0.75  x  2700  =  2000  ft/sec. 

Step  4 

Estimate  the  wave  attenuation.  At  300  psi  the  air  shock  velocity  is 

4.8  ft/msec.  Therefore  the  first  estimate  of  the  Mach  angle  of  the  wave  is 

,  .  -1  c  .-12  o!i  ,-0 

$  =  sin  y  =  sin  -jj— ^  =  24.6 

and  sec  </>  =1.1. 

The  wave  path  length  from  surface  to  structure  is  about  38  ft.  and  from 
surface  to  bedrock  to  structure  is  about  226  ft.  (see  sketch). 


oy 


4.8  ft/ms 


The  subtangent  duration  of  the  wave  is  (Figure  4.l4)  118  msec.  Therefore 

the  subtangent  length  L  is  about 

s 

L  =  c  t  =  2  x  118  =  236  ft. 
s  s 

The  wave  path  lengths  in  terms  of  the  sul  tangent  length  are 
for  the  direct  wave  n  =  38/236  =  0.l6 

for  the  reflected  wave  n  =  226/236  =0.95 

and  Sn/rr  is  (6  =  0.2) 


for  the  direct  wave 


5S  =  0.01 

7T 

for  the  reflected  wave  =  0.06 

7T 


00 


From  Figure  ^4- .  10  the  attenuation  factors  are 


for  the  direct  wave  95% 

for  the  reflected  wave  71$ 

Peak  attenuated  pressures  are  then 

for  the  direct  wave  285  psi 

for  the  reflected  wave  212  psi 


Step  5  -  Correct  Wave  Velocity 

Entering  uhese  pressures  into  Figures  ^ .13  a  and  h  wi  see  that  a  better 
average  for  c/c^  is  O.65.  Then 

c  =  O.65  x  2700  ft/sec  =  1750  ft/sec 

Recomputing  Lg,  n^,  attenuation  and  peak,  attenuated  pressure  (see  sketch)  we 
obtain 
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L 

s 

n 

6  n 

7T 

Attenuation 

Factor 

Attenuated  Pressure 

Direct  wave 

20 6' 

0.l8 

0.011 

95  $ 

284 

Reflected  wave 

206' 

1.07 

0.07 

69 i> 

207 

These  values  check  the  corrected  wave  propagation  velocity,  1750  ft/sec. 

Step  6  -  Estimate  Rise  Times 

The  rise  times  are  computed  from 

for  the  direct  wave  t  =  —  -  =  7.6  ms 

r  c  c2 

pi  o  pro 

for  the  reflected  wave  t  - - -  42  ms 

r  c  cx 

Step  7  -  Estimate  Arrival  Time  of  Reflected  Wave 


The  distance  traveled  by  the  reflected  wave  from  the  time  that  the  structure 
is  engulfed  by  the  direct  wave  until  the  reflected  wave  reaches  the  structure  is 
about  158  ft.  The  toe  of  the  wave,  moving  at  seismic  velocity  will  travel  this 
distance  in 

=  59  msec. 

Step  8 

Sketch  waveforms  (see  Figure  4.15). 

gteg,.  9 

Determine  maximum  impulse  in  soil  below  structure. 

Time  of  impulse  =  =  90  ms 


AT  STRUCTURE 


The  maximum  impulse  will  be  the  maximum  contained  within  any  90  msec 
period.  This  is  determined  from  Figure  *Kl5  on  the  pressure  sum  curve.  The 
maximum  impulse  is  shaded.  With j n  the  limits  of  accuracy  imposed  by  our  propaga¬ 
tion  velocity,  peak  displacement  is  given  by 

I 

7  Pc 

where 

y  =  peak  displacement 

I  =  impulse  (shaded  area  Figure  4.15) 

p  =  soil  density 

c  =  average  wave  propagation  velocity. 

Case  II  Stress  Strain  Curve  Concave  Problem 

A  structure  is  to  be  placed  35  feet  below  surface;  bedrock  is  120  feet 
below  surface.  The  design  overpressure  is  1000  psi  and  the  design  bomb  size 

O 

is  10  MT.  The  average  soil  density  is  115  lbs/ft  .  An  estimate  of  the  stress- 
strain  curve  below  1000  psi  is  shown  in  Figure  4.l6 

Develop  the  soil  displacement  and  pressure  time  curves  at  the  structure 
level . 

The  basic  one  dimensional  wave  equation  is  independent  of  material  stress 
strain  characteristics.  If  compressive  stresses  are  considered  positive,  the 


wave  equation  is 


=  "  P 


't* 


The  method  suggested  involves  a  series  of  steps  leading  to  a  plot  of 


(fc.7) 


for  a  particular  depth,  h,  as  a  function  of  time  between  t  *  and  t  *,  the  times 

o  s 

at  which  the  wave  front  and  shock  front  respectively  reach  the  point  of  interest 
at  depth  h.  Integration  of  this  plot  yields  particle  velocity,  (^f)  >  during 

the  time  interval  tQ*  <  t*  <  t  *.  A  special  calculation  is  required  to  find 
the  increase  in  momentum  due  to  the  shock,  and  the  influence  of  the  unloading 
portion  of  the  wave  can  be  evaluated  using  linear  elastic  theory,  since  a 
constant  propagation  velocity  is  assumed  in  the  unloading  region. 

Step  1  -  Construct  the  overpressure -time  curve  (See  Figure  4.17) 

Fi  m  Figure  4.1 


At 

1000 

psi  and 

1  MT 

At  1000  psi  and 

500 

psi 

bo' 

8.3  ms 

bo 

=  17.9  ms 
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21.8 

t30 

=  47 

200 

t20  = 

44 

t20 

=  95 
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11 

0 

1 — 1 
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131 

t10 

=  283 
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On  the  soil  stress-strain  diagram,  locate  the  peak  air  overpressure  and 
the  critical  stress  (i.e.,  the  stress  below  which  a  shock  will  not  propagate, 
where  the  stress  strain  curve  has  a  point  of  inflection).  Also  locate  the 
upper  and  lower  points  of  several  intermediate  shocks.  These  steps  have  been 
carried  out  in  Figure  4,16. 


The  following  results  are  obtained  from  data  on  this  figure . 

2 


115  ,  _  ,„-4  lbs  sec 

p  =  wrrm  =  1,72  x  10  — rr“" 


in 


2 

c  = 

o 

4 

10  1000  on 
1.72  X  .002  -  29  X 

a  .  2 

sec ' 

co  = 

4500  ft/sec 

Shock  Peak 

Precursor  Peak 

AP 

Ac 

Shock  Velocity, 

1000  psi 

105  psi 

895  psi 

.047 

870  ft/sec 

750 

135 

615 

.042 

762 

500 

155 

345 

.032 

650 

400 

3  75 

225 

.025 

595 

300 

180 

120 

.018 

515 

200 

200 

000 

.000 

370 

Step  Jj 

Plot  the 

secant  shock  velocity  as  a  function  of  the 

upper  shock  stress 

level.  This  is  shown  in  Figure  4.18. 


Firure  4j  8  RATIO  OF  GEOANT  SHOCK  VEIDCITY  TO  c 


VERSUS  STRESS 


Step  4 


Plot  the  air  blast  overpressure  on  a  function  of  cot.  The  plot  of  pressure 

versus  c  t  is  the  waveform  that  would  exist  in  the  soil  if  the  entire  wave 
o 

propagated  at  velocity  cq.  This  is  shown  in  Figure  4.19. 

Step  5 

Determine  the  time  t*  required  for  the  shock  front  to  penetrate  a  d'  'tance 

x.  Approximate  the  shock  front  velocity  during  a  time  interval  At*  with  the 

average  of  the  secant  velocities  at  the  initial  and  final  shock  pressures  for 

the  interval  At*.  During  the  time  interval  At*  the  toe  of  the  wave  moves  a 

distance  c  At*  and  the  shock  front  moves  a  distance  U.  .  At*  where  U.  .  is  the 
o  ij  ij 

average  secant  velocity  between  times  t*.  and  t*  where 

u  J 


At*  =  t* 


j 


t*± 


Thus  during  the  interval  At*  the  shock  front  moves  backward  with  respect  to 
the  toe  of  the  wave  a  distance 


(c  t.  -  c  t.)  =  A(c  t )  =  (c  -  U. .)At* 
'  o  i  o  j'  'o'  vo  ij 


distance  being  represented  as  cQt  as  used  in  Step  4. 


Solving  for  At* 

At*  = 


t .  -  t . 

_! _ i 

U.  . 
ij 


l  - 


For  any  point  of  closely  spaced  time  instants  t.  and  t.  the  corresponding 

^  J 

peak  pressures  can  be  read  from  the  graph  of  overpressure  versus  time.  The 

ratio  U..  ,/c  may  be  read  from  the  graph  of  Step  3.  Also  the  corresponding 
ij  o 


1.00 


lower  shock  pressure  can  he  read  from  the  graph  of  Step  2,  Thus  a  graph  of 
upper  and  lower  shock  stress  values  plotted  against  time  t*  can  be  constructed. 
The  data  is  shown  in  the  following  table;  the  plot  is  shown  in  Figure  4.20. 


Peak 

t 

At 

Average 

1 

At* 

t* 

Ax 

X 

Pressure 

(psi) 

(ms) 

(ms) 

Shock 

Velocity,  U 
(ft/sec)  ^ 

c 

0 

(ms) 

(ms) 

(ft) 

(ft) 

T50 

9 

9 

8l6 

1.22 

11.0 

11.0 

9.0 

9.0 

500 

13 

9 

706 

1.19 

10.7 

21.7 

7.6 

16.6 

400 

28 

10 

623 

1.16 

11  6 

33.3 

7.2 

23.8 

300 

47 

19 

555 

l.l4 

21.7 

55.0 

12.0 

35.8 

200 

95 

48 

443 

l.ll 

53-5 

108.5 

23.6 

59.4 

During  the  interval  At*  the  shock  front  progresses  a  distance 


Ax.  .  =  c.  .  At*.  . 
ij  10  10 

Therefore  a  graph  of  the  upper  and  lower  shock  stress  values  can  be  plotted 
against  the  distance  t  ie  shock  has  penetrated.  This  is  done  in  Figure  4.21. 


Step  6 

Draw  the  waveforms  corresponding  to  the  time  t  *  =  h  ,/c  (the  time  the 
wave  front  first  reaches  tne  structure)  and  x  =  h,  and  several  otuer  inter¬ 
mediate  waveforms.  This  is  shown  in  Figure  4.22 


Step  7 

Find  the  slope  of  the  instantaneous  stress  waveform  at  the  depth  h,  for 
each  time,  t*,  and  plot  this  slope  against  t*. 


7.8  ms  (the  time  wavefront  first  reaches  structure) 


=  15  ms 


Figure  4.22  WAVEFORMS  AT  DIFFERENT  VALUES  OF  t* 
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velocity  immediately  before  the  shock 
velocity  immediately  after  the  shock 
mass  density  (assumed  constant) 
stress  immediately  before  the  shock 
stress  immediately  after  the  shock 


strain  immediately  before  the  shock 
strain  immediately  after  the  shock 


During  the  unloading  portion,  beginnirg  immediately  after  the  shock,  elastic 


theory  can  be  used 


,  so 

^t  =  jk-  <5t>x  =  (^) 
1  0 


X 


.106 


where 


velocity  during  the  unloading  phase 
velocity  immediately  after  the  shock 
stress  immediately  after  the  shock 


stress  during  unloading 
density  (assumed  constant) 

velocity  of  propagation,  computed  from  the  initial  tangent 
modulus  and  assumed  constant,  thereby  implying  linear  elasticity 

4.1.5  Vertical  and  Horizontal  Free  Field  Inputs 

Strictly  speaking  the  vertical  velocity  would  be  obtained  by  multiplying 
the  velocity  curves  developed  for  Cases  I  and  II  of  Section  4.1.4  by  the  cosine 
of  the  Mach  angle.  The  Mach  angles  are  small  however  and  this  refinement 
therefore  is  smaller  than  the  errors  already  included  in  the  curves. 

The  horizontal  motion  components  are  more  difficult  to  estimate.  It  has 
been  pointed  out  already  that  test  records  diverge  grossly  from  the  predictions 
of  elastic  theory;  the  shear  wave  does  not  terminate  horizontal  velocity  as 
predicted  by  elastic  theory. 

If  suspension  systems  are  used  we  have  found  that  the  horizontal  systems 
should  be  of  low  frequency.  For  this  case  the  only  significant  input  is  the 
total  horizontal  displacement. 

If  shock  isolation  is  not  used,  dependence  being  placed  on  equipment  rugged¬ 
ness  and  the  inherent  "isolation"  provided  by  soil  structure  interaction  only  the 
initial  portion  of  the  horizontal  velocity  curve  is  required.  If  ruggedized 
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equipment  can  resist  the  onset  of  the  ground  shock  wave  It  can  subsequently 
resist  the  onset  of  the  reflected  wave. 


The  initial  horizontal  velocity  can  be  obtained  by  multiplying  the  predicted 
vertical  velocity  by  the  sines  of  the  appropriate  Mach  angles.  Since  the  Mach 
angles  are  small  at  the  higher  pressures  the  resulting  velocities  will  be  small. 

At  some  pressure,  not  yet  known,  the  directly  transmitted  shock  will  dominate 
the  horizontal  component  of  air  induced  motion. 

The  total  horizontal  displacement  is  equal  to  the  area  under  the  horizontal 
velocity  curve.  However,  it  must  be  remembered  that  the  horizontal  component  of 
reflected  wave  velocity  adds  to  the  direct  component.  In  a  purely  elastic  medium 

Wo.  uii  perrce^  ruxlcvwiun  at  DcCirtTClv  ctzxtL  HO  o  1  gfii X  Icaub  wave  dSCaiy  Xnc  liOirxZQEtucLL 

displacement,  x^,  would  be  given  by 

*h  -  *  sin  $ 

where 

I  is  the  impulse  of  the  aeroshock  wave, 

0  is  the  Mach  angle  of  the  ground  shock  compression  wave. 

It  does  not  seem  intuitively  reasonable  though  to  consider  that  after 
vertical  motion  of  the  structure  reaches  its  maximum  (at  arrival  of  the  reflected 
wave)  horizontal  motion  continues  for  any  length  of  time.  For  this  motion  to 
continue,  the  column  of  soil  between  structure  and  bedrock  would  have  to  continue 
to  shear,  a  process  that  rapidly  dissipates  energy  in  real  soils. 

Therefore  we  take  as  a  reasonable  approximation  of  peak  horizontal  displacement 

z,  =  sin  0  z 
h  max  v  max 


4ft 
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where  z,  and  z  are  the  maximum  components  of  horizontal  and  vertical 

h  max  v  max 

motions  respectively.  This  formula  is  equivalent  to 

_  c 

^h  max  “  U  v  max 

where 

U  is  the  aero shock  velocity 
and  c  is  an  average  wave  velocity. 

A  value  of  c  can  be  computed  from 


This  approximation  yields  a  horizontal  to  vertical  displacement  ratio  on 
the  order  of  1/3  for  Nevada  Test  Site  soil  and  overpressures  in  the  100  to  300 
psi  region.  --  ■  -  . . .  --  - 

4.2  Approximate  Estimates  of  Input 

If  preliminary  estimates  of  rattle  space  requirements  must  be  carried  out 
without  the  benefit  of  site  information  it  is  suggested  that  the  following  approach 
be  used. 

1.  Estimate  bedrock  depth  from  topographic  maps.  If  the  site  is 
unknown  guess  bedrock  to  be  150  feet  below  surface. 

2.  If  the  general  class  of  soil  (clay,  sand)  is  known  estimate 
compression  wave  velocity  (see  Section  8).  If  the  soil  is 
unknown  guess  a  wave  propagation  velocity  of  1000  ft/sec. 

3.  Estimate  times  required  for  wave  to  travel  from  surface  to 
structure  and  surface  to  bedrock  to  structure. 
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4.  Estimate  wave  decay  according  to  percyclic  damping  theory.  A 
conservative  value  of  the  damping  coefficient,  6,  can  he  taken 
as  0.2.  (it  will  be  r.igher  for  most  real  soils.) 

5.  Estimate  the  maximum  impulse,  I,  that  can  be  encompassed  between 
structure  and  bedrock.  This  will  occur  at  about  the  time  the 
reflected  wave  reaches  the  structure. 

6.  Estimate  maximum  vertical  displacement,  z  ,  from  the  formula 

v  max 

I 

z  =  — 

v  max  pc 


Estimate  maximum  horizontal  displacement  z,  ,  from  the  formula 

n  HL8X 

=  i 

^  max  pU 


where  U  is  the  aeroshock  front  velocity. 


Silos  are  sometimes  founded  on  bedrock.  If  the  silo  shell  does  not  contain 
girth  expansion  joints  the  silo  vertical  motion  will  be  equal  to  bedrock  vertical 
motion.  Hie  latter  can  be  estimated  conservatively  as  follows: 


1.  Estimate  bedrock  seismic  velocity,  c  .  If  bedrock  is  unknown 

s 

guess  it  to  be  10,000  ft  per  second. 

2.  Estimate  bedrock  density.  If  bedrock  is  unknown  guess  the 
density  to  be  150  lbs  per  cubic  foot. 

3.  Estimate  the  impulse  of  the  overpressure  wave,  I. 

4.  Compute  vertical  displacement 

_  1  I_ 

Zv  max  2  Pc 

s 

The  factor  ^  is  included  to  account  roughly  for  geometric  dispersion  of  +he 
wave.  Its  actual  value  generally  will  be  more  nearly  1/3  or  l/4. 
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SECTION  5 


SOIL  -  STRUCTURE  INTERACTION 


A  simplified  analysis  of  soil- structure  interaction  is  given  in  Appendix  A. 
Even  the  simplified  analysis,  however,  is  quite  complex  to  apply.  In  this 
section  the  main  results  of  the  analysis  and  estimation  procedures  accounting 
for  the  principal  influences  of  the  interaction  will  y  »  given. 

The  most  important  prediction  of  the  interaction  study  is  that  the  significant 
modes  of  most  structures  in  real  soils  will  be  damped  near  critical;  the  lower 
modes  will  be  over  critically  damped.  This  then  implies  that  the  structure  will 
follow  the  soil  motion  quite  closely  except  within  an  interval  equal  approximately 
to  the  first  half  or  three  quarters  period  of  undamped  vibration  of  the  structure. 

A  second  consequence  of  the  interaction  study  is  that  the  rise  time  for  the 
forcing  functions  of  the  lower  modes  (translational  rigid  body  mode  and  first 
bending  mode)  is  on  the  order  of  3/4  to  1  t  mes  the  transit  time  of  the  ground 
shock  wave  over  the  structure.  Further,  it  can  be  deduced  by  inspection  of 
Equations  (A-19)  and  (A-24)  (Appendix  A)  that  the  initial  rates  of  rise  of  the 
forcing  functions  for  the  higher  order  modes  are  no  greater  and  generally  less 
than  the  rate  of  rise  of  the  first  sending  mode  and  rigid  body  mode  (the  initial 
rates  of  rise  of  these  are  identical).  Very  shortly  after  initial  impingement 
of  the  ground  shock  wave  on  the  structure  the  rates  of  rise  of  the  higher  order 
forcing  functions  become  less  than  the  rates  of  rise  of  the  rigid  body  and  first 
bending  mode  forcing  function.  Further,  the  maximum  amplitudes  of  the  higher 
order  forcing  functions  decrease  rapidly  with  increase  of  modal  index.  (A 


review  of  Appendix  A. 4  with  attention  directed  toward  the  physical  meaning  of  the 
quantities  involved  will  make  this  evident  without  mathematical  manipulation. ) 

The  response  of  the  structure  modes  to  their  forcing  function  also  decreases 
rapidly  with  increase  of  modal  index.  Thus  we  intuitively  expect  that  the 
significant  motion  of  the  structure  interior  will  depend  only  upon  the  rigid 
body  mddes  and  the  first  few  distortion  modes. 

It  becomes  a  practical  impossibility  to  draw  general  conclusions  from  purely 
mathematical  deductions  because  of  the  complexity  of  the  problem.  Accordingly, 
a  computer  program  was  set  up  and  the  relative  responses  of  the  various  modes  for 
a  silo-like  structure  were  determined  for  a  range  of  parameter  variations.  (The 
computer  program  and  detailed  results  obtained  are  included  in  the  first  interim 
report  on  this  project.  They  are  not  repeated  in  this  report.) 

The  results  of  this  investigation  indicate  that  only  two  modes  need  be 
considered  for  developable  structures: 

1.  The  rigid  body  mode. 

2.  The  first  deformational  mode. 

An  elementary  manipulation  of  the  formulas  for  the  soil  resistance  factor 
K  and  the  modal  damping  factor  6^  shows  that  for  reinforced  concrete  shells  the 
damping  factor  for  the  rigid  body  mode,  6  ,  lies  between  the  limits 

0.6yf .  ar  <  0.T5yi  (5.D 

R/h  is  the  radius  to  thickness  ratio  of  the  shell.  (For  steel  shells  the  limits 
of  the  coefficient  ait  O.33  to  0.4?.) 


Thus  even  for  a  truly  massive  concrete  shell  having  an  R/h  ratio  of  only  4 


the  rigid  body  mode  is  overcritically  damped. 

If  a  structure  is  so  stiff  that  the  first  defonnatioual  mode  is  less  than 
critically  damped  the  deformation  of  the  mode  is  insignificant  when  compared 
with  the  rigid  body  motion  and  may  safely  be  neglected. 

5.1  Summary  of  Interaction  Theory 

The  theory  developed  in  Appendix  A  will  be  summarized.  Since  this  summary 
will  be  concerned  only  with  the  rigid  body  mode  and  first  deformational  mode  the 
subscripts  r  and  d  respectively  will  be  used  to  denote  them.  A  subscript  "i" 
will  be  used  to  denote  a  mode  in  general. 


The  amplitude,  T. ,  of  the  response  of  the  ith  mode  is  given  by  the  equation 


T.  +  25, ft.  T.  +  ft.  T.  =  t. 
1  iix  11  1 


(5-2) 


where 


Is  the  forcing  function  for  the  ith  mode,  to  be  defined  subsequently 


ft^  is  the  circular  frequency  of  the  ith  mode 
and  is  the  damping  factor  for  the  ith  mode. 


For  the  rigid  body  mode  the  circular  frequency  ft  x  given  by 


ft 


„  K 
2  r 


(5.3) 


and 


r  _  Pc 

Kr  '  n  wr 


(5=4) 


where 


m 

c 

P 

x> 


=  unit  mass  of  the  shell 

=  soil  compression  wave  velocity 

=  soil  unit  mass 

=  Poisson’s  Ratio  of  soil 

X^  =  a  characteristic  length  of  the  shell. 

The  K  factor  is  difficult  to  estimate.  Though  the  simplified  equations 
introduce  a  constant  to  account  for  the  soil  displacement  resistance  it  is  in 
fact,  even  for  a  linear-,  elastic  medium,  a  variable  depending  upon  the  medium 
characteristics,  structure  geometry,  and  manner  of  loading  the  structure.  For 
a  spherical  shell,  which  can  be  investigated  analytically,  the  K  obtained  for 
uniform  radial  expansion  of  the  sphere  is  twice  that  obtained  for  a  rigid  body 
translation  of  the  sphere. 


The  average  K  factor  computed  for  the  first  deformational  mode  lies  between 
these  limiting  cases. 


It  is  suggested  that  in  actual  computation  two  K  factors  be  used,  even 
though  the  derivation  of  the  basic  equations  is  based  on  a  single  constant  value. 

Equation  (5.^-)  is  based  on  the  theory  of  elasticity  so  it  can,  at  best,  be 
only  approximate  for  real,  nonlinear  soils.  The  characteristic  length  X  is 
given  approximately  by  the  equation 


where 

V 

A 


X 

r 


3V 

A 


structure  volume 
structure  surface  area. 


(5.5) 
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Equation  (5.5)  is  accurate  within  the  approximations  of  the  theory  of 
elasticity  for  spherical  shells  and  any  other  closed  shells  so  proportioned  that 
a  small,  rigid  lateral  displacement  of  the  shell  induces  a  uniform  stress  in  the 
direction  of  displacement  at  the  soil-shell  interface. 

If  experience  computing  static  foundation  displacements  can  be  considered 

indicative  of  the  accuracy  of  Equations  (5.3)#  (5-*0#  and  (5-5)  then  they  must 

+ 

be  considered  to  have  an  overall  uncertainty  factor  of  -  50  percent  or  so. 
However,  we  simply  have  not  enough  experimental  evidence  available  yet  to  justify 
development  of  more  precise  formulas. 

The  K,  factor  for  the  first  deformational  mode  can  be  taken  50  percent 
d 


higher  than  Kr* 


K.  =  1.5  K 
d  r 


(5.6) 


For  the  bending  mode 


fi2  =  CD2  +  J* 
d  dm 


(5.7) 


where 


tb  =  circular  frequency  of  first  bending  mode  in  vacuum 


(See  Section  8) 


The  damping  factor  is  defined  for  both  modes 


5  =—  c 

i  2  mj^ 


(5.8) 


(A. 24) 


The  forcing  function,  t  is  given  by 


the  norm,  N^,  is  given  by 

N,  dA,  (A. 25) 

1  A  1 

and  the  input  vector,  p  ,  by 

p.  =  n  •  p  +  Kz  +  pc  z.  (A. 19) 

In  these  formulas 
A  =  structure  area 

-  vector  mode 
n  =  normal  unit  vector 

p  =  free  field  pressure  tensor 

z  =  free  field  particle  displacement. 

Evidently,  any  system  of  solution  that  retains  even  the  major  influences 
of  these  equations  is  bound  to  be  complex.  It  is  possible  though  to  deduce  the 
principal  influences  of  the  interaction  from  critical  examination  of  the  equations 
and  to  frame  a  much  simpler  method  of  solution  that  is  accurate  enough  for  most 
shock  isolation  computations. 

5.2  Simplification  of  Interaction  Procedure 

We  begin  by  considering  the  input  vector  p^. 

Appendix  Section  A. 4  is  devoted  to  computation  of  p^.  Following  the  method 
presented  there  wo  envision  the  structure  to  be  enclosed  by  an  imaginary  rectangular 
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parallelepiped  having  two  faces  parallel  to  the  "plane"  of  the  advancing  ground 


shock,  wave.  Then  the  vector  p^  can  be  represented  as  six  scalar  forces  perpendi¬ 
cular  to  the  faces  of  the  parallelepiped. 

For  the  rigid  body  motion  the  forces  on  the  four  surfaces  that  are  not 
parallel  to  the  plane  of  the  shock  wave  contribute  nothing  to  the  forcing  function. 
On  the  two  remaining  surfaces  the  scalar  pressures  p  can  be  simplified  to: 

For  the  windward  surface 

P1  =  P  +  pcz  +  Kyz  (5.9) 

and  on  the  leeward  surface 

Px  =  p  -  pcz  -  Krz  (5.10) 

where  p  is  the  free  field  wave  pressure  and  z  is  the  free  field  velocity.  The 
free  field  displacement  z  can  he  obtained  by  integration  of  the  z  curve. 

At  any  instant  pcz  is  equal  to  p.  However,  in  order  to  carry  out  the 
interaction  computations  an  average  value  of  pc  must  be  chosen.  Since  the  shapes 
of  the  p  and  z  curves  may  differ  markedly  the  distinction  between  the  terms  has 
been  retained. 

For  the  deformational  mode  the  scalar  pressures  on  the  four  parallelepiped 
surfaces  perpendicular  to  the  ground  shock  wave  are 

P-L  =  -  PP  (5-11) 

where 

u  is  Poisson's  Ratio. 
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Now  let  us  examine  the  result  of  evaluation  of  the  forcing  functions  T 

and  The  vector  mode  0^  consists  simply  of  a  unit  vector  in  the  direction 

of  the  advancing  ground  shock  wave.  Therefore  its  scalar  projection  on  the 
windward  face  of  the  parallelepiped  is  simply  +1  and  on  the  leeward  face  is  -1. 
Thus  if  the  transit  time  of  the  wave  moving  over  the  structure  were  small  and 

the  rate  of  decay  of  p^  not  too  great,  after  engulfment  the  forcing  function 


would  he  given  to  a  very  close  approximation  by 


■  (pcz  +  K  z) 
r  N  r  ' 

r 


(5-12) 


For  the  rigid  body  mode 


N  2m 
r 


T  -  i  (pc£  +  K  z) 
r  m  r 


(5.13) 


(5.14) 


We  note  also  that  since  is  an  integral  over  the  surface  of  the  structure, 
peaks  and  steps  will  not  appear  in  T  The  abrupt  onset  of  pressure  cannot  be 
applied  simultaneously  to  the  entire  surface  of  the  structure  or  even  to  a  major 
fraction  of  the  surface.  The  effect  of  the  integral  is  to  smooth  out  irregularities 
within  an  interval  equal  to  the  transit  time. 

Since  the  component  of  contributed  by  p  cancels  out  at  any  time  when  the  p's 
on  opposite  sides  of  the  structure  are  equal,  the  only  contributions  of  p  to  the 

forcing  function  occur  when  there  is  a  rapid  change  in  the  magnitude  of  p  within 

an  interval  equal  to  the  transit  time.  If  such  rapid  changes  exist  they  are 

averaged  over  the  transit  time  interval  and  make  small  contributions  to  r  . 
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At  this  point  it  is  of  interest  to  compare  the  transit  time  with  the  period 
of  vibration  of  the  rigid  body  mode.  For  a  horizontal  cylinder  the  transit  time 


is  given  by 


t 

t  c 


The  period  of  vibration,  t  ,  of  the  rigid  body  mode  is 


"r  -  27TJr 


(5-15) 


(5.16) 


For  a  cylinder,  this  is  equivalent  to 


»  27®  (Jh1/2  3  _£ 

r  c  'R'  2  v  '  7 

'b 


(5-17) 


Thus  the  ratio  t  /t.  is 
r  t 


where 


(5.18) 


shell  thickness 
shell  radius 
shell  density 
soil  density 
Poisson's  Ratio  of  soil. 


’’•'or  a  concrete  shell  and  a  Poisson's  Ratio  of  0.3 


t  , 


hvl/2 


(5-19) 


Thus  for  a  cylindrical  shell  having  a  radius  to  thickness  ratio  less  than 

25  the  transient  time  is  less  than  a  period  of  vibration  in  the  rigid  body  mode. 

This  mode  is  overcriticat ly  damped.  Therefore  variations  in  the  t  curve  will 

r 

not  strongly  influence  the  resulting  response  curve.  Accordingly,  the  influence 


<  • 


of  p  on  the  rigid  body  mode  can  be  neglected  and  only  the  z  and  i  components 
considered. 


In  accordance  with  the  approximation  already  adopted  (Equation  (5*6)) 


Equation  (5*7)  can  be  rewritten 


2  2  2 
fl*  =  (l)  j  +1-5  ft 


(5.20) 


Substituting  the  expression  for  (Section  8)  and 

rv2  3  /4i\2  /  C\2  1  /R\ /C\2 

<E">  +  Ii«T  (h)(R) 


(5.21) 


where 


7g  =  soil  density 

yc  =  shell  density 

cc  =  sonic  velocity  of  shell  material 

and  other  symbols  are  as  previously  defined.  Then  the  ratio  t,/t  becomes 


*d)1/2 


1  2.  /  C\2/h..3 

(l+V>)  (7^  5  V  'R' 


(5.22) 


Taking  for  a  concrete  shell 
cc  =  11,000  ft/sec 
7  =  150  lbs/ft3 


and  for  soil 


c  =  1500  ft/sec 

7  -  115  lbs/ft3 


x>  -  0.3 
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(5.24) 


The  damping  factor  can,  after  some  manipulation,  he  expressed 


These  have  the  values 


R/h  - 

4 

7 

10 

20 

tA  ■ 

0.22 

0.38 

0.51 

0.71 

• 

0.118 

0.355 

0.68 

1.88 

Examination  of  the  table  shows  that  for  some  possible  values  of(R/h)  the 
damping  of  the  first  vibrational  mode  is  considerably  less  than  critical. 

However,  we  note  also  that  for  these  cases  the  period  of  vibration  is  only  a 
fraction  of  the  transit  time  and  therefore  the  response  will  closely  follow  the 
input . 

The  vector  mode  is  symmetrical.  Therefore,  the  velocity  and  displacement 
terms  cancel  out  except  when  these  have  a  rapid  change  of  value  within  the  interval 
of  the  transit  time.  Even  these  rapid  changes  will  not  markedly  influence  the 
resulting  response  curves. 


Taking  p  to  be  0.3  the  forcing  function  r ^  becomes  after  engulfment 

(5.25) 


~  0.27  p 

7 -  -  - - 

d  m 


During  the  initial  half  of  the  engulfment  time  all  three  components  will 
contribute  to  both  forcing  functions.  Therefore,  if  the  initial  rise  proves  to 
be  significant  it  is  suggested  that,  the  amplitudes  of  and  7 ^  be  computed  at 


1/2  t  reflecting  the  three  components.  Smooth  curves  starting  at  zero,  passing 
through  the  computed  points  at  l/2  t^  and  smoothly  merging  with  the  values  ac  t^ 
then  can  be  sketched. 


It  is  noted  that  if  Equation  (5.2)  is  divided  through  by  Q  then  the 
term  ^  has  the  dimensions  of  length  and  at  some  distance  from  the  origin 

will  represent  the  terminal  displacement  that  the  damped  response  asymptotically 
approaches.  In  this  form  the  relative  importance  of  the  rigid  body  and  deformational 
modes  can  be  easily  judged. 

Then,  subsequent  to  engulfment,  the  equation  for  the  rigid  body  mode 
forcing  functions  becomes 


i  3(l+h)  v  i  + 

tt2r  c  A 


(5*26) 


where 


structure  volume 


structure  surface  area 


Poisson's  ratio 


average  compression  wave  velocity 
free  field  particle  velocity. 


Now  the  expression 


3(l+o)  J 


is  of  the  same  order  of  magnitude  as  the  structure  "radius"  regardless  of 
structure  complexity.  Thus  the  first  teru  of  Equation  (5-26)  is  approximately 


equal  to 
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From  a  typical  example  it  can  be  shown  that  the  maximum  contribution  of 
the  velocity  term  to  the  displacement  is  on  the  order  of  1  inch  for  a  30  ft 


diameter  cylinder  and  1500  ft  per  second  soil.  The  maximum  displacement  is  about 
18  inches  for  the  same  structure;  thus  the  velocity  term  contributes  a  displace¬ 
ment  increment  only  5  percent  as  great  as  the  soil  displacement  term,  z. 

The  deformational  motion  is  more  diff-'  cult  to  investigate  because  the 
modal  geometry  0^  varies  between  structures.  For  a  horizontal  cylinder  the 
equation  becomes 

-§  - - f  B - -^j—  (5.2T) 

*  7c  52  [|  |)3  +  <&>*#] 

For  a  reinforced  concrete  shexl  of  R/h  =  7.5  and  c  =  1500  ft/sec,  the  term  in 
brackets  becomes  about  0.66.  Therefore 


for  a  30  ft  diameter  shell. 


(5.28) 


This  is  a  small  displacement  in  comparison  with  the  rigid  body  displacement. 
Another  fact  emerges,  however,  that  is  of  practical  importance.  The  displacement 
indicated  by  Equation  (5.27)  is  almost  exactly  the  displacement  suffered  by  the 
structure  surrounded  by  the  soil  and  statically  loaded. 

In  Section  5.3  a  simplified  design  method  based  on  these  observations  is 


1  ,'U 


given 


5*3  Simplified  Design  Method 


It  is  again  pointed  out  that  this  design  procedure  is  necessary  only  if 
isolators  are  to  be  omitted  or  if  the  more  desirable  and  economical  low  frequency 
systems  are  not  used.  A  method  for  determining  when  isolators  can  be  omitted  is 
given  in  Section  7* 


The  interior  structure  motion  consists  of  two  components: 

1.  A  rigid  body  motion  Tr 

2.  A  deformational  motion  T^ 


The  rigid  body  motion  is  the  solution  of  Equation  (5-29) 


T  +  T 


r  r 


z 


where 


and 

z 

m 

V 

A 

P 

c 


_2  _  Pc2 A 

r  IH^TinV 


_££_ 

2nfir 


26r  (l+v>)3V 

_n_  =  — - 

r  c  A 


free  field  displacement 
structure  mass  per  unit  area 
structure  volume 
structure  surface  area 

soil  unit  mass 

average  compression  wave  velocity. 


(5.29) 


If  t  is  the  transit  time  of  the  ground  shock  wave  over  the  structure  and 


27T 

trmTT 

r 


(5-30) 


the  response  Tr  will  be  practically  identical  with  z  after  a  time  period  t* 
where 


t*  **  t  +  t . 
r  t 


(5.31) 


If  necessary  Tr  can  be  determined  within  the  interval 


0  <  t  <  t# 

by  solution  of  Equation  (5-29) 


The  deformational  motion  can  be  approximated  by  Equation  (5.32) 

«  • 

T 

4-  +  (l+u)  P  T,  +  T  =1  (5-32) 

S?  CA  d  d 

r 

where 

n2d  -  a>2a+1.5«2r  (5-33) 

O)^  =  circular  frequency  of  structure  vibration  in  vacuum 

A 

z  =  static  deformation  of  structure  in  soil  under  the 
imposed  pressure. 

A 

Foe  computing  z  the  actual  soil  p  ssure  at  the  structure  depth  should  be 
averaged  out  within  an  interval  equal  to  the  transit  time  t  . 

The  net  interior  structure  motion,  T,  is  then 

T  =  Tr  +  Td  (5  .34  ) 

One  of  many  methods  for  solving  Equations  (5*29)  and  (5-32)  is  given  in 
Appendix  B. 
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In  this  section  the  pressures  and  velocities  used  have  been  the  net  values 
acting  in  the  direction  of  soil  particle  motion.  These  are  the  appropriate 
inputs  for  "floating”  structures  completely  surrounded  by  soil. 

If  a  structure  founded  on  bedrock  is  used,  e.g.,  some  silos  already  built, 
only  the  horizontal  component  of  motion  need  be  considered.  The  equations  of 
this  section  are  valid  for  the  horizontal  component  if  horizontal  pressures  and 
velocities  are  used  as  inputs  and  if  substantial  vertical  motion  is  prevented 
by  the  bedrock. 


127 


SECTION  6 


NONLINEAR  SYSTEMS 


Nonlinearities  of  shock  isolation  system  must  be  considered,  from  two 
standpoints . 

1.  All  real  systems  contain  nonlinear  elements  even  though  for 
practical  design  purposes  we  may  simplify  the  analysis  by 
"linearizing"  the  systems.  In  most  instances  the  errors  so 
introduced  are  of  no  practical  importance  and,  therefore, 

the  procedure  is  perfectly  acceptable.  However,  occasionally, 
the  nonlinearitiec  introduce  qualitative  changes  in  the 
behavior  of  the  real  system  which  are  of  sufficient  importance 
that  they  must  be  considered  in  the  design. 

2.  For  various  reasons  it  may  be  desirable  to  deliberately  design 
nonlinear  isolation  systems  in  order  to  gain  specific  advantages. 

The  theory  of  nonlinear  vibrations  is  extremely  complex  and  still  in  its 
infancy.  Further,  by  far  the  largest  portion  of  the  work  already  carried  out 
understandably  has  been  directed  towards  solution  of  steady-state  problems  since 
these  comprise  the  largest  proportion  of  nonlinear  problems  of  practical  importance. 
Because  the  principle  of  superposition  cannot  be  used  for  the  solution  of  nonlinear 
mechanical  problems  the  variety  of  existing  steady-state  solutions  are  of  very 
limited  value  to  designers  interested  in  the  transient  response  of  nonlinear 
systems.  No  analog  of  the  Duhamel  integral  exists  for  nonlinear  vibration 
problems . 
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At  the  present  state  of  our  knowledge  no  general,  well  developed  theory 
of  transient  nonlinear  vibration  exists .  The  mass  of  theory  and  data  available 
is  directed  toward  solution  of  certain  specific  problems  or  of  problems  that 
can  be  characterized  by  certain  specific  nonlinear  differential  equations. 

Further,  even  at  best,  application  to  particular  problems  of  the  methods  that 
have  been  developed  is  generally  quite  tedious  and  time-consuming  and,  if 
carried  out  analytically,  almost  invariably  requires  the  use  of  a  high  speed 
computer. 

It  was  believed  that  no  purpose  would  be  served  by  burdening  this  report 
with  a  compendium  of  solutions  and  techniques  for  analyzing  specific  nonlinear 
problems.  Determination  of  the  effects  of  nonlinearities  on  nearly  linear 
systems  is  appropriate  to  the  refined  final  design  of  shock  isolation  systems 
which  is  not  within  the  scope  of  this  report.  However,  preliminary  design  of 
nonlinear  systems  that  might  be  used  in  lieu  of  the  more  conventional  linear 
types  is  within  the  scope  of  this  report.  Therefore,  on  the  project  some  effort 
was  expended  to  develop  criteria  that  would  point  out  the  advantages  or  dis¬ 
advantages  of  nonlinear  systems  under  particular  circumstances  and  to  develop 
methods  to  expeditiously  determine  rattle  space  requirements,  acceleration  limits, 
etc .  for  nonlinear  systems .  An  attempt  was  made  to  do  this  for  the  general  non¬ 
linear  system  as  opposed  to  cataloging  specific  types  of  nonlinearities. 

Obviously,  the  results  presented  are  far  from  complete  even  for  cases  to  which 
they  most  strictly  apply.  Again  we  emphasize  that  the  purpose  of  this  section 
is  to  afford  the  designer  fairly  rapid  methods  that  he  can  use  to  expeditiously 
establish  limits  for  the  important  parameters  of  rattle  space  requirements  and 


peak  accelerations  to  which  the  isolated  body  will  be  subjected.  Even  more 
important,  perhaps,  is  presentation  of  means  by  which  he  can  determine  whether 
a  nonlinear  system  has  any  clear  cut  advantages  or  disadvantages  over  a  more 
conventional  linear  system. 

Nonlinear  systems  can  be  characterized  in  a  number  of  ways.  In  this 
section  we  will  find  it  convenient  to  regard  the  isolator  as  an  energy  storing 
device.  Two  general  classes  are  considered: 

1.  Conservative  systems  in  which  all  the  energy  delivered  to  the 
isolator  by  mechanical  motion  and  stored  within  the  isolator  as 
potential  energy  ultimately  is  released  by  the  isolator  as 
mechanical  motion.  A  conservative  system  is  the  nonlinear 
analogue  of  undamped  linear  isolators. 

2.  Nonconservative  systems  within  which  a  portion  of  the  energy 
delivered  to  the  isolator  by  mechanical  motion  is  dissipated 

as  heat  or  otherwise  and  is  not  subsequently  released  as  mechanical 
motion,  the  analogue  of  damped  linear  isolators. 

6.1  Conservati-vs  Nonlinear  Systems 

Consider  a  single  degree  of  freedom  nonlinear  isolator.  J.t  can  be  symbolized 
by  the  scheme  shown  in  Figure  6.1.  The  problem  is  to  determine  the  displacement 
x(t )  and  the  force  F(t)  in  terms  of  the  known  input  z(t)  and  the  isolator 
characteristics.  A  variable  y(t),  corresponding  to  the  rattle  space,  is  defined 

y(t)  ---  x(t)  -  z(t)  (6.1) 

Note  that  y(t)  is  defined,  positive  for  a  decrease  in  the  space  between  the  mass 
and  Lke  isolator  support. 

.1  v> 


Figure  6.1  SCHEMATIC  OF  CONSERVATIVE  NONLINEAR  SYSTEM 


b.3  TOROE-STnOKE  DIAGRAM  OF  ISOLATOR 


■  ^  I  *  1  s  f  I  •  •  f  v  f 


Now,  regardless  of  the  mechanical  complexity  of  the  isolator,  if  it  is 
conservative,  its  mechanical  characteristics  can  be  represented  by  a  curve 
similar  to  Figure  6.2  if  the  mass  of  the  isolator  is  taken  to  be  zero.  If  the 
masses  of  the  isolator  elements  are  considered  the  curve  of  Figure  6.2  broadens 
somewhat,  the  amount  of  broadening  being  proportional  to  the  ratio  of  the  kinetic 
energy  of  the  masses  to  the  potential  energy  stored.  This  condition  is  roughly 
analogous  to  a  spring -mass  system  having  a  spring  so  heavy  its  mass  must  be 
included  in  the  analysis.  Even  for  the  linear  system  this  greatly  complicates 
analysis.  Therefore,  in  this  section  -the  mass  of  the  isolator  will  be  presumed 
small  enough  to  neglect. 

For  the  massless  isolator  the  force  developed  is  equal  to  the  derivative 
of  potential  energy,  i.e., 

dE 

(6-2) 

Assume  now  that  the  isolator  has  been  placed  on  its  support.  Its  compression 

y  and  potential  energy  E  are  zero.  Then  as  the  supported  mass  is  gradually 

P 

lowered  on  the  isolator  it  compresses  an  amount  y,  stores  potential  energy  E^, 
and  develops  a  force  F. 

At  equilibrium  the  force  F  must  be  equal  to  the  weight  mg. 

Figure  6.4  illustrates  two  differing  energy- displacement  curves  and  their 
resulting  equilibrium  positions.  A  practical  value  of  nonlinear  systems  is 
immediately  evident.  The  equilibrium  energy  of  system  2  is  much  smaller  than 
that  of  system  1  for  the  mass  chosen.  Now  as  a  first  approximation  the  cost  of 
an  isolation  system  is  roughly  proportional  to  the  equilibrium  energy.  Therefore 
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Figure  fi.4  COMPARISON  OF  E^JT.UKBIUM  ENERGY  OF  TWO  NOHIOTEAR  SYSTEMS 


if  system  2  were  of  the  same  order  of  mechanical  c onplexity  as  syrten  1,  it 
would  be  considerably  cheaper  than  system  1. 


Ve  also  note  that  if  the  two  systems  were  designed  so  that 


d? 

dy 


(6.3) 


were  nearly  equal  and  constant  within  the  stroke  range  R  centered  on  the 
equilibrium  position  the  systems  would  be  nearly  linear  within  R  and  of  the 
same  frequency.  A  bi -linear  system  such  as  system  2  could  be  designed  for  a 
protective  structure  by  making  R  equal  to  twice  the  relative  displacement 
taken  from  a  shock  spectrum  at  frequency 


(6A) 


A  close  approximation  of  the  rattle  ^pace  requirements  can  be  quickly  made 
using  the  Ep  versus  y  and  F  versus  y  graphs. 

The  time  required  for  the  vertical  component  of  structure  motion  to  reach 
its  maximum  is  about  equal  to  the  time  required  for  the  ground  pressure  wave  to 
travel  from  the  structure  to  bedrock  and  reflect  back  to  the  structure .  Now, 
generally,  the  structure  will  be  within  100  feet  of  bedrock  and  the  ground 
pressure  wave  propagates  at  500  to  1000  ft/sec.  minimum.  Thus  the  time  to 
maximum  downward  displacement  of  the  structure  is  on  the  order  of  250  msec,  or 
less.  The  magnitude  of  the  structure  displacement  is  given  with  fair  accuracy 
by 

z  =  (6.5) 

max  pc 


c 


(6.6) 


where 

l(t*)  -  impulse  of  blast  wave  (ground  wave)  between 

time  zero  and  time  t* 

D  *  Depth  of  bedrock  below  structure 

TT  *  Soil  compression  wave  propagation  velocity 

p  *  soil  mass  density 

As  soon  as  the  floor  of  the  structure  begins  to  descend  the  isolator  begins 
to  expand  and  deliver  some  of  its  stored  energy  to  the  floor  and  soil  below. 

In  expanding  the  force,  F,  exerted  against  the  mass  decreases  somewhat  so  the 
mass  begins  to  descend,  thereby  gaining  kinetic  energy  due  to  its  velocity  and 
losing  potential  energy  as  measured  in  a  motionless  coordinate  system.  However, 
the  potential  energy  of  the  mass  with  respect  to  the  floor  increases  because 
the  floor  descends  faster  than  the  mass  and  therefore  the  mass  rises  with 
respect  to  the  floor. 

Since  isolator  systems  generally  are  of  long  period  (a  second  or  so)  and 
maximum  vertical  displacement  of  the  floor  occurs  generally  in  less  than  l/4 
second,  as  a  first  approximation  we  can  assume  that  the  mass  remains  motionless 
during  floor  descent  to  maximum.  Then,  if  necessary  we  can  correct  for  initial 
mass  motion. 

If  the  mass  remains  motionless  during  floor  descent  a  distance  Az  the 
isolator  opens  an  amount  -  /5y  and  the  potential  energy  of  the  mass,  with  respect 
to  the  floor,  increases  an  amount  mg^y. 


Figure  6.5  shows  normalized  Ep  versus  y  and  F  versus  y  curves  for  a  hypotheti¬ 
cal  system.  The  curves  have  been  normalized  by  dividing  both  E^  and  F  by  the 
weight  of  the  isolated  mass .  The  rest  position  is  indicated  by  point  A  on  both 
curves.  Now,  if  the  floor  instantaneously  dropped  a  distance  Az  the  isolator 
would  expand  a  distance  -  Ay  and  the  energy  stored  in  the  isolator  would  decrease 
from  point  A  to  point  B  on  the  E/mg  curve.  However  since  the  potential  energy 
of  the  mass  with  respect  to  one  floor  increased  an  amount  m^y  the  total  energy 
in  the  mass-isolator  system,  measured  with  respect  to  the  floor  would  be  equivalent 
to  that  at  point  C. 

Now,  since  we  have  postulated  a  conservative  system,  the  total  energy  of 
the  system  would  remain  constant  if  the  flooi  suffered  no  further  motion.  However, 
point  C  is  not  a  position  of  equilibrium  because  the  force  represented  by  point  M 
is  less  than  the  weight  of  the  isolated  mass.  The  mass  therefore  begins  to 
descend  losing  potential  energy  and  gaining  kinetic  energy.  Also  the  compression 
of  the  isolator  increases  so  the  isolator  gains  energy.  During  the  motion  the 
potential  energy  of  the  mass  would  be  represented  by  line  CH,  the  potential  plus 
kinetic  energy  of  the  mass  by  curve  CJH  and  the  potential  energy  of  the  isolator 
by  the  curve  BAD. 

At  the  point  of  maximum  displacement  the  mass  would  be  motionless  and  the 
kinetic  energy  zero.  Therefore  the  loss  in  potential  energy  of  the  mass  would 
be  just  equal  to  the  potential  energy  of  the  isolator.  Thus  point  D,  at  maximum 
displacement  can  be  established  by  projecting  line  CD  at  45°,  point  D  being  the 
intersection  of  the  line  with  the  isolator  energy  curve. 
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if  point  ^  and  point  D  are  connected  wioh  a  straight  line  the  ordinate 
between  line  ED  and  curve  BAD  is  the  kinetic  energy  of  the  mass  at  the  corres¬ 
ponding  displacement,  y. 

The  period  of  vibration  of  the  system,  T,  can  be  obtained  by  integrating 
the  reciprocal  of  velocity. 


t  -j^y 


dy 

E,_  1/2 


(6.7) 


For  a  real  system  the  floor  does  not  descend  instantaneously;  a  time  we 
have  denoted  t*  is  required.  During  time  t*  the  mass  descends  slightly  and 
gains  kinetic  energy  E^*. 

If  we  approximate  the  descent  of  the  floor  as  a  linear  function  of  time 
during  the  time  t*  and,  further  approximate  the  arc  MA  on  the  force  versus  y 
curve  by  a  straight  line  the  loss  of  potential  energy,  AE^,  is  given  by 


Ae  =  i  Ap  t*2g 
m  4 


_JL  „  1  (®)(t«2g) 
mg  I  rag' v 

and  the  gain  in  kinetic  energy  by 

k  o  m 


^  .  1  (4E)2  (t.2g) 

mg  H  vmg'  '  * 


(6.8) 


(6.9) 


(6.10) 


(6.11) 
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We  note  that  both  oi  these  corrections  will  be  quite  small  for  real  condi - 
tions  similar  to  those  illustrated  in  Figure  6.5. 

Figure  6.6  reproduces  the  energy  curve  of  Figure  6.5  and  illustrates  a 
method  for  correction  for  t*. 

Subsequent  to  time  t*  the  floor  of  the  structure  will  rise.  However,  since 
the  rise  is  due  to  the  decay  of  the  ground  pressure  wave  the  rate  of  rise  is 
slow,  a  second  to  several  seconds  depending  on  bomb  size.  This  rise  will  modify 
somewhat  the  response  of  the  system  but  the  slow  rise  alternately  adds  and 
subtrees  energy  from  the  oscillating  system  so  the  net  change  over  a  period  of 
time  is  small.  In  any  event,  for  a  forcing  function  of  small  duration  with 
respect  to  the  system  period,  neglect  of  the  rise  is  on  the  conservative  side, 
i.e.,  y  is  greater  for  no  rise. 

An  alternate  method  o.  plotting  the  energy  curves  may  be  preferable  for 
actual  design  purposes.  In  Figure  6.7  curve  (A)  is  the  normalized  isolator 
strain  energy  curve,  analogous  to  curve  BAD  of  Figure  6.5.  The  static  equilibrium 
position  of  the  mass  is  at  the  point  where  the  slope  of  curve  (A)  is  unity.  Curve 
(B)  is  the  normalized  mass  potential  energy  curve  and  curve  (C)  is  the  toual 
potential  energy  curve,  the  sum  of  curves  (A)  and  (B). 

Since  the  system  is  conservative  the  sum  of  potential  and  kinetic  energies 
must  remain  constant,  and  can  be  represented  by  the  horizontal  line  (D).  The 
difference  between  (c)  and  (d)  is  the  kinetic  energy  of  the  mass.  Points  (l) 
and  (2)  for  which  the  kinetic  energy  is  zero  are  points  of  extreme  relative 
displacement  between  mass  and  foundation. 
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6.6  RATTLE  SRACE  REQUIRED  UY  NONLINEAR  SYS  TIM;  RISE  TIME  CORRECTION 


Figure  6.7  ISOLATOR  ENERGY  CURVES 
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from 


The  rise  time  correction  can  be  made  by  computing  -  and  - 

mg  mg 

Equations  6.9  and  6.11.  Then  the  normalized  increment  of  total  energy  is 

AE  ^Em 
mg  mg  mg 

This  can  be  plotted  as  shown  on  Figure  6.7.  Then  the  range  of  mass  motion 
with  respect  to  the  floor  is  between  the  limits  of  points  1*  and  2‘  on  the  graph. 

Should  it  be  necessary  to  determine  the  response  of  a  nonlinear  system 
to  a  forcing  function  having  a  duration  of  a  quarter  period  of  more  one  has 
little  recourse  ot't  >r  than  to  carry  out  a  step-by-step  numerical  solution.  For 
a  complex  sy«+ -m  this  is  best  done  with  the  aid  of  a  digital  computor. 

6.2  Monconservative  Nonlinear  Systems 

If  the  general  conservative  nonlinear  system  shown  schematically  in  Figure 
6.1  is  modified  to  include  energy  dissipative  elements  it  can,  with  perfect 
generality,  be  represented  by  Figure  6.3.  In  Figure  6.8,  Ic  represents  an 
energy  storage  device  having  the  properties  of  a  conservative  isolator.  Block 
1^  represents  an  energy  dissipative  device.  No  energy  communicated  to  1^  is 
returned  to  the  system.,  n,v>\r  xhr  -t.-i-v  F.  always  resists  the  motion  y  regardless 
of  whether  y  and  its  derivatives  are  positive  or  negative. 

It  may  not  be  apparent  that  any  nonconservative  single  degree  of  freedom 
system,  regardless  of  its  complexity,  can  be  represented  by  a  scheme  as  simple 
as  Figure  6,8.  Energy  is  a  scalar  variable.  In  classical  mechanics  all  energy 
components  of  a  system  are  linearly  related.  Therefore  all  conservative  energy 
components  can  be  lumped,  if  we  so  choose,  and  all  dissipative  components  can 
also  be  lumped. 
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It  is  pointed  out,  however,  that  the  forces  Fc  and  F^  probably  would  not 
exist  separately  in  individual  physical  members. 

We  denote  the  energy  stored  in  isolator  I  by  E  .  The  curve  of  E  versus 

c  pc  pc 

y  can  be  deduced  from  the  physical  characteristics  of  any  real  system.  It 

simply  is  the  energy  within  the  isolator  at  position  y  with  all  time  derivatives 

of  y  zero.  The  F  versus  y  curve  is  then  simply  dE  /dy. 

c  pc 

No  energy  curve  analogous  to  Figure  6.2  exists  for  the  dissipative  element; 
energy  dissipated  depends  on  the  history  of  y  and  its  derivatives  as  well  as 
their  instantaneous  values.  However,  the  F^  curve  does  exist  but  is  much  more 
complicated  than  Figure  6.3.  F^,  in  many  real  system  is  a  function  of  at  least 
y  and  y.  It  may  depend  also  on  V  ,  the  absolute  acceleration  of  the  isolator. 

For  instance,  the  friction  against  an  unlubricated  shaft  turning  in  a  reamed 
hole  would  depend  to  a  degree  upon  the  acceleration  force  exerted  by  the  shaft 
against  the  hole.  Each  element  of  a  complex  system  would  have  a.  separate  time 
dependent  position  coordinate  so  that  if  we  attempted  to  lump  them  into  some 
generalized  coordinate,  w,  as  shown  on  Figure  6.8,  w  would  be  a  function  of  x,  z, 
and  all  of  their  derivatives. 

If  the  mass  of  the  isolator  can  be  considered  small  enough  to  neglect  in 
comparison  with  the  supported  mass,  m,  then  the  force  F^  is  a  function  of  y  and 
y  only.  The  dissipative  force  can  be  quite  accurately  approximated  by  a  function 
of  y  only  for  many  real  systems. 

Because  the  energy  dissipated  during  any  particular  time  interval  is  a 
function  of  the  previous  history  of  motion  the  technique  of  Section  6.1  cannot 
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be  extended  to  cover  the  dissipative  case.  This  is  because  the  simplicity  of 
the  technique  is  due  ultimately  to  the  fact  that  only  the  terminal  points  of 
motion  are  used.  If  the  technique  is  extended  to  account  for  the  effects  of 
intermediate  velocities  and  displacements  it  loses  its  simplicity. 

If  it  were  necessary  to  consider  damping  for  an  impulse  excited  nonlinear 
system  the  simplest  approach  would  be  to  use  the  phase  plane  method  which  is 
well  presented  in  the  literature*.  Actually,  for  a  system  of  a  fair  degree  of 
complexity  probably  the  best  procedure  would  be  to  program  the  problem  for  a 
digital  computer. 

Unless  the  nonlinear  system  is  quite  heavily  damped  acceptable  design 
estimates  of  rattle  space  and  acceleration  can  be  obtained  by  neglect  of  the 
damping  altogether  and  the  method  of  Section  6.1  used.  These  conditions 
generally  will  be  true  for  the  low  frequency  systems  used  in  shock  isolation  of 
sensitive  elements  of  hard  systems. 

Since  ground  shock  isolators  are  invariably  subjected  to  an  initial  high 
velocity  z  shortly  after  onset  of  the  ground  shock  wave  (within  several  milli¬ 
seconds)  velocity  sensitive  dampers  should  not  be  used;  they  would  communicate 
a  high  pulse  of  acceleration  to  the  isolated  equipment.  Coulomb  dampers,  however, 
are  insensitive  xo  velocity  and  further  have  the  desirable  property  that  the 
force  applied  to  the  isolated  equipment  can  be  made  nearly  constant. 

A  Coulomb  damped  nonlinear  system  proportioned  so  that  the  ratio  of  energy 
d4ssipatea  to  circulating  energy  for  the  first,  or  largest  ampxitude  cycle  of 

*See,  for  instance,  "Introduction  to  Nonlinear  mechanics" ,  IT.  Minorsky, 

J.  W.  Edwurds,  Ann  Arbor. 


vibration,  is  10  percent,  will  apply  a'  force  to  the  isolated  element  at  least 
10  percent  of  the  peak  conservative  element  force  Fc  but  less  than  20  percent 
of  the  peak  force  Fc<  Even  for  a  linear  conservative  element  (spring)  the 
peak  force  ratio  generally  is  less  than  15  percent. 

Such  a  system  would  ring  through  no  more  than  ten  periods  following  the 
decay  of  the  ground  shock  wave. 

6.3  General  Consideration  of  Nonlinear  Isolation  Systems; 

Their  Advantages  and  Disadvantages 

An  important  advantage  of  some  nonlinear  systems  was  pointed  out  in 
Section  6.1;  for  a  given  supported  weight  the  equilibrium  energy  of  a  "bi -linear” 
type  system  can  be  much  smaller  than  an  equivalent  linear  system  having  approxi¬ 
mately  the  same  frequency  of  perturbations  about  the  equilibrium  point.  Stated 
another  way,  a  truly  bi-linear  system  needs  far  less  energy  capacity  at  a  given 
frequency  and  excursion  range  than  does  a  linear  system  of  the  same  frequency. 

A  second  advantage  of  Coulomb  damped  systems  discussed  in  Section  6.2  is 
that  reasonable  damping  of  the  ringing  of  the  system  following  ground  shock 
excitation  can  be  obtained  without  the  consequent  high  accelerations  that  would 
be  communicated  to  the  isolated  equipment  by  velocity  sensitive  dampers.  Further, 
Coulomb  (constant  dry  friction)  dampers  (e.g.,  brake  drums)  are  probably  simpler 
and  more  reliable  than  any  other. 

The  characteristics  of  nuclear  blast  induced  ground  motion  are  such  that 


below  a  certain  vibration  frequency  the  excursion  of  the  isolated  mass  (rattle 
space)  remains  constant.  (This  is  represented  by  the  left,  displacement,  branch 


of  the  shock  spectrum. )  However,  as  the  frequency  of  the  system  is  reduced,  the 
acceleration  of  the  isolated  mass  is  reduced  and,  per’aps  of  more  importance, 
the  kinetic  energy  delivered  to  the  mass  is  reduced  proportionately.  Therefore 
its  tendency  to  oscillate  after  the  ground  shock  wave  has  passed  is  reduced. 

In  Section  6.k  some  practical  systems  that  can  he  built  to  support  sizeable 
masses  at  frequencies  two  or  even  three  orders  of  magnitude  lower  than  the  1  c"s 
range  of  most  present  hard  isolation  systems  are  briefly  discussed. 

The  disadvantages  of  nonlinear  systems,  with  one  notable  exception,  are 
three. 

1.  The  mechanical  complexity  of  nonlinear  systems  generally  is  considerably 
greater  than  that  of  linear  systems. 

° .  The  mathematical  problems  in  designs  of  nonlinear  systems  are  quite 
formidable  as  soon  as  the  single  degree  of  freedom  elements  are  departed  from. 

Even  for  the  elementary  single  degree  of  freedom  systems  the  computations  may 
require  a  digital  computor  for  any  other  than  impulse  loading. 

3.  The  coupling  between  elements  of  a  complex  nonlinear  system  can  be 
very  strong,  introducing  detrimental  motions  to  the  isolated  element  as  well  as 
grossly  complicating  the  analysis. 

The  exception  mentioned  is  the  Coulomb  damped  mass-spring  system.  Its 
mechanical  complexity  is  generally  less  than  that  of  any  other  damped  system 
and  its  repons e  to  an  arbitrary  input  can  be  determined  by  simple  graphical 
or  numerical  techniques  (Reference  8). 
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6.4  Negative  Spring  Isolators 

If  the  slope  of  the  energy  curve  of  Figure  6.2  "became  negative  then  the 
force  of  Figure  6.3  would  "become  negative  and  within  the  negative  region  the 
isolator  would  have  a  variable  negative  spring  constant. 

A  toggle  provides  a  simple  physical  example  of  such  a  device.  Figure  6.9 
schematically  shows  the  force -displacement  curve  for  a  toggle.  It  has  a  sub¬ 
stantial  negative  spring  region  centered  on  the  snap  through  position.  Further, 
by  proper  design  a  nearly  linear  region  of  about  1/3  h  (see  Figure  6.9)  can  be 
obtained. 

Negative  spring  devices  are  not  limited  to  toggles.  Any  energy  storage _ 

device  having  a  negative  energy  versus  displacement  curve  following  an  initially 
positive  curve  has  a  negative  spring  characteristic.  Mechanical  devices  can 
consist  of  cams  sliding  against  loading  rollers,  (Figure  6.10)  the  rollers  being 
loaded  by  any  appropriate  means.  Note  also  that  by  preloading  the  cam  followers 
the  lower  portion  of  the  cam  can  be  truncated. 

A  negative  spring  connected  in  parallel  with  a  positive  spring  having 
nearly  the  same  spring  constant  (slightly  higher  for  stability)  yields  a  spring 
system  of  substantial  energy  capacity  but  practically  zero  spring  constant  at 
the  operating  range  (Figure  6.11). 

Reference  5  describes  a  negative  spring  system  actually  installed  that 
isolated  a  ten  ton  suspended  mass  to  10  ^  g. 
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Figure  6.0  TOGGLE  NEGATIVE  SPRING  CHARACTERISTIC 


To  Isolated  Equipment 


Figure  6 . 10  schematic  of  cam  operated  negative  swing 
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The  principal  adverse  criticism  of  a  negative  spring  system  is  that  it 
is  bulky.  However  it  is  believed  that  with  sufficient  design  effort  the  bulk 
could  be  greatly  reduced.  If  this  conjecture  proves  to  be  true  they  would 
then  provide  particularly  desirable  elements  for  hard  installation  shock 
isolation.  The  only  input  parameters  required  are  stroke  and  suspended  weight. 
(Even  for  the  toggle  system  the  suspended  weight  can  be  varied  within  limits 
for  any  particular  isolator  and  still  maintain  the  frequency  at  near  zero). 

In  operation  the  isolated  mass  moves  only  minutely;  therefore  no  significant 
kinetic  energy  is  gained  by  the  mass  and  there  is  no  coupling  between  isolation 
elements  nor  ringing  subsequent  to  ground  shock  excitation.  The  isolated  ma~. 
remains  at  1  g  vertical  throughout  ^he  attack. 
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SECTION  7 


ISOLATOR  ELIMINATION 

The  possibility  that  isolators  can  be  eliminated  from  the  system  when 
ruggedized  equipment  is  used  should  be  checked. 

Any  criterion  for  isolator  elimination  must  be  approximate;  variations  in 
equipment  characteristics  are  too  great  to  allow  development  of  an  accurate 
general  method.  It  is  believed  that  the  method  presented  here  for  development 
of  a  criterion  for  any  particular  installation,  though  quite  conservative  still 
will  serve  to  indicate  that  isolators  can  be  omitted  in  many  appcl-Pic  instances-. . — 

Before  proceeding  with  development  of  the  criterion  an  aspect  of  the  shock 
isolation  problem  that  has  not  been  touched  in  other  parts  of  this  report  will 
be  considered. 

It  is  a  matter  of  experience  that  certain  peculiar  effects  take  place  in 
a  true  shock  environment  which  theory  is  hard  put  to  explain.  For  instance  it 
has  been  observed  that  if  a  large  charge  of  explosive  is  detonated  under  water 
and  close  to  a  ship,  though  the  ship  structure  aod  machinery  may  suffer  no  observable 
damage,  certain  brittle  elements  are  destroyed.  The  ship  may  be  left  undamaged 
but  in  darkness  because  the  filaments  of  incandescent  lamps  are  broken. 

In  an  attempt  to  "explain"  these  phenomena,  Weiner  (Reference  10)  postulated 
that  such  damage  is  caused  by  very  sharp  rises  of  stress,  rises  occurring  in  one 
or  two  microseconds.  Though,  to  our  knowledge,  Weiner's  predictions  and  conclu¬ 
sions  have  not  been  tested,  his  hypothesis,  if  accepted,  does  account  for  many 
shock  peculiarities. 
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Further,  Weiner  shows  that  conventional  shock  tests  would  not  disclose  such 
equipment  weaknesses  because  the  xise  times  of  the  tests  are  too  long.  In  addition, 
Weiner  points  out  that  many  materials  used  in  construction  have  concave  stress- 
strain  curves  and  that  sharp  pulses  would  be  generated. 

His  report  also  points  out  that  if  the  sharp  rise  postulated  is  the  reason 
for  failure  of  small,  more  delicate  elements  of  equipment  then  such  damage  can  be 
prevented  by  a  simple  technique:  mount  the  equipment  on  brackets  so  constructed 
that  no  stress  can  be  communicated  from  input  to  equipment  by  any  path  that  does 
not  include  a  flexural  element.  Weiner  showed  that  the  postulated  damaging  pulses, 
though  propagated  considerable  distances  as  tension  or  compression  waves,  are 
attenuated  to  the  vanishing  point  when  propagated  as  bending  through  a  distance 
of  one  or  two  beam  thicknesses. 

Now  the  required  "beams"  could  be  simply  wall  brackets  made  of  bent  plates 
and  having  elements  in  three  mutually  perpendicular  planes.  Most  bent  plate 
brackets  would  qualify. 

In  the  following  development  it  is  assumed  that  sharp  spikes  of  stress 
developed  inside  the  structure  and  mounting  equipment  are  avoided  by  use  of 
simple  mounts,  or  brackets,  so  constructed  that  there  are  no  direct  compression 
or  tension  paths  between  equipment  and  input;  somewhere  along  the  line  all  re¬ 
actions  are  resisted  by  bending  moments.  It  imposes  no  great  hardship  on  the 
designer  to  specify  that  the  natural  frequency  of  any  bending  element  in  the 
bracket  be  no  higher  than  500  cps. 


In  brief,  the  criterion  for  elimination  of  isolators  is  expressed  as  an 
equivalent  drop  to  a  hard  surface.  As  is  repeatedly  inadvertently  demonstrated 
in  electronics  laboratories,  even  sensitive  test  equipment  can  withstand  a  drop 
of  a  few  inches  to  a  hard  floor.  On  a  simple  bracket  or  flexible  mount  much 
equipment  can  stand  a  drop  of  two  or  more  feet. 

The  question  is  "How  high  a  drop  corresponds  to  the  interior  motion  of  a 
hard  shelter?" 

This  question  can  be  answered  with  fair  accuracy  by  the  five  following  steps 

1.  Estimate  the  free  field  motion  as  accurately  as  possible  in 
accordance  with  the  methods  of  Section  4. 

2.  Estimate  the  resulting  interior  structure  motion  ir.  accordance 
with  the  methods  of  Section  5-3*  Carry  out  at  least  a  rough 
solution  of  Equations  (5.29)  and  (5.32)  to  obtain  the  early 
history  of  interior  motion. 

3.  Using  the  Newmark -Hansen  Report  (Volume  in  of  this  report)  sketch 
the  peak  relative  response  spectrum  for  a  combination  of  simple 
pulses  that  approximate  the  early  history  of  interior  structure 
motion  obtained  in  Step  2. 

4.  Determine  the  peak  pseudo  velocity  from  the  spectrum.  Denote  it 

A 

V 

max 

5.  Compute  the  equivalent  drop,  x  ,  from 

i  ^  2 

1  max 
Xeq  “  2  g 

where  g  is  the  gravity  constant. 


No  exact  rules  for  the  maximum  drop  that  can  he  withstood  by  various  items 
of  equipment  can  be  given.  A  typical  value  of  maximum  interior  velocity  is  4  ft  ! 

per  second  having  an  equivalent  drop  height  of  only  3  inches.  Most  of  the  j 

equipment  in  a  hard  structure  could  resist  a  drop  of  3  inches  without  damage.  | 

1 

* 

It  may  be  desirable  to  have  available  a  rough,  conservative  estimate  of  the  j 

equivalent  drop.  This  can  be  obtained  for  vertical  motion  from  Figures  7.1  j 

through  7.3.  The  equivalent  drop  for  the  horizontal  component  can  be  taken  as  » 

i 

the  product  of  the  vertical  drop  and  the  tangent  of  the  Mach  angle  of  the  j 

compression  wave.  ! 

i 

( 

i 
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SECTION  8 

DETERMINATION  OF  INPUT  PARAMETERS 


Even  though  some  of  the  input  parameters  entering  into  the  ground  shock 
isolation  problem  can  only  be  estimated  it  is  necessary  in  order  to  carry  out  a 
design  to  have  firm  procedures  for  making  the  estimates.  In  this  section  formulas 
and  data  useful  in  making  these  estimates  are  compiled.  An  attempt  has  been 
made  to  simplify  some  of  the  formulas  at  the  expense  of  theoretical  accuracy. 

The  accuracy  remaining  however  is  more  than  is  generally  significant  in  view 
of  the  inherent  errors  in  input  data. 


8.1  Simplified  Structure  Frequency  Formulas 


Notation 


f 

f 

n 

n 

c 

h 

R 

l 

e 


frequency  of  vibration 

[cycles/sec] 

frequency  of  vibration  of  nth  mode 

^ycles/sec ] 

modal  index 

[dimensionless] 

sonic  velocity  of  shell  material 

[ft  per  sei ] 

shell  thickness 

frt] 

shell  radius 

[ft] 

distance  between  circumferential  nodal 

lines  [ft] 

half  central  angle  of  arch 

[radians] 

All  frequencies  of  vibration  are  computed  for  structures  vibrating  in  vacuum. 


8.1.1  Domes 

For  all  modes  c 

f  =  &rR 


(8.1) 


J  t  ' 


8.1.2  Circular  Cylinders 


For  all  extensional  (breathing)  modes 

f  = 

2rrR 


(8.1) 


For  deformational  (bending)  modes  having  no  circumferential  nodal  lines  and 
free  ends 

(8.2) 


»  „  _ 1_  n(n  -1)  /h,*]  _c_ 

n  Uv3  i>2+iii/2  R  J  ^ 


These  modes  have  2n  longitudinal  nodal  lines. 

For  deformational  modes  having  circumferential  nodal  lines  spaced  a  distance 
r  _  _  -s 


rni  *  f  n  <  1  + 


7TR 

-4- 

n(n  -l) 


[(2n2-l)(n2-l)  +  n2  (^ f  +  6  (f)2]  f-(8.3) 


fn  is  computed  by  Formula  (8.2). 


8.1,3  Circular  Arches 


For  deformational  modes  having  no  circumferential  nodal  lines  and  free  ends 


„  T  i  rn(rD  -i)  *  -I 

n '  L  ^ 


c 

27TR 


where 


P  _  n  7T 
rn  -  2  & 


(8.4) 


(8.5) 


These  modes  have  n  nodal  lines  between  the  spring  lines . 


For  deformational  modes  having  circumferential  nodal  lines  spaced  a 
distance  1 


f  ,  =  f  -< 
nj5  n 


1  + 


ttR 

"I 


r  f  -i) 

n  *  n 


2  r  1/2 

(£Tn2-l)(rn2-l)  +i'12(y)2+  6(|)2  y  (8.6) 


f  and  r  are  computed  by  formulas  (8.4)  and  (8.5)  respectively, 
n  n 


8.1.4  Rectangular  Slabs 

Fo-r  bending  modes  having  nodal  lines  spaced  at  "a"  in  one  direction  and 
spaced  at  "b "  in  the  perpendicular  direction 


or  f  =  ———  ch  (8.8) 

4  Vi  A" 

where  i  is  the  rectangle  diagonal  and  A  is  its  area. 

8.2  Normal  Mode  Geometry 

Unless  for  some  reason  a  detailed  analysis  of  an  interaction  problem  is 
to  be  carried  out  the  normal  mode  geometry  will  not  be  required.  If  such  an 
analysis  is  to  be  carried  out  the  modes  can  be  obtained  for  most  structures  from 
the  literature. 

If  the  method  of  Section  5-4  is  used  modal  geometry  is  not  required;  the 
statically  computed  deflection  is  used  instead. 


8.3  Soil  Parameters 

8.3*1  Bedrock  Elevation  and  Material 

These  data  should  be  obtained  from  borings.  If  not  available  they  can  be 
estimated  from  topographic  maps.  If  these  are  not  available  estimate  bedrock 
to  be  150  ft  below  surface  and  having  a  seismic  velocity  of  10,000  ft  per  second. 


8.3.2  Stress-Strain  Parameters 


If  at  all  possible  a  confined  compression  test  stress-strain  curve  should 
be  obtained  from  soil  samples. 


If  a  design  must  be  carried  out  to  be  applicable  to  a  number  of  unspecified 
sites  assume  an  "elastic"  soil  weighing  100  lbs/ft3  having  a  compression  wave 
velocity  of  1000  ft  p^r  second,  and  a  percyclic  damping  factor  of  0.25.  (The 
latter  number  is  quite  conservative;  a  more  realistic  value  would  be  O.35.) 

If  a  minimal  description  of  soil  conditions  at  a  site  is  unavailable,  elastic 
theory  and  the  following  approximate  values  can  be  used. 


Soil  Type 


Approximate  Density 


Top  sof  light 
dry  or  ujist 
loamy  silt 

Clayey  top  soil,  semi- 
consolidated  sandy  clay, 
nubble  or  gravel, 
loose  rock  talus,  wet  loam 

Cemented  sand, 
sand  and  clay 

Saturated  consolidated  clay 


100  lbs/ft3 
100  lbs/ft3 


100-110  lbs/ft3 
120  lbs/ ft3 


8.3.3  Percyclic  Damping  Factors  for  Real  Soils 


Approximate 
Compression  Wave 
Velocity 

650  ft/ sec 


1200  ft/sec 


2000  ft/ sec 
3000  ft/sec 


The  following  values  have  been  extracted  from  the  geophysical  literature. 
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Material 
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Fine  sand  and  silt 


Course  sand 


Shale 


Limestone 


Sandstone 


Granite 


Percyclic  Damping  Factor 


O.35  (O.25-0.50) 


0.0^3 

O.Ql+7 


0.060 


0.030 


1 
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APPENDIX  A 

DEVELOPMENT  OF  MOTION-TIME  METHOD  FOR  DISTRIBUTED  SYSTEMS 


A,1  Statement  of  the  Problem;  Method  of  Approach 

Assume  a  position  in  the  free-field  and  the  time  history  of  free-field 
pressure  at  the  salected  position.  The  time  history  is  presumed  to  include  all 
effects  that  influence  the  pressure  waves  arriving  at  the  selected  position, 
such  as,  the  incident  compression  wave,  the  incident  shear  wave,  the  directly 
transmitted  wave,  and  any  reflected  or  defracued  components  of  any  of  the 
foregoing  waves.  In  this  section  a  method  will  be  developed  for  determination 
of  the  motion-time  history  of  the  interior  of  a  structure  placed  at  the  selected 
position  in  the  free-field.  The  word  "method"  in  the  previous  sentence  is 
underlined  to  emphasize  that  a  general  procedure  is  developed,  not  a  series  of 
formulas  to  be  applied  in  specific  instances.  The  method  presented  admits  of 
gross  approximation  or  of  application  to  problems  deserving  of  a  high  degree  of 
accuracy  for  which  an  electronic  digital  computer  would  be  necessary  to  carry 
out  the  computations.  Thus,  if  only  the  free-field  compression  wave  can  be 
estimated  with  any  degree  of  accuracy  a  solution  can  be  obtained  within  this 
degree  of  accuracy.  On  the  other  hand,  if  site  conditions  are  sufficiently 
well  known  so  that  all  of  the  pressure  components  stated  above  are  known  within 
tolerable  accuracy  then  this  information  can  be  reflected  in  the  computed  interior 
structure  motion,  if  desired. 
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A  difficulty  with  all  engineering  methods,  and  a  difficulty  that  becomes 
more  pronounced  with  ir  'easing  generality  of  the  methods,  is  that  they  must  be 
understood  for  their  successful  application.  It  is  not  enough  merely  to  determine 
which  numbers  should  be  substituted  into  a  series  of  formulas  in  order  to  compute 
a  resulting  number  desired.  For  this  reason  in  the  following  subsection  the 
mathematical  expressions  U6ed  will  oe  developed  in  considerable  detail  and  the 
emphasis  will  be  placed  upon  the  physical  meaning  of  the  mathematical  expressions 
used.  This  is  considered  to  be  very  important  because  in  actual  application  most 
of  the  precise  mathematical  expressions  will  be  approximated  by  their  physical 
counterparts.  The  ease  and  facility,  therefore,  with  which  a  given  degree  of 
accuracy  can  be  achieved  in  the  results  will  hinge,  to  a  great  extent,  on  the 
thorough  understanding  of  the  processes  used  by  the  designing  engineer  and  good 
intuitive  estimates  of  the  accuracy  of  approximations  used. 

The  general  methods  of  mathematical  applied  mechanics  are  so  well  understood 
that  their  application  has  become  almost  standardized  in  technique  amounting,  in 
a  sense,  to  something  of  a  cult.  Without  meaning  to  disparage  the  considerable 
mass  of  .  ccellent  solutions  that  have  been  obtained  by  this  approach  it  is  pointed 
out  that  sometimes  the  formal  mathematics  tends  to  obscure  the  physics  of  the 
problem  being  investigated.  The  development  of  the  next  subsection,  in  particular, 
does  not  approach  the  problem  in  the  classical  manner.  Considerable  physical 
insight  is  gained  by  this  departure  from  custom. 

A. 2  The  Differential  Equation 

Consider  a  linear,  elastic,  statically  loaded  shell  structure*.  Its 

*  Since  shells  are  the  structures  of  interest  this  development  is  carried  out  in 
terms  of  shell  characteristics.  Appropriate  modifications  render  it  applicable 
to  any  linear,  elastic  structure. 


A  -£ 


equation  can  be  represented 

D  L  (u)  =  p  (A.l) 

where 

D  is  a  scalar  constant 

u  is  the  absolute  vector  displacement  of  a  point  on  the  shell 

p  is  the  vector  pressure  at  a  point  on  the  shell 

L(u)  is  a  linear  vector  operation,  a  function  of  u  and  its 

spacial  coordinates.  Equation  (A.l)  can  be  taken  as  the 

definition  of  L(u) . 

It  should  be  appreciated  that  the  pressure  vector  p  of  equation  (A.l) 
includes  the  forces  that  we  normally  term  reactivi  .  in  other  words,  if  the 
pressure  p  was  summed  over  the  entire  structure  the  resultant  would  be  zero, 
and  the  moments  of  the  pressure  p  summed  over  the  entire  structure  must  be  zero. 
The  following  elementary  examples  are  included  to  fix  the  meaning  of  the  linear 
vector  operator  and  the  circumstances  under  which  the  vectorial  nature  of  u  and 
p  become  significant. 

For  a  flat  plate  laterally  loaded  the  analog  of  Equation  (A.l)  is 
D  V2  ^(u)  =  p 

where 

V2  = 


Laplacian  operator 


E  ■  modulus  of  elasticity 

h  *  plate  thickness 

t>  *  Poisson's  ratio 

It  should  be  observed,  however,  that  since  the  displacement  and  pressure 

vectors  are  in  the  same  direction  over  the  entire  surface  of  the  plate,  the 
vector  character  of  displacement  and  pressure  do  not  contribute  anything 
significant  to  this  particular  problem  and  therefore  they  might  as  well  be 
written  as  scalars  (as  they  usually  appear  in  the  literature).  If  the  flat 
plate  were  of  nonuniform  thickness  then  the  equation  would  be 

1  vWro  -  5 

(l-» ) 

For  an  inextensible  circular  arch  of  central  angle  20  and  radius  R  the 
analogue  of  Equation  (A.l)  is 


where 

El  is  the  stiffness  of  the  arch 
p  is  the  vector  pressure  per  unit  length  of  rib. 

Now^for  the  arch 

u  =  uy  n  +  up  l 

where  n  and  I  are  unit  vectors  in  the  radial  and  tangential  directions 
respectively . 


Also,  for  the  inextensible  theory 


u 


r 


d 


Since  each  component  of  the  vector  displacement  must  satisfy  the  differential 
equation 


p 


r 


where 

P  *  Pr  n  +  ^  1 

we  see  that,  although  the  vector  L(u)  must  be  colinear  with  the  vector  p,  in 
general  the  vector  u  is  not  colinear  with  p.  Physically  the  situation  will 
arise  for  any  structure  that,  under  the  influence  of  a  concentrated  load,  moves 
laterally  with  respect  to  the  load . 


Equation  (A.l)  is  not  sufficient  in  itself  to  determine  an  actual  deflection 
even  though  the  loads  are  given.  It  only  expresses  the  relation  between  the 
deflection  at  one  point  to  that  at  another.  In  order  4_o  get  absolute  values  for 
the  displacements  we  must  prescribe  a  set  of  boundary  conditions  to  the  problem, 
i.e.,  a  set  of  points  where  the  displacements,  moments,  shears,  or  slopes  are 
zero.  The  complete  set  of  equations  for  a  particular  problem  then  consists  of 
Equation  (A.l)  plus  the  appropriate  boundary  conditions. 


It  is  well  to  emphasize  at  this  point  that  for  any  particular  problem  an 
equation  analogous  to  (A.l)  always  exists  buu  iu  may  be  Loo  compjex  mathematic 'illy 


to  use  or  even  to  write.  An  example  of  an  essentially  simple  physical  system 
that  does  not  have  a  tractable  equation  analogous  to  (A.l)  is  a  steel  building 
frame;  however,  the  general  conclusions  that  can  be  drawn  by  operating  with 
Equation  (A.l)  are  applicable  to  the  steel  building  frame. 

The  equations  of  free  vibration  of  the  above  system  in  vacuum  are  then 
(damping  neglected) 

•  • 

D  L  (u)  +  mu  =  0  (A. 2) 

Bn  (u)  -  0  (n-1,  2,...)  (A. 3) 

where 

m  is  the  unit  mass  of  the  shell,  a.  function  of  the  spacial  coordinates 
Bn(  )  are  the  boundary  conditions 

a  is  an  index  denoting  the  individual  boundary  conditions. 

Equations  (A. 2)  and  (A. 3)  can  always  be  solved  in  the  form 
1=00 

u  -  Z  A.  </>..  sin  (co.t  -  9.)  (A. 4) 

i_l  i  ■  ±  x  l  '  ' 

where 

A^  are  arbitrary  constants 

and  9^  are  arbitrary  phase  angles 

<b ^  are  a  set  of  orthogonal  vector  functions  of  the  spacial  variables  only, 
satisfying  the  boundary  conditions  and  the  orthogonality  relation 

/  <T.  •  0.  m  dA  =  0  ...  if  i  f  j  (A. 5) 

A  1  J 

=  N.  ...  if  i  =  j. 

i 
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The  area  integral  is  taken  over  the  entire  surface  of  the  shell. 


The  Nj  are  constants,  termed  the  norms,  of  the  orthogonal  system.  Then 
the  circular  frequencies,  satisfy  the  relations 

D  L  (0i)  =  moy2i  (A. 6) 

The  physical  meaning  of  Equations  (A. 4)  th  'lugh  (A. 6)  is  that  any  steady 
state  vibration  u  can  be  represented  as  the  Siam  of  an  infinite  number  of  rather 
special  displacement  configurations,  0^  vibrating  sinusoidally  at  specific 
frequencies  each  frequency  being  related  to  the  associated  displacement 
configuration  by  Equation  (A. 6).  Particular  displacement  configurations,  u, 
are  obtained  by  adjusting  the  values  of  the  arbitrary  constants  and  @ 

Physically  it  is  by  no  means  obvious  that  this  situation  should  be  true. 

2 

If  we  regard  co  1  in  Equation  (A. 6)  as  simply  an  arbitrary  number  then  certainly 

the  equation  could  be  solved  for  some  function  We  might  make  a  guess  as  to 

—  2 
the  general  shape  of  0 ^ substitute  it  into  the  linear  operator,  divide  by  m  cb 

and  obtain  a  better  approximation  for0i-  If  this  procedure  were  carried  out  an 

infinite  number  of  times  we  would  expect  to  find  seme  particular  displacement 

configuration  satisfying  Equation  (A.6)  for  the  particular  value  ofo^  used. 

• 

In  this  manner  could  be  generated  for  any  arbitrary  value  of  od^;  however,  it 
would  be  found  that  only  those  functions  0  corresponding  to  certain  particular 
values  of  cb^,  would  also  satisfy  th'e  boundary  conditions.  We  would  still  have 
an  infinite  number  of  0^  left  satisfying  the  boundary  conditions  but  this  infinity 
would  be  smaller  than  the  number  that  did  not.  A  fairly  long  mathematical 
development  then  shows  that  any  two  functions  0 .  and  0,  satisfying  both  Equation 
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(A. 6)  and  the  boundary  conditions  also  satisfy  Equation  (A. 5)*. 

Evidently  since  the  geometry  of  the  <j^  are  influenced  by  the  boundary 
conditions  this  influence  is  reflected  in  the  associated  frequencies  of  vibration 
03^ .  Thus  if  two  structures  identical,  except  as  to  boundary  conditions,  were  to 
be  placed  in  steady  state  vibration  the  vibrations  of  the  two  structures  might 
be  markedly  different.  The  most  familiar  example  of  this  situation  is  the 
difference  in  vibrations  between  a  hinged  end  and  a  fixed  end  beam  that  are 
otherwise  identical. 

Since  the  element  of  mass,  m  dA,  is  always  positive  it  can  be  concluded 
from  Equation  (A.  5)  that  in  general  0^  and  0^  will  have  positive  and  negative 
portions,  otherwise  their  product  would  not  average  out  to  zero  over  the  surface 
of  the  structure.  In  most  treatises  on  orthogonal  functions  the  functions  are 
considered  to  be  scalars  rather  than  vectors  and  to  be  functions  of  one  independent 
variable  only  rather  than  two  required  for  a  general  description  of  a  shell 
structure,  the  two  dimensional  functions  for  the  shell  being  represented  as  the 
product  of  two  one  dimensional  functions  each  of  one  of  the  coordinates  of  the 
shell  surface.  Such  one  dimensional  functions  and  two  dimensional  functions 
comprising  a  product  of  two  one  dimensional  fractions  have  a  series  of  nodal 
points  or  nodal  lines,  i.e.,  points  or  lines  at  which  the  value  of  the  function 
is  zero.  For  the  two  dimensional  case  the  nodal  lines  divide  the  function  into 
alternate  positive  and  negative  zones  somewhat  similar  to  a  checkerboard.  As 
the  index  number,  i,  increases,  the  number  of  nodal  lines  on  such  modes  increases 
also. 

*  See,  for  instance,  Air  Force  Special  Weapons  Center  Technical  Report  TR-Sy-2 
"Protective  Construction1',  T.G.  Morrison,  Part  III  (u). 


The  vector  functions  ¥  in  general  do  not  exhibit  this  characteristic. 

The  scalar  value  of  the  vector  displacement  is  always  positive  or  possibly  zero 
and,  in  general,  if  the  vector  is  resolved  into  its  comnonents,  at  a  point  where 
one  component  is  zero  the  other  components  will  not  be  zero;  therefore,  the 
vector  modes  do  not  have  clearly  defined  nodal  lines .  The  components  of  the 
modes,  however,  do  have  nodal  lines  but  the  nodal  lines  of,  say,  the  radial 
component  and  the  lateral  component  of  the  vector  modes  of  a  circular  cylinder 
occur  at  different  positions.  This  observation  is  quite  important  when  setting 
up  the  technique  for  numerical  evaluation  in  a  later  section. 

If  the  structure  is  subjected  to  a  dynamic  load  p  then  the  equations  of 
motion  become 

D  L  (u)  +  m  u  =  p 

Bn  (u)  *  0.  (A. 7) 

Equations  (A. 7)  have  a  solution  in  the  form 

00 

u  -  Z  T.  (t)  ¥,  (A- 8) 

.  ,  i  i 

i=l 

where  T^(t)  are  the  solutions  of  equations  of  the  form 

T.  +  05. 2  T.  =  r.it)  (A. 9) 

i  li  i 

The  functions  t.  arc  given  by 

i 

T At)  -J-J  ?.J  dA.  (A. 10) 

i  A 
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In  subsequent  sections  considerable  space  will  be  devoted  to  methods  for 
computing  the  value  of  -r^t).  It  is  worthwhile,  therefore,  to  get  some  feeling 
for  the  physical  significance  of  this  expression.  Since  the  dimensions  of  p 
are  pounds  per  square  inch,  the  dimensions  of  ^  are  inches,  and  those  of  dA  are 
square  inches,  the  integral  In  Equation  (A. 10)  simply  represents  the  total  energy 
that  would  be  delivered  to  the  mode  by  the  pressure  P  if  the  pressure  were 
held  constant  and  the  mode  allowed  to  move  from  a  position  of  zero  displacement 
to  its  maximum.  Thus,  the  product  N^T^t)  represents  a  time  varying  function  at 
any  instant  equal  in  value  to  the  energy  that  could,  at  that  instant,  be  delivered 
to  the  mode  by  the  pressure  function  p. 

We  will  now  set  up  the  equation  of  motion  of  the  same  structure  immersed  in 
an  elastic  medium  that  offers  resistance  to  motion  of  the  imbedded  structure. 

We  will  denote  the  absolute  motion  of  a  particle  on  the  surface  of  the  structure 
by  u,  the  absolute  motion  of  a  free-field  particle  by  z,  and  the  relative  motion 
of  the  structure  with  respect  to  the  free-field  by  w.  Thus, 

u  =  z  +  w. 

The  medium,  in  general,  offers  three  components  of  resistance  to  motion  of  the 
structure.  These  will  oe  discussed  in  some  detail. 

The  first  component  we  will  term  the  displacement  sensitive  component. 
Consider  a  cylindrical  shell  immersed  in  an  elastic  medium.  Now,  if  the  shell 
is  uniformly  expanded  outward  against  the  elastic  medium  a  pressure  will  be 
developed  by  the  medium  against  the  shell  at  the  interface,  this  pressure  being 
proportional  to  the  amount  of  outward  expansion.  In  a  truly  elastic  medium  if 
the  shell  were  bonded  to  the  medium  a  tension  would  be  developed  at  the  interface 
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if  the  shell  were  uniformly  compressed  inward.  The  proportionality  constant  j 

for  both  the  compression  and  tension  cases  would  be  the  same.  Now  suppose  that 

the  shell  was  nonuniformly  expanded  outwards  against  the  elastic  medium.  Then 

the  interface  pressure  developed  would  vary  from  point-to-point  and  would  be  a 

function  not  only  of  the  elastic  constants  of  the  surrounding  material  but  also 

of  the  geometry  of  the  deformation.  In  general,  we  could  write  j 

k  F(v)  v  -  pd  (/ • 11) 

where 

Pd  *  interface  pressure  due  to  relative  displacement  of 

structure  with  respect  to  free-field  3 

1 

_  j 

F(w)  -  a  dimensionless  function  of  w 

k  *  a  proportionality  constant. 

A 

For  the  case  of  uniform  expansion  or  contraction,  previously  discussed,  the 
dimensionless  function  F(w)  would  be  identically  unity. 

The  relative  motion  of  the  structure  with  respect  to  the  free-field,  w, 
can,  in  general,  be  resolved  into  three  components,  one  of  these  normal  to  the 
surface  of  the  structure  and  two  more  mutually  orthogonal  with  each  other  and 
the  normal  component  but  tangential  to  the  surface  of  the  structure.  The  latter 
two  components  could  be  combined  into  a  single  tangential  component,  the  direction 
of  which,  in  general,  would  vary  from  point-to-point  around  the  surface  of  the 
structure.  We  will  denote  the  normal  component  w  •  n  and  the  tangential  component, 
w  •  i  where  n  is  the  unit  vector  normal  to  the  she1  /  and  J  is  the  unit  vector 
tangent  to  the  shell  and  lying  in  the  plan  defined  by  n  and  w.  Now  in  a  truly 
elastic  medium  two  components  of  pressure  are  developed  at  the  interface  due  to 
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the  time  derivative  of  these  components  of  motion.  The  normal  component  of 
pressure  which  we  will  denote  pvn  is  given  by 

Pvn  “  ^  ^  ‘  ^  (A. 12) 

and  the  tangential  component  is  given  by 

a 

Pve  "  c2p  J  w  •  J  (A. 13) 

where 

c^  is  the  compressional  wave  velocity 
Cg  is  the  shear  wave  velocity 
p  is  the  unit  mass  of  the  medium 

Equations  (A. 12)  and  (A. 13)  are  independent  of  ye  stress  strain  law  and 

are  accurate  for  any  material  if  the  instantaneous  wave  velocities  are  used. - 

Practically,  their  principal  value  is  for  materials  having  constant  wave 
velocities  (elastic,  viscoelastic,  etc.).  It  is  pointed  out  that  these  equations 
are  written  only  for  the  interface  pressures  developed  and  do  not  represent  the 
time  history  of  the  waves  radiated  from  the  structure  which  ultimately  result 
from  the  interaction. 

If  the  structure  were  not  present  in  the  free-field  but  the  position 
occupied  by  the  structure  was  engulfed  by  a  pressure  wave,  the  mathematical 
surface  defining  the  structure-medium  interface  would  be  subjected  to  a  component 
of  pressure,  the  free-field  pressure.  In  a  solid  medium,  elastic  or  otherwise, 
the  free-field  pressure  can  be  represented  by  a  second  order  tensor.  We  will 
denote  this  free-field  pressure  tensor  p.  On  any  surface  arbitrarily  described 
within  this  second  order  tensor  field  the  vector  pressure  on  the  surface,  which 
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we  will  denote  by  p,  is  given  by 


p  *■  n  .  f?  (A.l4) 

where  n  is  the  vector  unit  normal  to  the  surface. 

Now  if  an  actual  structure  is  immersed  in  the  free-field,  the  vector  pressure 
at  the  st ructure -medium  interface  is  equal  to  the  sum  of  the  four  components 
enumerated  above.  The  pressure  component  that  we  ordinarily  regard  as  the 
reflected  wave  or  the  doubling  of  free-field  pressure  against  a  motionless  surface 
is  simply  the  sum  of  components  and  p^  due  to  the  relative  velocity  of  structure 
and  free-field  particles.  For  equilibrium  this  interface  pressure  must  be 
balanced  by  two  pressure  components  developed  by  the  structure: 

1.  A  component  due  to  the  deformation  of  the  structure,  which 

we  have  symbolized  by  DL(u),  and 

•  • 

2.  The  inertia  of  the  structure,  mu. 

Thus  the  equation  of  motion  of  a  structure  immersed  in  elastic  material  is 

•  • 

DL(u)  +  m  u  =  n  •  p  -  k  F(w)w  -  c np  n  w  •  n  -  c^p  1  w  .  ~l  (A. 15) 

The  minus  signs  appear  on  the  right  side  of  Equation  (A. 15)  because  for  an 
inwardly  directed  pressure  to  be  developed  at  the  interface  due  to  structure 
deformation  the  deformation  u  would  have  to  be  positive  inward.  Since  u  is  the 
sum  of  the  free-field  motion  and  relative  motion  w,  a  positive  w  would  repre¬ 
sent  a  -tendency  for  the  shell  surface  to  separate  from  the  free-field  particlep 
yielding  a  tension  stress  or  negative  pressure. 

Substituting  w  =  u  -  z  into  Equation  (A. 15)  and  rearranging  we  obtain 
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(A.16) 


DL(u)  +  k  F(w)u  +  c^P  n  u  •  n  +  c^p  X  u  •  X  +  mu 

“  n  •  p  +  k  F(w)  z  +  c^P  n  z  •  n  +  c^>  J  z  •  T 

With  the  exception  of  the  term  k  F(w)z  the  right  side  of  Equation  (A.l6)  is  a 
function  only  of  the  free-field  pressure,  free-field  particle  motion,  and 
structure  geometry.  Therefore,  with  a  suitable  approximation  for  F(w)  the 
right  side  of  Equation  (A.l6)  is  computable  from  initial  data.  We,  therefore, 
define  p^ 

P1  s  n  •  p  +  k  F(w)z  +  c  p  n  z  •  n  +  c^o  i  z  •  J  (A. 17) 

It  is  instructive  to  observe  the  form  taken  by  Equation  (A.l6)  for  certain 
limiting  cases.  Consider,  first,  that  the  soil  motion  and  structure  response 
are  so  slow  that  the  process  becomes  quasi-static,  i.e.,  all  time  derivatives 
become  zero.  Then  we  see  that  the  interface  pressure  would  be  just  the  sum  of 
free-field  pressure  and  a  component  due  to  the  distortion  of  the  soil  around 
the  structure,  induced  by  the  disparity  between  free-field  motion  and  structure 
motion. 

If  the  struc'.ure  were  a  rigid,  immovable,  flat  surface  parallel  to  the 
advancing  wave  front  then  DL(u)  would  vanish,  z  would  be  equal  to  u,  and  \ 
would  be  equal  to  zero.  It  is  not  apparent, but  for  this  case  the  function 
F(w)  would  be  zero  and  also  the  dot  product  z  •  ~l  would  be  zero.  Then  since  ~z . 

the  particle  velocity,  is  just  equal  to  for  a  plane  wave,  the  right  side 

Pc 

of  Equation  (A.l6)  would  become  simply  2p  as  it  should  for  this  condition,  the 
2  representing  one  hundred  percent  reflection. 


Equation  (A.l6)  is  extremely  difficult  to  solve.  The  term  kF(w)u  renders 
the  equation  nonlinear.  Also,  the  two  first  derivative  terms  having  different 
constant  factors  c^  and  c  mitigate  against  the  use  of  a  normal  mode  solution 
even  without  the  complication  of  the  nonlinear  term.  At  this  point,  then,  there 
would  be  two  general  courses  open  to  us.  One,  we  could  insist  on  a  mathematically 
rigorous  solution  and  carry  out  the  evaluation  of  Equation  (A.l6)  by  means  of  a 
large  digital  computer;  two,  we  could  simplify  Equation  (A.l6)  and  perhaps  carry 
out  the  evaluation  by  simpler  means.  When  we  consider  the  degree  of  uncertainty 
that  must  exist  for  the  several  parameters  entering  into  Equation  (A.l6)  the 
latter  course  seems  to  be  the  most  reasonable  to  adopt.  A  rigorous  solution  of 
a  problem  for  which  the  input  data  is  probably  accurate  to  no  better  than  a 
factor  of  ten  or  twenty  percent  yields  only  illusory  accuracy  in  the  final 
results.  However,  we  must  be  careful  in  simplifying  the  equation  that  the 
elimination  of  a  complicating  factor  does  not  change  the  qualitative  behavior 
of  the  system  and  that  the  simplifications  do  not  introduce  numerical  inaccuracies 
greater  than  can  be  tolerated. 

It  will  be  found  that  the  nonlinear  term  has  small  influence  on  the  net 
response  of  the  structure.  Actually,  it  contributes  only  a  small  increment  to 
the  "spring  constant"  of  the  structure,  by  far  the  larger  portion  being  contributed 
by  the  factor  DL(u)  for  most  cases.  The  sin^l^  exception  occurs  for  rigid  body 
motion  of  the  structure  m  which  case  the  linear  operator  contributes  zero.  But 
for  this  case  the  function  F(w)  is  identically  unity.  We  will,  therefore, 
approximate  the  nonlinear  term  by  the  linear  term  Ku  and  develop  rational  methods 
for  estimation  of  an  average  value  of  K  to  use  for  any  particular  structure. 
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The  velocity  terms  are  a  little  more  difficult  to  handle.  We  begin  by 
determining  graphically  the  meaning  of  the  sum  of  the  two  velocity  terms.  We 
note,  first  of  all,  that  c^  is  always  greater  than  c^,  generally,  on  the  order 
of  two  or  three  times  as  great  as  Figure  A. la  shows  the  graphical  construc¬ 

tion  of  the  vector  representing  the  sum  of  the  two  velocity  terms.  It  is  seen 
that  in  general  this  sum  will  have  a  numerical  magnitude  different  than  that 
of  z  and  that  the  direction  of  the  resulting  vector  will  differ  from  z.  In 
Figure  A. lb  the  two  wave  velocities  c^  and  have  been  replaced  by  some  average 
wave  velocity  c  such  that  the  length  of  the  resulting  vector  has  the  same  length 
as  that  appearing  in  Figure  A. la.  However,  this  approximation  yields  a  resultant 
having  the  same  direction  as  the  vector  z.  Now,  in  general,  for  a  normal  mode, 
the  angle  p  may  vary  between  zero  and  90°  for  different  locations  on  the  mode; 
however,  at  the  points  where  the  angle  |3  is  90°  the  length  of  the  vector  7  iS  a 
minimum  and  at  the  point  where  /3  is  zero  the  length  of  z  is  a  maximum.  Further, 
a  little  trigonometry  discloses  that  the  maximum  value  possible  for  the  angle  7 

is  about  18-1/2°  for  which  its  cosine  is  about  .95*  This  occurs  for  a  ratio  of 
^1  *^1 

—  =  4.  For  the  more  likely  value  of  —  *■  2,  the  maximum  value  of  7  is  less  than 
c2  c2  ^ 

16*.  Therefore,  since  the  net  effect  of  the  z  terms  is  to  remove  e^erg,/  irom. 
the  vibrating  system  we  can  make  a  very  accurate  approximation  by  substituting 
c p  z  and  cp  u  for  the  first  derivative  terras  in  Equations  (A.l6)  and  (A.17).  In 
a  subsequent  section  the  value  of  the  adjusted  c  is  computed  so  that  the  energy 
radiated  by  the  approximate  term  is  equal  to  the  energy  radiated  by  the  two 
accurate  terms .  It  is  found  that  for  most  conditions  c  will  be  equal  to  more 
than  9<3  percent  of  c^ 
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With  these  approximations  our  working  equations  become 


DL(u)  +  Ku  +  c  p  u  +  mu  *  p^  (A.l8) 

— -  m  • 

=  n  •  p  +  Kz  +  c  pz  (A.  19) 

A. 2.1  Physical  Meaning  of  Terms  of  Fundamental  Equations 

The  terms  in  Equations  (A.l8)  and  (A.19)  have  been  grouped  into  components 

representing  the  motion  of  the  structure  and  components  representing  the  motion 

of  the  soil.  The  pressure  p^  is  a  fictitious  quantity,  it  does  not  represent 

the  pressure  at  the  interface  between  the  structure  and  soil.  The  latter  would 

be  just  the  sum  of  the  pressure  components  necessary  to  maintain  the  structure 

•  • 

in  its  deformed  condition,  DL(u),  and  the  inertia  component,  mu.  All  of  the 
other  terms  represent  the  effects  of  the  free-field  pressure  and  soil  structure 
interaction,  however,  since  P^  is  immediately  computable  once  the  pressure  wave 
p  has  been  established,  it  forms  a  convenient  forcing  function  for  the  equation 
of  motion  (a. 18). 

Actually,  the  term  P,is  not  quite  as  fictitious  as  it  appears.  If  Equation 
(A. 18)  is  multiplied  through  by  u  ,  the  velocity  of  the  structure,  and  then 
integrated  with  respect  to  surface  area  and  time,  each  of  the  resulting  terms 

has  physical  significance. 

t  _  ^  t  _ 

/  /  DL(u)  .  u  dA  dt  +  K  /  /  u  ’  u  dA  dt 

0  A  0  A  (A. 20) 

t  _  t  ..  t  _ 

+  c  pj  /  u*udAdt  +  m//  u  •  u  dA  dt  =  /  /  p.  •  u  dA  dt 
o  A  0  A  0  A 
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Thus,  in  Equation  (A. 20)  the  first  integral  represents  the  potential  energy 
stored  as  strair  energy  in  the  entire  structure  between  the  times  z^o  and  t. 

(All  energies  are  here  computed  with  respect  to  the  fixed  coordinate  system.) 

The  second  integral  represents  the  recoverable  potential  energy  stored  as  strain 
energy  in  the  soil  surrounding  the  structure.  The  third  integral  represents 
the  energy  dissipated  during  time  zero  to  t  in  the  form  of  a  wave  radiated  from 
the  structure.  The  fourth  integral  represents  the  kinetic  energy  delivered 
to  the  structure  within  the  interval  zero  to  t.  Therefore,  the  integral  on  the 
right  hand  side  of  the  equal  sign  represents  the  total  energy  delivered  to  the 
soil-structure  system  during  the  time  intervals  zero  to  t.  Thus,  the  pressure 
represents  the  component  of  pressure  around  the  structure  available  to  deliver 
energy  to  the  soil- structure  system  during  motion  of  the  structure. 

Immediately  we  see  that  if  u  and  u  have  been  estimated.  Equation  (A.20) 
provides  a  useful  check  of  the  results.  This  is  particularly  true  for  rigid 
body  motion  because  in  this  instance  u  and  u  are  invarient  over  the  surface  of 
the  structure  and  the  first  integral  becomes  identically  zero.  Then  any  over¬ 
estimate  of  u  would  indicate  a  larger  total  energy  on  the  left  side  of  Equation 
(A.20)  than  right,  clearly  a  physical  impossibility. 

A. 2. 2  Range  of  Applicability  of  the  Basic  Equations 

Define  p2  ■  p-^  -  Kz  (A.2l) 

Substituting  (A.2l)  and 

u  =  z  +  w 
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into  Equation  (A. 20)  there  results 

t  __L  t  _  ^  t  _  . 

/  /  DL(u)  •  u  dA  dt  +  /  /  K  w  •  u  dA  dt  +  c p  /  /  u  *  u  dA  dt 
o  A  o  A  o  A 

t  W  t  _ 

+  m  /  /  u  •  u  dA  dt  -  /  /  P  .  u  dA  dt  (A. 22) 

o  A  o  A  d 

MU 

We  see  now  that  p^  is  not  influenced  by  K  and  that  the  second  integral 
appears  in  terms  of  w  •  u,  where  w  is  the  departure  of  the  structure  motion 
from  the  free -field  motion,  z. 

Equation  (A. 22)  provides  a  useful  means  for  estimating  the  degree  of 
approximation  introduced  into  the  solution  when  nonlinear  or  nonelastic  soils 
are  approximated  by  an  elastic  soil  with  suitable  adjusted  constants.  We 
observe  first  that  the  factor  DL(u)  is  numerically  equal  to  the  pressure  that 
must  be  distributed  over  the  surface  of  the  structure  in  order  to  maintain  it 
in  its  deformed  configuration.  Also,  we  note  that  even  if  the  structure  were 
highly  nonlinear  there  would  be  some  other  nonlinear  opera  Lor,  N(u),  for  instance, 
that  would  replace  the  linear  operator  throughout  the  entire  argument  to  this 
point.  Also,  this  nonlinear  operator  would  be  equivalent  to  the  static  pressure 
necessary  to  maintain  the  structure  in  its  deformed  condition.  Now,  whether 
elastic  or  inelastic,  the  pressure  sc  computed  is  equal  to  the  design  pressure 
for  the  structure,  or  more  properly  the  equivalent  static  design  pressure.  For 
most  underground  structures  this  equivalent  static  pressure  is  on  the  order  of 
the  peak  overpressure  that  the  structure  is  designed  to  resist. 
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The  factor  Kw  also  represents  a  pressure.  Thus,  we  can  compare  the  relative 
influence  of  the  first  and  second  integrals  of  Equation  (A. 22)  by  comparing  Kw 
with  peak  design  overpressure,  inasmuch  as  all  f  the  other  factors  in  the 
integrals  are  identical.  If  the  real  soil  were  highly  nonlinear,  and  a  linear 
approximation  substituted,  probably  the  best  approximation  would  be  to  take  a 
value  of  K  that  would  yield  the  same  energy  content  at  a  specific  displacement 
level  as  did  the  real  soil.  In  other  words,  the  triangular  area  under  the  stress- 
strain  curve  of  the  approximating  elastic  soil  would  be  made  equal  to  the  area 
under  this  stress -strain  curve  of  the  nonlinear  soil  at  its  expected  maximum 
strain.  It  is  interesting  to  note  that  although  the  second  integral  in  Equation 
(A. 22)  represents  reversible  energy  for  the  elastic  case,  for  the  inelastic  case 
it  represents  both  reversible  and  irreversible  (hysteresis)  energy  stored  in  the 
surrounding  soil.  In  the  elastic  case  all  of  this  energy  may  be  delivered  to 
the  structure  at  some  point  during  its  vibration  but  for  the  inelastic  case  only 
a  portion  could  be  delivered.  Thus,  the  net  effect  of  substituting  an  elastic 
equivalent  for  a  real  inelastic  soil  i  ;  slightly  conservative  because  the  possible 
maximum  relative  displacement  and  maximum  absolute  acceleration  of  the  structure 
are  proportional  to  the  square  root  of  the  energy  available  to  the  structure. 


To  gain  an  estimate  of  the  relative  magnitude  of  the  first  two  integrals  we  form 


\l 


the  ratio,  k(t) 


R(t)  = 


/  /  Kw.u  dA  dt 

o  A 


/  / DL(u) •  u  dA  dt 
o  A 
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Since  we  desire  order  of  magnitude  values  only  we  define  R 


Kw 


/  /  u  dA  dt 


R  = 


av  o  A 


DL(u)  •  u  dA  dt 
av  J  \ 


where  w  and  L(u)  represent  average  values .  Since  the  integrals  are  identi- 
av  av 


cal  we  can  substitute 

K  w 


R  = 


av 


K  w 


av 


DL(u)  DL  (w)  +  ;,L  (z) 

'  'av  '  av  v 


av 


Each  of  the  three  terms  in  the  last  equation  represent  a  pressure:  The 
individual  terms  become  more  meaningful  if  we  separate  the  motions  represented 
into  components  of  the  structure  motion. 

The  total  structure  motion  u  can  be  represented  as  the  sum  of  four  components 

u  =  z  +  z  +  w  +  w 
o  c  o  c 

where 

z  =  free-field  motion  at  structure  centroid,  constant  for 

entire  structure 

z  =  z  -  z 

c  o 

wq  =  departure  of  structure  centroid  from  free-field  position 

of  structure  centroid.  This  would  be  zero  for  a  perfect 

impedance  match  between  structure  and  soil 

w  =  w  -  w 

c  o 
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Now 


DL  (vo)  -  DL(zq)  -  0 

because  these  represent  undistorted  motion. 

Thus  ___  __ 

.  Kw  +  Kw 
r  «  0  c  av 

DL(w  )  +  DL(z  ) 

'  c  av  c  av 

Now  the  pressure  component  DL(w  )  is  equal  to  the  static  pressure  required 

c  av 

to  deform  the  structure  an  amount  wc  from  its  rest  configuration.  This  is  on  the 
order  of  a  small  fraction  of  the  peak  overpressure  that  the  structure  is  designed 
to  resist.  The  pressure  component  DL(zc)qv  however  represents  the  pressure  that 
would  be  required  to  uniformly  compress  the  structure  the  amount  that  the  pressure 
wave  compresses  the  free-field.  This  is  a  very  large  pressure,  particularly  for 
soft  soils  having  a  small  K.. 

Now,  the  departure  of  the  structure  centroid  from  the  free-field  centroid  is 
very  small  and  becomes  zero  shortly  after  engulfment.  Thus  the  term  Kwq  is  small 
in  comparison  with  the  others . 


Therefore,  within  the  accuracy  being  attempted  we  may  write 
^  Kw. 


R  - 


c  av 

DL(z  ) 
v  c  av 


Now,  w^  and  z are  of  the  same  order  of  magnitude,  probably  within  a  factor 
of  1.0  to  2. 

~  Kz 


R  - 


c  av 


DL  (z  ) 

c  av 
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Thus,  the  ratio  R,  is  of  the  same  order  of  magnitude  as  the  ratio  of  the 
pressure  required  to  reduce  the  volume  of  displaced  soil  to  the  pressure  required 
to  reduce  the  structure  volume  by  the  same  amount. 

In  general,  this  is  a  small  number. 


For  instance,  for  a  circular  cylinder  the  pressure  ratio  is  of  the  order 


■soil 

3cyl 


E*  p 

<!> 

c 


where 

Eg  =  modulus  of  elasticity  of  soil 

*  modulus  of  elasticity  of  cyMnder 

(R/h)  -  radius  to  thickness  ratio  of  cylinder. 


For  a  thin  concrete  silo  in  fairly  stiff  soil 


-soil  .  30000  v  pn 
Pcyl  4000000 

For  a  typical  existing  silo 


1 

W 


Psoil  _  20000  _  1 

p  ,  4000000  x  '  =  TJo 

ycyl 

This  lengthy  discussion  has  been  included  in  order  that  the  reader  might 
apt.'  eciate  that  the  conclusion  that  net  response  is  insensitive  to  K,  is  quite 
general,  and  does  not  hinge  on  numerical  values  chosen  for  a  few  specific 
structures . 


Therefore,  we  can  use  quite  approximate  values  of  K  without  seriously 
influencing  the  accuracy  of  results  obtained.  An  exception  to  this  statement 
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is  that  K  should  be  estimated  as  accurately  as  practical  when  computing  the  rigid 
body  motion  of  the  structure.  However,  we  are  fortunate  in  that  K  can  be  estimated 
with  beet  accuracy  for  this  mode  of  motion. 

It  is  not  generally  appreciated  that  the  formula 
P  »  p  c  v 

is  independent  of  any  assumed  stress-strain  law  if  the  c  value  used  is  the 
instantaneous  wave  velocity.  Of  course  the  wave  velocity  does  in  general  depend 
upon  the  stress  and  stress-strain  laws.  However,  for  both  the  linear  elastic 
and  linear  viscoelastic  stress-strain  laws  c  is  constant.  Therefore,  if  the 
correct  value  of  c  for  a  particular  soil  is  determined  the  only  errors  introduced 
by  the  velocity  terms  are  those  following  from  neglect  of  the  inclination  of  the 
pressure  vector  with  respect  to  the  velocity  vector.  This  angle  (7  of  Figure 
A. la)  is  generally  quite  small.  Further,  s^nce  the  value  of  c  in  the  equation 
is  adjusted  to  radiate  the  correct  total  amount  of  energy  during  one  cycle  (the 
rate  during  the  cycle  will  not  be  quite  correct,  in  general)  the  influence  un 
structure  motion  is  negligible. 

In  conclusion,  it  is  believed  that  Equations  (A.l8)  and  (A.19)  provide  nearly 
tie  best  accuracy  attainable  with  linear  equations.  A  slight  improvement  could  be 
made  without  loss  of  linearity  if  the  constant  K  were  replaced  with  a  linear 
operator  but  this  would  greatly  complicate  computation.  Any  further  improvement 
would  introduce  the  considerable  difficulties  associated  with  nonlinear  equations. 
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A. 3  Formal  Solution  of  the  Fundamental  Equations 


The  fundamental  equations  are 

•  M 

DL(u)  +  Ku  +  cP  u  +  mu  ■  p^  (A.l8) 

-  n  •  p  +  Kz  +  cpz^  (A.  19) 

If  we  have  any  set  of  orthogonal  functions  0^  satisfying  Equation  (A. 6) 

DL  (0^  )  -  m  0"  (A. 6) 

and  some  set  of  boundary  conditions  and,  further,  if  the  mass  per  unit  area  of 
the  structure  is  constant  then  the  fundamental  equations  can  be  solved  by  the 
following  set  of  equations 


1“°® 

U  -  2:  T,  (t)  $ 
i-1 

(A. 23) 

i  '  s:  {  V  ¥1  “ 

(A. 24) 

N.  *  |  m  f  .  0".  dA 

1  A  1  1 

(A. 25) 

T.+  —  T.  +  fh?  +  — )  T.  =  T . 
i  m  i  1  m  1  1 

v'A.26) 

If  the  mass  per  unit  area  were  not  constant  we  see  that  the  velocity  term  in 
Equation  (A. 26)  would  render  the  equation  meaningless  because  all  of  the  other 
terms  in  the  equation  are  functions  of  time  only. 

It  should  be  noted,  further,  that  although  we  have  chosen  to  carry  the  soil 
spring  constant  K  as  a  separate  factor  and  that  although  this  appears  in  Equation 
(A. 26)  as  the  ratio  —  there  is  no  necessity  that  K  be  taken  as  an  average  constant 


(A. 27) 


factor.  If  we  were  to  define  a  new  linear  operator 

✓N 

D  L  (u)  =  D  L(u)  +  K  u  , 

then  a  somewhat  different  set  of  orthogonal  functions  would  be  generated  by 

K 

the  equation  analogous  to  (A. 6)  and  the  factor  —  in  (A. 26)  would  vanish.  Thus, 
the  only  intrinsic  limitation  on  the  soil  resistance  factor  is  that  it  be  given 
by  some  linear  operator  operating  on  the  displacement  u. 

There  is  an  element  of  arbitrariness  in  the  solution  as  presented  to  this 
point  because,  as  has  been  stated,  the  orthogonal  functions  0  employed  can  be 
that  set  of  functions  generated  by  any  arbitrarily  assumed  set  of  boundary 
conditions.  One  would  suspect,  however,  that  of  all  of  the  possible  sets  that 
might  be  used  in  a  given  instance  some  particular  set  will  result  in  the  least 
amount  of  computational  effort  required.  This  aspect  of  the  theory,  though  of 
considerable  practical  Importance  to  the  engineers  who  must  carry  out  the  computa 
tions,  generally  is  glossed  over  in  treatises  on  orthogonal  functions.  Therefore 
1L  will  be  developed  in  detail  in  the  following  paragraphs  as  it  applies  to 
solution  of  equations  (A.l8)  and  (A. 19). 

Our  general  method  of  procedure  is  to  divide  the  total  motion  u  into  a  rigid 
body  component  u^  plus  a  deformational  component  u^ .  We  then  expand  both  and 
u^  in  terms  of  the  orthogonal  functions  0^  and  determine  boundary  conditions  for 
^resulting  in  the  simplest  overall  solution.  Thus, 


u  +  u , 
r  d 

(A. 28) 

i=  oo 

i=oo 

2  S.  0.+ 
i=l 

>:  t.  0. 

.  ,  ii 

1=1 

(A. 29) 
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Since  ur  is  a  rigid  body  motion 


DL  (uj.)  *  0 

Then,  substituting  (A. 24)  and  (A. 29)  into  (A.18) 
<s'i  +  V  +  %  (S,  +  V  +  §  (S1  +  V  +  «*  T, 


(A. 30) 


(A. 31) 


Suppose,  however,  that  it  were  permissible  to  set 

ur  -  Sr  ^  ,  (A. 32) 

that  is  the  entire  rigid  body  motion  could  be  represented  by  the  rth  mode  of 
the  system.  We  know,  for  instance,  that  this  is  true  for  a  complete  cylinder 
where  the  first  mode  represents  the  rigid  body  motion.  Then  for  all  values  of 
the  index  i,  other  than  r.  Equation  (A. 31)  would  become 

*i  +  *S  *t  +  +  Ti  !i  + r)  <*•»> 

In  view  of  Equation  (A. 30)  the  frequency  as  given  by  Equation  (A. 6)  would 


be  sex  u  fux  the  x'th  mode,  i.C., 

u^  =  0 


Thus  for  the  rth  mode  equation  (A. 31)  becomes 


(s  +  T  )  +  ££.  (3  +  T  )  •.  -  (S  +  T  )  -  T 
v  r  r  m  '  r  r'  m  '  r  r'  r 


(A. 34) 


(A. 35) 


where  we  see  that  for  the  rth  mode  the  distinction  between  S  and  T  is  arbitrary 

r  r 

and  T^  may  be  taken  equal  to  zero,  in  conformity  with  definitions  (A. 28)  and 
(A. 29)  which  restricted  u^  and  the  associated  time  functions  T^  to  deformational 
motion . 
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Thus,  if  a  rigid  body  mode  0r  satisfying  Equation  (A. 32)  can  be  shown  to 
exist  then  all  of  the  rigid  body  motion  would  result  from  a  single  solution  of 
Equation  (A. 26);  no  summations  of  the  type  indicated  by  Equation  (A. 29)  would  be 
necessary.  Inasmuch  as  the  rigid  body  motion  comprises  by  far  the  largest  per¬ 
centage  of  the  total  motion  for  small  stiff  structures,  such  as  solid  propellant 

/  —  7 

missile  silos,  a  criterion  for  establishing  the  rigid  body  mode  is  of  considerable 
practical  importance. 


If  a  rigid  body  mode  exists  it  will  satisfy  the  equation 


$  =  Constant  Vector 

r 


then  from  the  orthogonality  definition 


/  *  ?,  m  dA  =  0, 

A 


$  •  /  f  m  dA  =  0 

r  A  * 


Denote 


if/ .  =  J  0.  m  dA 
A 


(A.36) 


(A. 37) 


(A. 38) 


(A. 39) 


Thus,  the  vector  0 A  is  either  equal  to  zero  or  perpendicular  to  the  vector 
Since  in  general  the  infinite  set  of  vectors  will  not  all  have  the  same  direc¬ 


tion  cosines  we  can  conclude 


J  m  dA  =  0 


(A.40) 


Now  consider  the  structures  to  be  in  space,  not  acted  upon  by  any  exterior 
forces,  and  in  steady  state  vibration  in  the  0^  mode  about  the  mass  center.  The 
mass  center  under  these  conditions  is  motionless.  Then  the  structure  displacement 
is  given  by 


A -29 


V 


.Jr 


w  *  AJ  cb.  sin  0).  t 
i  i 


(A. 41) 


vhere 


A.^  is  arbitrary. 

The  velocity  and  structure  momentum  are  given  by 


V  “  ^i  Ai  cos  ^i  t 


m  *  mv  *  too.  A.  cos  cd.  t 
i  i  i  i 


where  m  is  the  momentum  of  an  element  of  the  structure. 


(A. 42) 


(A.l+3) 


Since  the  centroid  of  the  structure  is  motionless  the  total  momentum  of 


the  structure,  M,  must  be  zero  at  all  times.  Thus, 

M  =  /  in  dA  =  0 
A 

A  CD.  cos  CD  t  /  m  0  dA  =  0 
11  1  A  i 


Therefore 


/  m  0.  dA  *  0 
A  1 


(A. 1+1+) 


(A.  1+5) 


(A. 46) 


Equation  (A. 46)  is  identical  with  Equation  (A.40);  thus  the  criterion  for 
the  existence  of  a  rigid  body  mode  is  that  the  mode  should  be  computed  for 
structure  vibration  about  its  centroid  without  influence  from  external  forces 
or  reactions. 

For  self-closing  shells,  such  as  cylinders  or  toroids,  normal  modes  are 
customarily  computed  for  the  motionless  center-of-mass  condition.  However,  for 
open  shells,  such  as  domes  and  arches,  normal  modes  generally  are  computed  for  a 
position-fixed  boundary,  the  spring  line  of  the  shell,  that  is  not  in  the  plane 


A-30 


of  the  mass  centroid.  For  domes  the  normal  modes  computed  for  the  motionless 
mass  centroid  will  not  differ  greatly  from  those  appropriate  to  the  position 
fixed  spring  line;  however,  for  cylindrical  arches  having  no  tie  between  the 
foundations  the  difference  between  these  two  conditions  will  result  in  marked 
divergence  between  the  mode  geometries. 

Use  of  the  centroidal  normal  modes  imposes  one  complicating  factor  on 
computation  of  the  forcing  function.  Since  the  footing  bases  move  with  respect 
to  the  centroid  during  vibration  of  the  structure  the  generalized  forcing 
functions,  ,  must  include  the  refect  of  foundation  reactions,  as  defined 
by  Equation  (A.19)  does  not  include  this  term.  Essentially  the  intensity  of 
the  reaction  at  any  given  instant  of  time  would  be  equal  to  the  static  reaction 
necessary  to  balance  the  forces  distributed  over  the  surface  of  the  structure 
at  that  same  instance  of  time.  These  forces  include  the  three  components  on 
the  right  side  of  Equation  (A.19). 

As  a  practical  matter  when  open  structures,  domes,  arches,  toroidal  arches, 
etc.,  are  used,  generally  they  are  very  large  structures,  aerodynamic  systems  shelters, 
command  and  communication  centers,  etc.  For  the  large  structures  the  rigid  body 
motion  may  be  qui t  e  small  in  comparison  with  the  deformation  motion.  When  this 

> 

is  true  the  functions  will  be  considerably  smaller  than  the  functions  T  of 

i  Equation  (A. 3l) •  If  the  functions  are  small  enough  to  neglect  then  Equation 

| 

t  (A. 31)  becomes  identical  to  Equation  (A. 25)  and  if  position-fixed  boundary 

I  normal  modes  were  to  be  used  the  boundary  reaction  forces  would  contribute  no 

energy  to  the  system  and,  therefore  would  not  have  to  be  included  in  the  computa¬ 
tion  of  the  generalized  forcing  function.  Unfortunately,  no  general  imles  can 
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be  laid  down  for  determination  of  the  negligibility  of  with  respect  T^; 
however,  if  the  rigid  body  motion  were  to  be  computed  using  the  centroidal 
rigid  body  mode  (which  is  just  a  lateral  translation  of  the  rigid  structure) 
and  the  deformational  motion  computed  for  position-fixed  foundation  modes  the 
results  obtained  would  be  conservative,  i.e.,  the  deformational  motion  would  be 
somewhat  overestimated  though  the  degree  of  conservatism  might  be  difficult 
to  establish . 


An  important  characteristic  of  the  modal  motion  can  be  determined  from 
the  free  vibration  solution  of  Equations  (A. 26)  (and  equivalently  Equations 
(A. 33)  and  (A.35)) .  If  we  define  £2^ 


(A. 47) 

and  6  the  damping  coefficient  expressed  as  a  fraction  of  critical  damping 

26.^  =  ^  (A. 48) 

then  the  free  vibration  corresponding  to  Equation  (A. 26)  becomes 


T,  +  26ft.  T.  +Q2  T.  =  0  (A. 49) 

1  1  1  1  i  '  ' 

This  is  the  familial  equation  of  the  damped  oscillator.  Solving  for6i  we 

obtain 


6  =  £?- 
i  32  ^  m 


(A. 50) 


If  6^  is  greater  than  one,  the  mode  is  over-critically  damped.  For  the 
uver-critically  damped  modes  the  response  function  asymptotically  approaches 
the  forcing  function  .  If  the  latter  contains  no  spikes  -or  jumps  the  response 
function  generally  remains  less  than  the  forcing  function.  In  fact,  these 
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conclusions  are  substantially  correct  even  for  oscillatory  motions  occurring 
when  6  ^  is  less  than  one  but  greater  than  the  square  root  of  one-half . 

Many  methods  are  known  for  solution  of  equations  of  the  general  form  of 
Equation  (A. 26).  One  method  is  presented  in  Appendix  B  of  this  report;  however, 
it  should  be  appreciated  that  a  designer  has  his  choice  of  the  various  analytic, 
phase  plane,  computer,  or  other  methods  for  solution  of  the  equation. 

One  problem  in  the  solution  of  Equation  (A. 26)  is  evaluation  of  the 
forcing  function  r. .  A  detailed  general  method  of  approach  is  presented  in 
the  next  section. 

A. 4  Evaluation  of  Generalized  Forcing  Functions 

In  this  section  a  general  method  for  evaluation  of  the  forcing  functions 
r.  will  be  presented.  The  method  is  applicable  to  hand  numerical,  graphical, 
and  digital  computer  computation. 

In  the  development  that  follows  it  is  presumed  that  the  geometry  of  the 
normal  modes  0^  is  known.  Eor  most  structures  likely  to  be  used  the  modal 
components  are  available  in  the  literature. 

In  general  the  ground  shock  wave  will  not  be  a  plane  wave .  However,  by 
far  the  largest  portion  of  the  energy  delivered  to  the  structure  will  come  from 


the  air-induced  ground  shock  wave.  This  wave,  strictly  speaking,  is  not  plane, 
but  the  portion  of  the  wave  actually  impinging  against  the  structure  can  be 
taken  as  plane  with  negligible  error  even  for  structures  of  a  couple  of  hundred 
feet  span  and  waves  resulting  from  very  small  bombs .  It  should  be  noted  that 
the  air-induced  wave,  in  general,  will  have  at  least  six  components,  two  direct 
and  four  reflected  from  bedrock.  The  six  components  of  the  wave  are: 
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1.  The  incident  compression  wave 

2.  The  incident  shear  wave 

3.  The  compression  wave  resulting  from  reflection  of 
the  compression  wave. 

4.  The  shear  wave  resulting  from  reflection  of  the  compression  wave 

5.  The  compression  wave  resulting  from  reflection  of  the  shear  wave 

6.  The  shear  wave  resulting  from  reflection  of  the  shear  wave. 

These  six  components  are  shown  in  figure  A. 2.  It  should  be  noted  that  for 
a  substantially  level  bedrock  layer  that  the  unit  vectors  perpendicular  to  the 
advancing  wave  fronts  all  lie  in  the  same  vertical  plane,  but  that  the  inclination 
of  the  unit  vectors  with  respect  to  horizontal  varies  among  the  components. 

In  a  nonhomogeneous  soil  consisting  essentially  of  differing  horizontal 
strata  with  fairly  well  defined  interfaces  the  number  of  components  of  the  ground 
shock  wave,  each  still  substantially  plane,  arriving  at  the  structure  could  be 
considerably  increased.  In  addition  to  those,  at  the  position  of  the  Soi-ucturc, 
there  generally  will  be  a  large  number  of  randomly  oriented,  low  intensity  waves 
resulting  from  diffraction  of  the  air-induced  ground  shock  wave  around  boulders 
and  other  local  discontinuities  and  a  component  of  the  directly  transmitted 
shock  originating  at  the  crater  and  propagating  outward  through  the  soil  to  the 
structure.  Since  this  latter  component  must  travel  through  a  much  longer  soil 
path  than  does  the  air-induced  shock  to  reach  the  structure  the  number  of  random 
defractions  that  have  occ’irred  is  considerably  greater  for  the  directly  transmitted 
component  of  shock  than  for  the  air-induced.  In  the  100  to  1000  psi  range  the 
directly  transmitted  component  generally  appears  essentially  as  hash  consisting 
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Figure  A. 2  GROUND  WAVES  INDUCED  BY  BLAST  WAVE  AND  BEDROCK 
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of  spikes  of  acceleration  of  fairly  high  intensity  but  short  duration,  a  fraction 
of  a  millisecond  to  a  few  milliseconds . 

It  is  not  our  purpose  in  this  report  to  present  detailed  methods  for  deter¬ 
mining  the  free-field  ground  motion  at  the  structure  position.  This  information 
is  contained  in  a  number  of  lengthy  reports  devoted  exclusively  to  that  topic. 
Indeed,  this  is  one  of  the  major  areas  of  research  still  being  investigated. 
However,  in  Section  4  of  this  report  we  do  present  some  general  rules  for 
estimating  the  order  of  magnitude  of  importance  of  the  various  components  of 
shock  that  might  reach  the  structure. 

Whatever  methods  or  data  are  ed  to  establish  the  free-field  motion-time 
history  at  the  structure  position  these  data  will  include  pressure  in  three 
directions,  particle  velocity,  and  particle  displacement,  all  given  as  functions 
of  time.  Generally,  accelerations  associated  with  each  wave  component  also  are 
available. 

The  equation  for  the  generalized  forcing  function,  7\,  is 

n 

i.  f 

Ti  =  n7  Ja  pr  0.  dA  (A. 24) 

where 

n  •  p  +  kz  +  pcz  (A.  19) 

It  would  be  possible  of  course  to  expand  these  expressions  by  the  methods 
of  vector  analysis  and  obtain  formulas  for  their  evaluation  in  terms  of  the  scalar 
magnitudes  of  the  various  vectors  and  the  angles  between  these.  However,  a  much 
simpler  system  of  evaluation  can  be  developed  from  a  consideration  of  the  physical 
meanings  of  the  terms  and  introduction  of  an  artifical  device,  an  imaginary 
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rectangular  parallelepiped  circumscribed  about  the  structure.  Though  this 
artificiality  does  introduce  one  avenue  of  possible  confusion  (which  is  discussed 
subsequently)  it  materially  simplifies  the  computations  required  to  evaluate  the 


ri  integrals,  particularly  the  component  due  to  the  free  field  pressure  tensor, 

Now,  we  note  that  in  integral  (A. 24)  the  product  of  the  pressure  vector  p^  j 

and  the  differential  area,  dA,  is  a  differential  force.  The  dot  product  of  this  * 

force  with  the  displacement,  then  is  the  differential  energy  delivered  to  the 
structure  by  the  force  p^  dA  moving  through  a  distance,  0^.  Thus  the  integral 
in  Equation  (A- 24)  is  the  total  energy  that  could  be  delivered  to  the  structure  i 

by  the  force  field  p^  if  the  latter  remained  constant  while  the  structure  deformed 
into  the  displacement  configuration,  0^. 

Consider  any  shell  structure.  About  the  structure  describe  an  imaginary 
rectangular  parallelepiped.  It  is  not  essential  that  the  parallelepiped  hive  any 
pciru  i  c  ulai’  oiuemauion^  ncvevcr,  computalions  will  bt:  simplest.  ii  iuur  faces  of 
the  parallelepiped  are  vertical,  two  of  these  parallel  to  the  direction  of  the 
advancing  ground  shock  waves.  Figure  A. 3  illurtriboo  such  a  parallelepiped 
circumscribed  about  a  vertical  cylinder.  Three  coordinate  directions  £,  r),  and  £ 
and  three  unit  vectors  f,  7],  and  £  are  defined  ir.  the  figure.  In  Section  AA  of 
Figure  A. 3  a  differential  element  of  area,  dA,  has  been  emphasized.  A  pressure 
vector  p^  and  a  displacement  vector  0  act  on  the  area  dA.  For  clarity  these  are 
shown  in  the  f,  rj  plane  though  this  is  not  essential. 

Now  if  both  the  force  vector,  p^,  and  the  displacement  vector,  0,  are 
resolved  into  their  components  in  the  £,  r),  and  £  directions  the  results  are 
as  shown  in  the  figure  (only  the  £  and  77  components  are  shown).  Thus  the 
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energy  delivered  to  the  element  dA  is  equal  to  the  sums  of  the  energies  delivered 
by  the  components  of  the  pressure  vector  perpendicular  to  the  faces  of  the 
parallelepiped  acting  through  the  projections  of  the  displacement  volume  on  the 
same  faces . 

The  vectors  kz  and  pcz  present  no  problems  in  resolution  into  their 
components  on  the  faces  of  the  parallelepiped. 

The  tensor  p  presents  a  somewhat  different  problem.  Though  the  vector 
n  •  p  is  very  simple  to  represent  symbolically  its  actual  numerical  computation 
can  become  quite  tedious.  It  turns  out  that  this  computation  need  not  be  carried 
out.  Figure  A. 4  shows  a  wave  front  and  also  a  differential  element  of  soil  behind 
the  wave  front.  Tk  components  of  pressure  acting  on  the  differential  element  of 
soil  are  shown.  Generally  the  components  shown  as  pfe  and  pc  of  Figure  A. 4  are 
equal  but  this  is  not  necessarily  true.  The  Mohr's  cycle  for  the  pressure  in  the 
7 1,  f  plane  also  is  given  on  Figure  A. 4,  and  the  component  of  pressure  p  appearing 
on  the  vertical  face  uf  the  parallelepiped  ABCD  is  shown  as  p.  The  horizontal 

component  of  .this  in  the  r\  direction  is  shown  as  p  •  rj.  A  factor  X  is  defined  on 

\ 

v.  J  «—  -- 

the  Figure  where  XpQ  is  equal  to  p  •  7).  Now  since  behind  the  wave  front,  pQ,  p^, 

and  pc  bear  a  constant  relationship  to  each  other  the  factor  X  is  a  constant 

applicable  to  planes  ABCD  and  EFGH.  It  should  be  noted  that  the  pressure  p  •  T) 

appearing  on  plane  in  ABCD,  is  the  same  in  numerical  magnitude  as  that  appearing 

on  the  plane  EFGH,  other  conditions  being  equal.  Thus  the  pressure  at  any  point 

on  plane  ABCD  or  EFGH  is  just  the  wave  pressure  p  multiplied  by  a  constant  factor 

8 

X^  which  is  easily  determinable  in  any  specific  instance.  Further,  if  the  ground 
pressure  wave  decays  as  a  function  of  time  the  pressure  on  a  vertical  plane 


A- 39 


Figure  A.k  PRESSURES  INDUCED  BY  PLANE  COMPRESSION  WAVE 
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perpendicular  to  the  T)  axis  at  anv  given  instant  is  equal  to  the  pressure  p 


a 


existing  at  that  instant  multiplied  by  X  .  The  X  factor  is  neither  position 


V 


nor  time  dependent  for  any  given  plane  wave.  Factors  X4  and  X^  can  similarly 
be  defined  for  the  other  four  faces  of  the  parallelepiped. 


Figure  A, 5  shows  a  similar  construction  for  a  shear  wave.  It  should  be  noted 
that  the  shear  wave  has  no  component  on  planes  ABFE  and  DCGH.  \  factors  for 
these  waves  can  be  defined  similarly  to  those  for  the  compression  waves. 


The  time  factor  involved  in  the  wave  motion  has  not  yet  been  considered. 
Figure  A. 6  shows  a  section  of  the  cylinder  of  Figure  A. 3  uoken  parallel  to  the 
rjj  £  sixes.  A  compression  wave  is  shown  during  the  period  of  engulfment.  Nov, 
it  should  be  carefully  observed  that  although  the  computations  are  carried,  out 
by  applying  pressure  components  to  the  surfaces  of  the  parallelepiped  the  varia¬ 
tion  of  pressure  component  versus  time  on  each  element  of  the  modal  displacement 
component  on  the  parallelepiped  is  identical  with  the  pressure  versus  time  at  the 
element  of  structure  surface  projected  on  the  faces  of  the  parallelepiped,  e.g., 
pressure  on  element  da.  at  time  t  is  equal  to  the  component  of  pressure  on 

A 

element  dA^  at  time  t  projected  perpendicular  to  plane  ABCD. 


Thus  for  the  instant  depicted  in  Figure  A. 6  if  the  pressure  versus  time 
curve  is  represented  p(t).  where  time  is  taken  as  zero  at  the  wave  front,  and 

A 

if  we  denote  the  component  of  pressure  at  element  da^  by  p^  we  have 


A 


Projection  P(0) 
=  Projection  P(0) 


(A. 51) 
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Figure  A.  5  PRESSURES  INDUCED  BY  PLANE  SHEAR  WAVE 
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and 


*  Projection  P(^-) .  (A. 52) 

The  distance,  s,  is  measured  perpendicular  to  the  plane  of  the  wave  front. 

In  summary,  the  integral  (A. 24)  can  be  evaluated  as  follows: 

1.  For  the  ith  mode  compute  the  norm  Ih  . 

2.  Circumscribe  an  imaginary  parallelepiped  about  the  structure. 

3.  Impose  the  displacement  of  a  normal  mode  upon  the  structure 
and  translate  components  of  this  deformation  onto  the  six  faces 
of  the  parallelepiped.  These,  then,  will  appear  as  relief  maps. 

4.  Divide  the  faces  of  the  parallelepiped  into  a  grid  of  small 
elements  (acceptable  accuracy  generally  can  be  obtained  with  a 
course  grid  of  four  or  five  elements  in  each  direction) . 

5.  Determine  the  time  lag  for  each  element. 

6.  Determine  the  product  01  pressure  components  and  volume  for  each 
element  and  each  of  three  components  of  pressure. 

7.  Gum  the  mines  for  all  elements  at  specific  time  instants  and 
divide  by  the  norm  N^. 

8.  Plot  the  results  as  a  function  of  time.  This  is  the  function 
T 

i* 

The  actual  computation  outlined  in  the  previous  eight  steps  becomes  a  matter 
of  bookkeeping  and  for  different  waves  and  structure  characteristics  differing 
techniques  will  be  most  expeditious.  Several  general  characteristics  that  can  be 
used  to  simplify  the  overall  procedure  will  be  enumerated. 
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1. 


For  plane  waves  the  factor  p cz  is  equal  in  numerical  magnitude  to 
Pa  of  the  pressure  tensor;  however,  whereas  the  latter  is  compression 
on  both  the  windward  and  leeward  faces  of  the  structure,  i.e.,  plane 
ABCD  and  plane  EFGH,  the  pcz  term  is  compression  on  the  windward  face 
and  tension  on  the  leeward  face.  The  direction  of  does  not  change 
as  it  passes  the  structure.  Also  it  should  be  noted  that  the  factor 
pcz  contributes  no  component  of  pressure  to  planes  ABFE  and  DCGH  if 
they  have  been  oriented  as  suggested. 

The  factor  Kz  produces  a  compression  stress  on  planes  ABCD  and  ABHD. 

It  produces  a  tension  stress  on  planes  EFGH  and  BFGC;  it  produces 
no  stress  on  planes  ABFE  and  DiCGH. 

3-  If  the  structure  is  subjected  to  more  than  a  single  wave  generally 
the  simplest  procedure  is  to  compute  the  effects  of  each  wave 
separately  and  then  sum  the  results  in  the  last  step. 

4.  If  mor  than  a  single  mode  is  to  be  investigated  all  factors  for  the 
various  modes,  except  the  displacement  factors  (<5\  •  r;  da,  for  instance) 
are  identical  for  the  several  modes. 

5.  For  rigid  body  motion  the  "relief  map"  on  planes  ABCD,  EFGH,  AEHD, 

and  DFGC  represent  just  displacements  of  the  planes  along  their  normals. 
The  planes  ABFE  and  DCGH,  contribute  zero  to  the  rigid  body  motion. 

6.  Each  component  of  the  rigid  body  motion  is  orthogonal  to  all  of  the 
other  modes  of  the  system  and,  therefore,  may  be  handled  separately. 

The  reader  may  well  question  whether  the  detailed  calculations  outlined  above 
are  justified  in  view  of  the  inherent  inaccuracy  of  input  data.  Generally  they 
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will  not  be  justifiable.  However ,  the  general  characterisitics  of  the  inodes  and 
essential  geometric  relations  presented  can  be  used  in  specific  instances  to 
delineate  which  of  the  input  wave  components  are  significant,  which  can  be 
neglected  (generally  all  three  shear  components  and  the  shear  reflection  induced 
compression  components),  and  to  establish  the  cardinal  characteristics  of  the 
forcing  functions  t  such  as  rise  time,  peak  value,  and  decay  rate.  Further, 
rough  estimates  of  the  peak  amplitudes  and  rise  times  of  the  higher  order  modes 
can  be  determined  generally  with  the  end  result  that  the  higher  order  modes  can 
be  neglected. 

A. 5  Coupled  Modes 

In  reality  the  modes  of  the  system  that  we  have  approximated  by  Equations 
(A. 18)  and  (A. 19)  are  coupled,  the  coupling  being  introduced  by  the  two  first 
time  derivative  terms  of  Equation  (A.l6). 

In  this  section  we  will  develop  the  theory  of  coupled  moles  for  this  general 
system.  The  results  indicate  that  the  simplifications  made  are  acceptable  from  a 
practical  standpoint.  Further,  the  theory  provides  the  necessary  corrections  for 
the  rare  instances  where  the  modal  coupling  must  be  considered. 

If  we  make  the  simplification  that  the  nonlinear  resistance  term  kF(w)  can 
be  replaced  by  the  term  K  Equation  (A.l6)  becomes 

•  •• 

KL(u)  +  K(u)  +  c  p  n  u  •  n  +  c  p  h  u  •  1  +  mu 

_  =  _  (A. 53) 

=  n  •  p  +  kz  +  c-jp  nz  •  n  +  c 2p  i  z  •  i 
•  • 

We  now  define  two  vectors  ih  and  ib 

ru  Yz 

where  ^  c _  c 2 •_  _ 

"  —  nu-n  +  —  £u-£-u  (A.  54) 
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Referring  to  Figure  A.l  we  see  that  pc  ip  and  pc  ip  are  the  vector  differences 

11  z 

between  Equation  (A.l8)  and  (A. 19)  and  the  accurate  Equation  (A. 53)' 

Making  this  substitution  and  shifting  all  correction  terms  to  the  right  side 

of  the  equal  sign,  Equation  (A. 53)  becomes 

•  •  •  • 

DL(u)  +  Ku  +  cpu  +  mu  *  n  •  p  +  Kz  +  cpz 

•  • 

+  -  c^  (A. 56) 


DL(u)  +  Ku  +  qou  +  mu  -  ' 


(A. 57) 


Since  ip  is  defined  in  terms  of  the  input  parameters  it  is  immediately 
z 

computable  so  we  can  simplify  the  appearance  of  Equation  (A. 57)  slightly  by 
defining 

Pi'  *  Px  +  *p\  (A. 58) 

Then 


DL(u)  +  K(u)+  cpu  +  mu  =  p^'  -  c pip  (A. 59) 

Equation  (A. 59)  is  of  the  same  form  as  Equation  (A.l8)  with  the  term  added. 
Define  '  *  ir  /  P, '  •  &  dA  (A.60) 

Then 
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(A. 61) 


Substitute  Equations  (A. 6),  (A. 23),  (A. 54)  and  (A.6l)  into  Equation  (A. 59)* 
Divide  through  by  m.  The  result  is 

i=l 


l”oo 
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(A. 62) 


Equations  (A. 66)  can  be  solved  by  successive  approximations.  At  least  five 
conclusions  can  be  drawn  f :  om  Equations  (A. 65)  and  (A. 66). 

1.  The  coupling  factors  remove  energy  from  each  mode  and  deliver  it 
to  all  other  modes. 

2.  The  coupling  between  modes  is  a  function  of  the  lateral  (as  opposed 
to  normal)  motions  of  the  inodes.  (The  amplitudes  of  these  are 
invariably  small  fractions  of  the  normal  components  of  motion.) 

3-  The  numerical  values  of  the  factors  within  the  brackets  are  small 
and  decrease  rapidly  with  increasing  modal  index.  The  values  are 
the  largest  for  right  circular  cylinders.  Even  for  this  case  the 
coupling  factor  between  the  first  and  second  deformational  modes 
is  less  than  3$. 

For  tne  T^  monotonically  approaching  zero  the  coupling  of  energy  to 
other  modes  is  very  small  (less  than  1$  total  for  the  first  10  modes). 

5.  If  damping  (£p)  is  very  small  and  if  the  T^  are  approximately  cyclic, 
over  a  period  of  time  considerably  energy  can  be  interchanged  between 
modes . 

We  conclude  therefore  that  the  simplification  introduced  in  the  previous 
subsections  will  not  materially  affect  computed  results  for  the  heavily  damped 
underground  structures . 

A  development  similar  to  the  above  can  be  carried  out  when  the  mass  of  the 
structure  varies  from  point  to  point. 


A.  6  Natural  Frequency 

In  order  to  use  the  material  presented  in  this  chapter  it  is  necessary  to 
know  or  determine  the  geometries  of  the  normal  modes  ^  and  the  associated 
natural  frequencies. 

In  theory  these  can  be  obtained  by  solution  of  Equation  (A. 6)  together 
with  the  boundary  equations 

DL(0i)  *  mSiJj2  0i  (A. 6) 

Bn(0i)  =0  (n  -  1,2..) 

These  equations  are  deceptively  simple.  Actually  the  operator  L(  )  may  be 
a  complex  differential  expression  and  in  fact,  since  actual  numerical  work  must 
generally  be  carried  out  for  the  vector  functions  in  terms  of  their  components 
Equation  (A. 6)  generally  must  be  replaced  by  two  or  three  scalar  equations  that 
are  interdependent.  The  boundary  equations  are  similarly  multiplied  in  number. 

We  are  fortunate  in  that  most  of  the  structure  types  of  interest  have 
already  been  subjected  to  such  analysis  and  as  a  result  the  frequencies  and 
modal  geometries  (generally,  the  components  rf  the  modes)  are  known. 

If  a  structure  having  geometry  significantly  different  tnan  those  for 
which  modal  data  exists  must  be  analyzed  by  the  rormal  mode  method  the  modal 
geometries  and  frequencies  can  be  determined  by  approximate  methods.  The 
determination  need  not  be  highly  accurate. 
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APPENDIX  B 

A  NUMERICAL  METHOD  OF  SOLUTION 
OF  THE'  RESPONSE  EQUATION 


A  method  for  solving  equations  of  the  form 

T  26  '  ~ 

—  +  7?-  T  +  T  =  r 

Cl2  “  *  " 


(B.l) 


is  desired  (Note  T  is  equivalent  to  rfo  of  Section  5  and  Appendix  A.  In 
this  appendix  it  is  presumed  to  be  a  known  function  of  time. 

We  have  the  initial  conditions  that 


T  (0)  -  0 

T  (0)  =  0 

A 

7  (0)  -  0 


(B.2) 


Integrating  once 


m  2ft  ^  t  ~ 

T  +  /  T  dt  =  /  rdt 
17  00 


(B.3) 


Now  suppose  that  the  time  axis  were  divided  into  a  number  of  equal  incre¬ 
ments  At.  We  order  the  increments  from  zero  to  m  with  zero  at  the  origin.  Also 

we  denote  the  values  of  T  and  t  at  the  end  of  the  nth  increment  by  T  and  r  . 

n  n 

With  these  definitions  we  can  approximate  Equation  (B-3)  by 


T  , ,  -  T  . 
n+1  n-1 

2fi2  At 


ir1  +  At 


i=n  I  r.  i=n 

2  T  -  At  JP+  2  T 

i-1  11  l_  i=1  1 


(B.4) 


If  we  denote 


flAt  «As 

we  have  for  all  n  greater  than  zero 

2  A 

Tn+1  =  (ri  ‘  Tl>  +  2^s)2 


-  46  As  I  +  T  . 

n  n-i 


For  n  *  0  we  obtain  for  T. 


i=n 

^  i*n 

2 

Ti  "  X 

i-l 

l-l 

T 


i  (As  y 


fis)*  +  46(As)  +  2 


(B.5) 


(B.6) 


(B.7) 


Equation  (B.6)  can  be  simplified  for  computation 
Denote 


n  x  t  -  T 
n+1  xn+l  n 


then 


<W  [Tn-In) 

-  («As  +  l)(Tn-Tn  l) 

+  <T„-i  -  VP 

This  can  be  rewritten 


A  A 


°n+l  ’2(is)  [Tn-Tnl 


-  (46As+l)  C  +  C  . 
v  '  n  n-1 


(B.8) 


(B.9) 


(B.10) 
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In  this  scheme 


^  P 

T,(^r 

Q  B  B  _ i _ _ 

1  1  (As)2  +  46  (As)  +  3 

c2  «  |  (As)2  (ri  -  cx)  -  (tfAe+l)  CL 


(B.ll) 

(B.12) 


Table  B.l  shows  one  setup  for  the  computations. 


TJlBLE 


A  -  2V*B' 


B  =  (U6  As+l) 
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